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Abstract 


Observational studies were performed using historical rocketsonde data to 
investigate long-term temperature trends, solar-cycle variations, and interactions between 
tropical and extratropical latitudes in the middle atmosphere. Evidence from tropical, 
subtropical, and midlatitude North American rocketsonde stations indicated a consistent 
downward trend over 25 years, with a solar cycle component superposed. The trend is 
about -1.4 to-2.0 K per decade and the amplitude of the decadal oscillation is about 1.1 K. 
Prior to trend derivation it was necessary for us to correct temperatures for aerodynamic 
heating in the early years. The empirically derived correction profile agrees well with a 
theoretical profile of Krumins and Lyons. A study was also performed of the correlation 
between equatorial winds and north polar temperatures in winter, showing that the entire 
stratospheric wind profile near the equator - including the quasi-biennial oscillation 
(QBO) and stratopause semiannual oscillation (SAO) - is important to the extratropical 
flow, not merely the QBO component as previously thought. A strong correlation was 
discovered between winter polar temperatures and equatorial winds in the upper 
stratosphere during the preceding September, suggesting a role for the second cycle of the 
SAO. 


Further investigation of the Arctic Oscillation demonstrated a significant 
downward propagation of this mode from the mesosphere to the troposphere in winter. 
The statistical significance of this relationship was established. Finally, a contemporary 
and comprehensive review of the quasi-biennial oscillation was prepared for publication, 
reporting results from a 1998 workshop on the QBO. 
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Middle atmosphere cooling trend in historical 
rocketsonde data 

Timothy J. Dunkerton, Donald P. Delisi, and Mark P. Baldwin 

Northwest Research Associates, Bellevue, Washington 


Abstract. Data from the historical rocketsonde net- 
work demonstrate that significant cooling of the up- 
per stratosphere and lower mesosphere (~30-60 km) 
occurred in northern midlatitudes of the western hemi- 
sphere and in the tropics during 1962-1991. The down- 
ward trend of temperature averaged over this layer was 
about -1.7 K/decade and temperatures were apparently 
modulated by the solar cycle with amplitude ~1 .1 K. 
The trend was a function of height and somewhat larger 
in the lower mesosphere relative to the middle and up- 
per stratosphere. 


Introduction 

Cooling of the middle atmosphere is expected in re- 
sponse to the increasing concentration of greenhouse 
gases of anthropogenic origin. This signature is the 
largest and theoretically the most obvious component 
of climate change. Detection of a temperature trend 
above 10 hPa, however, is compromised by a number 
of problems including poor spatial and temporal cover- 
age, instrument calibration, and inhomogeneities within 
each dataset. Nevertheless, there is general agreement 
that significant cooling has occurred over the last 2- 
3 decades. The observational evidence includes mea- 
surements from the Stratospheric Sounding Unit (SSU) 
[Nash and Edge , 1989], the historical rocketsonde net- 
work [Angell, 1987, 1991; Golitsyn et a/., 1996, and 
references therein], and Rayleigh lidar [Hauehecorne et 
a/., 1991; Keckhut et at, 1995]. There is also indirect 
evidence of cooling in sodium lidar data suggesting a 
hydrostatic contraction of the atmosphere below ~95 
km [Clemesha et ai } 1992], and a possible increase of 
noctilucent cloud sightings in the summer mesosphere 
[ Thomas , 1991]. 

Identification of a climate trend distinct from natu- 
ral interdecadal variability will require much more data 
than are currently available. Nevertheless, it is likely 
that observations in the modern era (1950s and there- 
after) will play a pivotal role in the detection of human 
influence, having been acquired in a time period when 
such an effect might be first observable. Data from the 
historical rocketsonde network provide the longest con- 
tinuous evidence of cooling in the upper stratosphere 
and lower mesosphere ( ~ 1 962—199 1 ) . In order to make 
optimum use of these data it has been necessary for 
us to assemble the dataset from various sources and to 
carefully analyze the data for possible biases and inho- 
mogeneities. 

Copyright 1998 by the American Geophysical Union. 
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Data Analysis 

Data from the historical rocketsonde network (World 
Data Center ‘A’) including North American stations 
and a few tropical sites for the period 1969-1991 were 
obtained from the National Center for Atmosphere Re- 
search, supplemented by additional data from the Cli- 
mate Prediction Center (Mel Gelman, personal commu- 
nication) and National Climatic Data Center. These 
data are comprised of individual rocket soundings ac- 
companied by rawinsonde data from the same stations. 
For the period before 1969 we obtained pre-processed 
data formerly utilized by scientists of the meteorology 
group at Control Data Corporation (David Venne, per- 
sonal communication). These data were originally ac- 
quired from rocket soundings but had been processed 


a 



K 


b 



Figure 1. (a) Average of empirical temperature cor- 
rection profiles (black line) plus or minus one standard 
deviation (shading). The standard correction is also 
shown (gray line), (b) Average profile of temperature 
trend, after the average correction was applied to sta- 
tions individually, plus or minus one standard deviation 
(shading). 
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Figure 2. Time-height cross section of monthly mean deseasonalized temperature at Kwajalein. (a) Uncorrected 
and (b) corrected data are shown with the same color scale. At Kwajalein, the transition to corrected data 
occurred in April 1975. Units: K. 


onto a regular vertical grid at 2 km intervals. The pre- 
processed data were obtained by us through 1982, and 
were found to agree with the original data in their pe- 
riod of overlap. Of the 22 stations included in this rock- 
etsonde network, 17 provide data for more than 10 years 
and 9 provide data for more than 20 years. Data at most 
stations were analyzed but in this letter we focus atten- 
tion on six of these lying either in western hemisphere 
midlatitudes or tropics: 

Ascension (14.8°W, 8.6°S), 

Kwajalein (167. 4°E, 8.4°N), 

Barking Sands (159.9°W, 22.0°N), 

Cape Kennedy (80.7°W, 28.3°N), 

Point Mugu (119. 5°W, 34.1°N), 
and 

Wallops Island (75.7°W, 37.5°N). 

Their location outside the winter polar vortex is such 
that interannual variability is relatively small. 

Angell [1987, 1991] compared rawinsonde and rock- 
etsonde temperatures at 26-35 km and concluded that, 
because of a sharp apparent cooling of rockets with re- 


spect to rawinsondes in the early 1970s, the rocket trend 
at that time was spurious. This is consistent with a sug- 
gestion by Johnson and Gelman [1985], who examined 
data from the North American network in the band 25- 
55° N. Our analysis shows that individual stations did 
not exhibit this cooling at the same time, and the spu- 
rious ‘trend’ really amounts to a near-discontinuity in 
measured temperature. We examined header informa- 
tion for soundings at tropical and midlatitude stations 
and determined that the spurious jump was due primar- 
ily to a change from uncorrected to corrected tempera- 
tures (to account for aerodynamic heating, etc.) rather 
than an instrumental change (e.g., Arcasonde to Data- 
sonde). In this letter we therefore display temperature 
data using all instrument types. Included in the dataset 
are a few measurements obtained from falling spheres. 
The sphere data are colder, generally falling outside the 
range of observed layer-mean temperatures, but are rel- 
atively few in number and do not affect the results sig- 
nificantly. 

In order to utilize uncorrected data before the early 
1970s it is necessary to subtract a positive tempera- 
ture correction profile. The standard correction is doc- 
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umented in Krumins and Lyons [1972]. For compari- 
son, we chose instead to derive empirical correction pro- 
files for each station in the following manner. Biennial- 
mean deseasonalized temperature profiles were obtained 
immediately before and after the switch from uncor- 
rected to corrected data at each station (the biennial 
interval chosen to minimize effects of the quasi- biennial 
oscillation). We then assumed that the two biennial- 
mean profiles should have been equal apart from a lin- 
ear trend. The empirical correction profile and final 
trend estimate were obtained at each level by iterating 
this procedure a few times until convergence. The aver- 
age correction profile is shown in Figure la, and agrees 
well with the standard correction, except at lower lev- 
els where the empirical correction is about 1 K. Both 
corrections increase with altitude near the stratopause 
in approximately exponential fashion; the scale height 
is about 7.9 km for the standard correction, close to a 
density scale height. For the purpose of plotting data 
above 60 km, where our derivation of empirical cor- 
rection would be unreliable, the correction profile was 
extrapolated using a constant scale height. 

Application of the average correction to uncorrected 
data reduces the net trend over the entire record but 
does not eliminate either the trend or an apparent solar 
cycle influence as demonstrated below. The average 
trend of the six stations as a function of height is shown 
in Figure lb: approximately -2 K/decade in the upper 
stratosphere, increasing to in excess of -4 K/decade in 
the lower mesosphere. 


Results 

Monthly mean uncorrected and corrected tempera- 
tures at Kwajalein are shown in Figures 2a and 2b, re- 
spectively. These figures display rawinsonde data below 
3.7 scale heights (~25 hPa), rocketsonde data above 4.7 
scale heights (~9 hPa), and a linear blend of the two 
in their region of overlap. Prominent features include 
descending phases of the quasi-biennial oscillation, an 
apparent decadal variation coinciding with the solar cy- 
cle (with maximum temperature near solar maximum 
and vice versa), and a cooling trend throughout the 
record. Note the apparent positive correlation with the 
solar cycle: e.g., solar maxima occurred around 1969, 
1981, and 1990; solar minima occurred around 1976 and 
1986. A transient effect of the El Chichon eruption in 
1982 was prominent at Kwajalein in the middle half of 
this year with warming in the lower stratosphere and 
cooling aloft as documented by Dunkerton and Dehsi 
[1991]. A cooling trend was present at most levels of 
the middle atmosphere, whereas the apparent solar in- 
fluence was seen primarily above 35 km. 

Figure 3a compares rawinsonde and corrected rock- 
etsonde temperatures averaged vertically in their re- 
gion of overlap, at five stations. At each location the 
two datasets agree reasonably well, especially the high- 
frequency variability; however, there was a brief period 
near 1970-71 when rocket temperatures were anoma- 
lously high by about 2 K, and another interval after 
1984 when rockets were slightly cool by about 1-2 K. 
These discrepancies affect the trend estimates, which 
range from -0.19 to -0.73 K/decade for rawinsondes, 
and from -0.81 to -1.33 K/decade for rocketsondes. The 
downward trend is real, but possibly exaggerated in 
rocketsonde data. At higher levels, there is good agree- 



Figure 3. (a) Layer-mean temperature at 3. 7-4. 7 scale 
heights for five stations, derived from corrected rock- 
etsonde data (black line) and rawinsonde data (gray 
line), (b) Layer- mean rocketsonde temperature at 4.1- 
7.9 scale heights. A linear least-squares fit to each series 
is also shown (light solid). 


ment between stations concerning both the downward 
trend and apparent solar influence, as shown in Figure 
3b. Here, temperatures are averaged vertically over the 
layer 4. 1-7.9 scale heights (approximately 28-56 km) . 
Over the last 25 years, downward trends of magnitude 
-1.38 to -2.01 K/decade are obtained after correction. 
Amplitude of the decadal oscillation is ~1.1 K. 

Figure 4 shows the average of the individual station 
data shown in Figure 3b. Data are plotted only when 
three or more stations have data in a given month. 
The dashed straight line is the linear least-squares fit 
through two solar cycles, starting at the solar minimum 
in 1965 and ending at the solar minimum in 1986. The 
slope of this line is -1.76 K/decade. The solid straight 
line is the linear least-squares fit through two solar cy- 
cles, starting at the solar maximum in 1969 and ending 
at the solar maximum in 1990. The slope of this line 
is -1.69 K/decade. The 10.7 cm solar flux is shown in 
the bottom part of this figure. Several interesting fea- 
tures are notable. For example, the data show a posi- 
tive correlation with the solar flux, indicating an appar- 
ent, solar cycle in the height-averaged temperature data. 
It should be recognized that the apparent solar signal 
is small and is best seen when the data are averaged 
vertically and data from several stations are blended 
together; the significance of this result should not be 
overstated. The data also show a reasonably constant 
trend over nearly three decades (in the sense that sim- 
ilar trends are obtained using subintervals within this 
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Figure 4. Average of layer-mean temperatures at 4.1- 
7.9 scale heights for the five stations of Figure 3. Trend- 
lines between solar maxima (solid) and solar minima 
(dashed) are superposed. 

record stretching over two solar cycles). This result 
appears to validate the temperature correction method 
even though the break points occurred at different times 
depending on location. For comparison, using only one 
solar cycle [e.g., OTP , 1988], the trends for solar min- 
imum to solar minimum are -2.53 and -1.92 K/decade, 
and the trends for solar maximum to solar maximum 
are -1.93 and -1.28 K/decade. 

Conclusion 

Data from the historical rocketsonde network demon- 
strate that significant cooling of the upper stratosphere 
and lower mesosphere (~30-60 km) occurred in north- 
ern midlatitudes of the western hemisphere and in the 
tropics during 1962-1991. The downward trend of tem- 
perature data averaged over this layer was about -1.7 
K/decade and temperatures were apparently modulated 
by the solar cycle with amplitude^l.l K. The trend was 
a function of height and somewhat larger in the lower 
mesosphere relative to the middle and upper strato- 
sphere. 

Uncorrected temperature data prior to the early 1970s 
require the user to correct for aerodynamic heating and 
other heat transfer effects; the date before which correc- 
tion is required varies from station to station. Correc- 
tion profiles derived empirically are similar to the stan- 
dard correction, increasing with altitude approximately 
in exponential fashion with a scale height close to that 
of atmospheric density. The magnitude of derived tem- 
perature trend depends on the application of a suitable 
correction to early data; this is the most salient issue 
for trend studies using these data. The overall behavior 
of temperature time series does not depend significantly 
on thermistor type (e.g., Arcasonde or Datasonde) since 
instrumental changes were relatively few and most vari- 
ability was observed in periods of time monitored by the 
same instrument. Other factors not discussed here such 
as measurement time of day and agreement with adja- 
cent rawinsonde might alter the trend estimate, but our 
investigation (to be reported separately) indicates that 
these are relatively minor effects. 

Putting our results in a larger context, although evi- 
dence of similar cooling has been found in other dat asets 
[e.g., Golitsyn et a/., 1996], it is likely that tempera- 
ture trends associated with climate change are not glob- 
ally or hemispherically uniform but contain a dynami- 
cal component associated with planetary waves and in- 


duced mean meridional circulation. Further analysis of 
existing ground-based and satellite datasets and synthe- 
sis with newer methods of observation (e.g., lidar) will 
help to clarify the nature of the climate change signal 
as affected by atmospheric dynamics. 
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Summary 

Equatorial winds in the stratosphere are known to influence the frequency of stratospheric sudden 
warmings. Sudden warmings, in turn, influence the Earth’s climate both through their direct influence on 
polar temperatures and through the temperature dependence of ozone depletion in the lower stratosphere. 
The conventional (Holton-Tan) explanation for the equatorial influence on stratospheric warmings is 
in terms of the equatorial winds in the lower stratosphere (~20-30 km) acting as a wave-guide for 
mid-latitude planetary wave propagation. This study employs stratospheric temperature analyses and 
equatorial rocketsonde wind data extending to 58 km to diagnose the relationship between the Northern 
Hemisphere polar temperatures and equatorial zonal winds at all height levels in the stratosphere. In 
addition to the recognised Holton-Tan relationship linking the polar temperatures to the quasi biennial 
oscillation in equatorial winds in the lower stratosphere, a strong influence from equatorial winds in 
the upper stratosphere is found. We suggest that this may be associated with the rate of onset of the 
westerly phase of the semi annual oscillation in the upper stratosphere, although the observations alone 
cannot provide a conclusive, causal relationship. The main diagnostic tools employed are correlation 
studies and composite analysis. The need for continued high quality, equatorial wind measurements at 
all stratospheric levels is also stressed. 

Keywords: Stratosphere Quasi biennial oscillation Semi annual oscillation Northern hemi- 

sphere winter Stratospheric warming Interannual variability 


1. Introduction 

The Northern Hemisphere (NH) stratospheric winter circulation displays substantial 
inter-annual variability (Labitzke 1982). Some winters are extremely disturbed and are 
accompanied by stratospheric sudden warmings, in which the polar temperature increases 
by 20° C or more in just a few days (Andrews et al. 1987). In contrast, other winters have 
fewer warming events and the polar vortex remains cold and undisturbed. Variability 
associated with stratospheric warmings is the primary source of variability in the NH 
winter lower stratosphere. The variability of NH winter stratospheric temperatures is 
also closely coupled to the Arctic Oscillation, the leading mode of variability in the 
troposphere (Thompson and Wallace 1998, Baldwin and Dunkerton 1999, Ramaswamy 
et al. 2000). 

Global, annual-mean temperatures in the lower stratosphere have fallen by ~0.6°C 
per decade over recent decades (e.g. Pyle et al. 1999). It is important to characterise nat- 
ural variability if we are to assess how much of the observed temperature trend is directly 
related to human activity. It is also important because temperature directly influences 
the processes that lead to ozone trends in this region. Stratospheric warmings result in 
a weak, disturbed vortex with warm temperatures and substantial downward transport 
of ozone-rich air. Conversely, a strong, stable vortex with cold temperatures will result 
in large chemical ozone loss (due to the temperature-dependence of the destruction reac- 
tions) and weak transport of ozone (Chipperfield and Jones 1999). A number of studies 
have found evidence for long-term changes in the stratospheric circulation, including a 

* Corresponding author; Rutherford Appleton Laboratory, Chilton, Didcot, Oxon., OX11 0QX, UK. 
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strengthening of the Arctic vortex (Zurek et al. 1996) and changes in the NH geopoten- 
tial heights (Labitzke and van Loon 1995). Any such changes would also influence the 
long-term trends in ozone and temperature (Hood et al. 1997). 

Although stratospheric sudden warmings have been observed and documented for 
many years, the factors that control their variability are not well understood. An essential 
requirement for their development is the presence of quasi-stationary planetary wave 
disturbances. These waves are generated in the troposphere and propagate vertically and 
horizontally into the stratosphere, where they can be of sufficient amplitude to displace 
the vortex away from the pole. These waves are the main mechanism for heat transfer to 
the polar region and associated with each vortex displacement there is a rapid increase 
in temperature at the pole. 

The dependence of the winter-time average polar temperature and vortex strength 
on the phase of the quasi biennial oscillation (QBO) in the lower stratosphere was noted 
nearly twenty years ago (Holton and Tan 1980, 1982). The QBO is an oscillation of the 
equatorial zoned wind in the lower stratosphere. They oscillate between easterlies and 
westerlies with an average period of 28 months. The maximum amplitude of the oscillation 
is around 25-30km. When the equatorial winds at 20-30 km are in an easterly phase, the 
northern polar vortex is generally warmer, more disturbed by waves and disruption of 
the vortex by major mid- winter warmings is more likely (see also Dunkerton and Baldwin 
1991, Baldwin and Dunkerton 1998). The conventional explanation for this equatorial 
influence on polar temperatures is that the QBO winds in the lower stratosphere influence 
the background mean flow, which then affects the propagation of planetary scale waves 
(e.g. O’Sullivan and Young 1992; O’Sullivan and Dunkerton 1994; Hamilton 1998; Niwano 
and Takahashi 1998; Baldwin et al. 2000). More specifically, the QBO winds in the lower 
stratosphere determine the position of the zero wind line near the equator, which acts 
as a wave-guide for the planetary wave propagation. This results in enhanced poleward 
heat transfer during an easterly QBO phase and weaker transfer during a westerly phase. 

While this mechanism is generally accepted, it has also been recognised that this 
pattern of ‘warm disturbed easterly phase’ / ‘cold undisturbed westerly phase’ does 
not always hold up, particularly in late winter (Hamilton 1998; Baldwin et al. 2000). 
Labitzke and van Loon (1988) and Naito and Hirota (1997) have noted that the periods 
when the Holton-Tan relationship holds up well coincide with periods when the 11-year 
solar sunspot cycle is at its minimum phase (defined, for example, by the 10.7cm solar 
flux). Conversely, it appears to substantially weaken or reverse during periods of solar 
maximum. While the solar cycle / QBO link and the solar cycle / NH polar temperature 
link has been studied extensively (e.g. Labitzke and van Loon 1988; Gray and Dunkerton 
1990; Teitelbaum and Bauer 1990; Salby and Shea 1991; Dunkerton and Baldwin 1992; 
Kodera 1991, 1993, 1995; Balachandran and Rind 1995; Rind and Balachandran 1995; 
Labitzke and van Loon 1996; Naito and Hirota 1997; Haigh 1999; Shindell et al. 1999), 
the precise mechanism that relates (1) the solar cycle, whose primary influence is in the 
region of the equatorial stratopause (^50 km), (2) the QBO, which has its maximum 
amplitude in the lower equatorial stratosphere (^20-30 km) and (3) the winter-time lower 
stratospheric polar temperatures, has yet to be clearly elucidated. 

In this paper, we aim to explore the possible links between these three phenomena 
by examining the relationship between polar temperatures and equatorial winds using 
an equatorial wind data set that measures much higher than the radiosonde dataset 
more usually employed. While the radiosonde observations extend to ~20-30km, the 
rocketsonde dataset employed in this study extends to 58km (Dunkerton and Delisi 1997; 
Dunkerton et al. 1998). The rocketsonde dataset thereby includes not only the lower 
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Figure 1. Correlation for the periods (a) 1964/65 - 89/90 and (b) 1955/56 - 1998/99 between JF North 
Polar (NP) temperature at 24km and monthly-averaged equatorial wind in each month from January of 
the previous year to February of the same year. Contour interval is 0.1. Dotted contours denote negative 
values. Shading indicates values are significant at the 95% confidence level (no values were significant at 
the 99% confidence level). For details of the significance test applied, see appendix A. NP temperatures 
are from the Berlin stratospheric analyses. Equatorial winds are from radiosonde ascents. 


stratospheric winds (20~30km) usually associated with the observed QBO modulation 
of the polar temperatures but also the higher altitudes that encompass the stratopause 
region, where the direct solar influence is greatest. In section 2 we briefly review the 
results of correlations between NP temperatures and the lower stratospheric winds from 
the radiosonde dataset. In section 3 the rocketsonde dataset is described and a detailed 
analysis is presented of the relationships between NP temperatures and equatorial winds 
at ail heights to 58km using the rocketsonde dataset. This analysis includes correlation 
studies, regression analyses, and composite studies. The main conclusions of the study 
are presented in section 4. 


2. Radiosonde Analysis 

Figure la shows the correlation over 24 winters (1964-90) between the January- 
February average (JF) North Polar (NP) temperature and the equatorial zonal winds 
(u eff ) from radiosonde observations in each month from January of the previous year to 
February of the same year. We choose to use JF temperature since most major midwinter 
warmings occur in these months. The polar temperatures are anti-correlated with the 
lower stratospheric equatorial winds in the previous months with the region of maximum 
correlation gradually descending through the atmosphere with time, corresponding to 
the descent of the QBO winds. This pattern is consistent with the mechanism described 
above (hereafter referred to as the Holton-Tan mechanism), with easterly (negative) wind 
anomalies associated with warm (positive) temperature anomalies and vice versa. For the 
26-year period 1964-90, the maximum correlation coefficients are of order -0.4. However, 
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when the correlation is repeated using data from 43 winters (1955-99, see figure lb), 
which is the full extent of the available radiosonde data, the maximum correlation drops 
to only 0.25. 

Figure 2a shows the actual JF north polar temperatures and corresponding lower 
stratospheric equatorial wind for the period 1964-90 that were used to determine the 
correlations of figure la. It shows that the Holton-Tan relationship of negative correlation 
between the two signals is present for certain periods but not others. For example, the 
two signals are anti-correlated during the periods 1965-68, 1971-78 and 1986-88, but in 
the intervening years either the opposite is true (e.g. 1969, 1980 and 1984-85) or there 
is no apparent relationship (e.g. 1982-83). It is these latter periods that degrade the 
correlation, resulting in the relatively low correlations seen in figure la. A scatter plot of 
the data (figure 2b) further demonstrates the weakness of the resulting correlation. 

It is possible that the difference in maximum correlation displayed by the analysis 
of the two different time periods in figure 1 can be understood as the difference in the 
relative number of solar minima and maxima in the two selected time periods (Labitzke 
and van Loon 1988; Naito and Hirota 1997). While the 1955/90 period encompasses an 
equal number of solar minimum and solar maximum periods, the 1964/90 period has a 
slightly larger influence from solar minimum than solar maximum periods. 


3. Rocketsonde Analysis 
(a) The data 

Regular rocketsonde ascents were made from Kwajalein (8°N) during the period 
1969-90 and from Ascension Island (8°S) during the periods 1962-77 and 1979-89. Occa- 
sional gaps of up to a few months duration are present in both records. In order to obtain 
a single, unbroken time series of equatorial winds, the monthly- averaged data from the 
two stations were first de-seasonalised. This was done by calculating and subtracting the 
climatological monthly mean at each height and at each station. Gaps of up to three 
months duration in the records for each station were then filled by linear interpolation. 
For those heights and months when data were available for both stations, the two data 
sets were then combined by averaging them together. When data were only available from 
one station, that data was used. In this way, an unbroken time series of de-seasonalised 
‘equatorial’ (8°S-8°N averaged) winds was achieved from January 1964 to February 1990 
for the height region 20-58 km. 

Figure 3a shows the time-height cross-section of the raw monthly-averaged data from 
Kwajalein. The QBO is the dominant signal in the lower stratosphere (20-35 km). The 
maximum amplitude of the oscillation is around 25 km, with the wind varying between 
-30 ms” 1 and -4-20 ms” 1 . Above 35 km the semi-annual oscillation (SAO), with a period 
of six months, dominates. Figure 3b shows the combined Kwajalein and Ascension Island 
data set, which is derived from the de-seasonalised data as described above. (Note that 
the SAO has been removed by the deseasonalising process). The QBO period has been 
further highlighted in figure 3c using a 9-60 month bandpass filter. 

(6) Correlation studies 

Using the combined de-seasonalised rocketsonde wind data (figure 3b), we have 
carried out a correlation study to examine the link between the equatorial winds and the 
Northern Hemisphere polar temperature in winter. Figure 4 shows the strength of the 
correlation between the observed JF North Polar temperature at 24km from the Berlin 
stratospheric analyses (courtesy K. Labitzke) and the rocket-sonde equatorial winds at 
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Zonal — wind anomaly (ms' 1 ) 


Figure 2 . (a) Time-series of JF North Polar temperature anomaly (K) at 24km 30hPa) (dotted line) 

and equatorial wind anomaly (ms' 1 ) at 24km in the preceding December (solid line), (b) Scatter plot 

of the data. 


each height and month leading up to January/February (i.e. identical to figure la but 
using the rocketsonde data). Note first the negative correlation of ~0.4 at 20-35 km which 
slowly propagates downwards with time. This was the correlation noted by Holton and 
Tan nearly 20 years ago and is virtually identical to the results of the correlations using 
radiosonde wind data (figure la). 

However, there is also significant correlation at higher levels. In the 30-50 km region 
there is a positive correlation of ~0.3-0.4, also descending slowly with time. This change 
in sign of the correlation with height reflects the next incoming phase of the equatorial 
wind QBO (see figure 3). Above 50 km, there is a region of relatively high correlation in 





Figure 3. (a) Monthly-averaged zonal wind (ms" 1 ) derived from rocketsonde ascents at Kwajalein 

(8°N). (b) De-seasonaiised monthly-averaged zonal wind dataset derived by combining data from both 
Kwajalein and Ascension Island as described in the text, (c) As in pane! (b) after passing through a 9-60 
month bandpass filter. Zero value is denoted by contour line. 


September/October. At 52 km in September, the correlation exceeds -0.6 and is signifi- 
cant at the 99% confidence level. We note here that there can be no suggestion that this 
high correlation is simply a response of the equatorial winds to the winter polar warming, 
since the strongest correlation is found with the previous September, 4-5 months before 
the major warming takes place. 

(c) Regression analyses 

Linear regression analyses were carried out in order to assess the effectiveness of 
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Figure 4. As figure 1 but using rocketsonde data that extends to 58km instead of radiosonde data. 
Light and dark shading denote regions where values are significant at the 95% and 99% confidence levels 

respectively. 


the rocketsonde equatorial winds (n eq ) at different heights and in different months as 
a predictor of the winter-time polar temperature. Initially, simple linear regression was 
used to fit the JF NP temperature to n eq at (a) 24km in December (as an indicator of 
the phase of the QBO in the lower stratosphere), and (b) 52km in September. The best 
fit to w eq at 52 km was found to describe 37% of the variance in the polar temperature. 
By contrast, the best fit to n eq at 24km was found to describe only 16% of the variance. 

A series of multiple linear regressions was then carried out in order to find the two 
parameters that best predict the polar temperature. In each analysis, u e)? for a given 
height and month was employed as a fixed parameter, while u eq for every other available 
height and month was used in turn as the second parameter. In this way, we sought to 
maximise the amount of the variance in the polar temperature that could be described. 
In the first such analysis, n eq at 24km in December was employed as the fixed parameter. 
It was found that u eq at 52 km in September was most effective as the second parameter, 
with a total of 41% of the variance in the polar temperature being described. 

Repeating the analysis with u eq at 52 km in September as the fixed parameter, 
it was found that the most effective second parameter was \i eq at 48 km in the same 
month. A total of 51% of the variance in the polar temperature was described. This is 
an interesting result on two accounts. Firstly, the greatest additional contribution did 
not come from the lower stratosphere - the region in which the Holton-Tan mechanism is 
at work. Secondly, the two parameters consist of the equatorial wind immediately above 
and immediately below 50 km in September. In figure 4 it can be seen that there is a 
sharp change in the sign of the correlation at this point. This suggests that the vertical 
gradient of the wind at — 50 km in September may be in some way linked to the North 
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Polar temperature in the following winter. 

In order to investigate this link further, the vertical gradient of zonally-averaged 
zonal wind dn eq /dz was derived from the raw rocketsonde data using centred differ- 
ences. These values were then de-seasonalised and combined in the same manner as that 
described above for the equatorial winds. The correlation between the JF North Polar 
temperature and du eq f dz was found to peak at 50km in September, with a correlation 
of -0.68. Using simple linear regression it was found that du eq /dz at 50 km in September 
could describe 46% of the variance in the JF North Polar temperature. The above mul- 
tiple linear regression was repeated, using du cg /dz at 50 km in September as the fixed 
parameter and all possible values of both n €q and du eq /dz as the second parameter. It 
was found that n eq at 44km in October was most effective as the second parameter, with 
a total of 62% of the variance in the polar temperature being described. 

These results suggest a link between the JF polar temperature and the equato- 
rial wind and its vertical gradient in the upper stratosphere in the preceeding Septem- 
ber/October. This strong link is evident in figure 5a, which shows the JF polar temper- 
atures for 1965-90, overlaid by d\i eq /dz at 50 km in the preceeding September. The two 
signals are reasonably well anti-correlated over most of the period. A comparison of the 
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scatter-plot of the data (figure 5b) with the corresponding plot in figure 2b for the lower 
stratospheric winds (u e9 at 24km) clearly demonstrates that the relationship between the 
JF polar temperature and du e(? /dz at 50km is stronger than with \i €q at 24km. 

(d) Composites 

In figure 6 we show a series of latitude - height 'composite differences’ of JF winds 
using daily wind data derived from the National Centers for Environmental Prediction 
/ National Center for Atmospheric Research (NCEP/NCAR) 1200 UTC heights and 
temperatures (Kalnay et al. 1996). Firstly, composites were derived by averaging together 
the JF wind distributions from the five warmest and five coldest Northern Hemisphere 
winters between 1964/65 and 1989/90. The appropriate years were selected on the basis 
of the North Polar temperatures at 24km. We do this to highlight the differences in the 
JF wind distribution between years in which major warmings have occurred and those in 
which they have not. The ‘cold’ composite minus the ‘warm 1 composite is shown in figure 
6a. During cold winters the westerly polar vortex is less disturbed by planetary waves and 
their accompanying warming events and hence there is a strong, high-latitude westerly 
anomaly, with maximum composite differences of 35 ms” 1 at 32 km (~ 10 hPa). There 
is a small westerly anomaly at the equator, consistent with the Holton-Tan mechanism, 
peaking at around 11 ms" 1 . 

Secondly, the wind distributions were composited on the basis of the direction of the 
equatorial wind at 24 km in January. These composites test the Holton-Tan mechanism 
whereby the polar temperatures and vortex are influenced by the equatorial wind in the 
lower stratosphere. The ‘west QBO’ minus ‘east QBO 1 composite differences in figure 6b 
show a high latitude westerly wind anomaly peaking at ~ 15 ms" 1 at 32km, consistent 
with a stronger westerly polar jet during the westerly phase of the QBO at 24km. 

Finally, composites were constructed on the basis of du e<7 /dz at 50km in September. 
The difference between the ‘large shear’ and ‘small shear 1 years has a maximum amplitude 
of -25 ms -1 at high latitudes (figure 6c). Hence, years with stronger than average vertical 
wind shear near the equatorial stratopause tend to be characterised by a cold, more stable 
westerly polar vortex, and vice versa. These composites confirm that, while by no means 
accounting for the full amplitude of the differences in figure 6a, the composite differences 
based on du eq /dz at 50km in September are substantially larger than the composite 
differences constructed using the 24km equatorial winds. 

In figure 7 we show the height-time evolution of the equatorial wind from the NCEP 
analyses for the ‘warm’ and ‘cold 1 composites described above, and the difference between 
them. The major feature of the equatorial wind series between 40-60 km is the semi- 
annual oscillation (SAO), in which the equatorial wind direction reverses from easterly 
to westerly and back again on a time-scale of approximately 6 months. September is the 
onset of a westerly phase at around 50km. There is a notable difference in the nature of 
this westerly phase between the ‘warm 1 and ‘cold 1 composites. In the ‘warm 1 composite, 
the westerly phase extends deeper in the atmosphere, with the zero wind line at its lowest 
extent in November reaching —30km compared with only ~ 40km in the ‘cold 1 composite. 
This is a result of a more rapid descent of the westerly phase during September and 
October. The vertical wind shear at — 50km in September reflects this difference, being 
weaker in the ‘warm 1 composite than in the ‘cold 1 composite (as we have already noted 
in the discussion of figure 6c above). 

Figure 7 suggests that the difference between the warm and cold composites may 
be associated with the relative strength and penetration of the SAO westerly phase 
around September/October. The westerly phase of the SAO arises as a result of energy 
deposition by westerly-phase vertically propagating equatorial waves e.g. gravity and 




Figure 6. Latitude- height distributions of JF zonally-averaged zonal wind (ms' 1 ) from the NCEP 
analyses for the period 1964/65 - 1989/90. (a) composite difference of 5 coldest minus 5 warmest NH 
winters as determined by the 30 hPa JF NP temperatures, (b) composite difference of 5 most westerly 
minus 5 most easterly QBO phase based on equatorial winds at 24km in the previous December, (c) 
composite difference of 5 strongest minus 5 weakest du C(J /dz at 50km in the previous September. Contour 
interval is 5 ms -1 . Negative (easterly) values are denoted by dotted contours. Light and dark shading 
denote regions where values are significant at the 95% and 99% confidence levels respectively. The 
significances of the values were tested by performing a Student’s t-test on the difference between the two 

mean wind distributions. 
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Figure 7. Time-series of monthly- averaged, zonally- averaged zonal wind at the equator (ms -1 ) derived 
from the NCEP analyses, (a) Composite of 5 warmest NH winters, as determined using 30 hPa NP 
temperatures, (b) As for panel (a) but for the 5 coldest winters, (c) The difference between the two 
composites i.e. panel (a) minus panel (b). Contour interval is 5 ms -1 . Negative (easterly) values are 
denoted by dotted contours. Light and dark shading denote regions where values are significant at the 
95% and 99% confidence levels respectively. The significances of the values were tested by performing a 
Student’s t-test on the difference between the two mean wind distributions. 
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Kelvin waves (Andrews et al 1997). These waves are able to propagate relatively easily 
through lower stratospheric easterlies. However, when the waves reach heights (~ 50km) 
where the background zonal winds are of similar magnitude and direction to their own 
phase speeds, the waves are unable to propagate further. They are dissipated and their 
energy is transferred to the background flow, thus forcing the local background flow 
even more westerly. In this way, the westerly SAO phase gradually descends through 
the atmosphere with time. The rate of descent of the westerly phase, and hence the 
vertical wind shear at 50km, is therefore dependent on the QBO in the lower stratosphere 
(Lindzen and Holton, 1968; Gray and Pyle 1989; Hirota et al 1991; Kennaugh et al 
1997), since the relative ease of vertical propagation of the westerly- phase waves depends 
on the wind speed and direction in the lower stratosphere. Strong easterlies at 20-30 km 
allow enhanced penetration and hence more effective westerly forcing higher up (and vice 
versa). The resulting QBO modulation has been observed at mesospheric levels in the 
HRDI wind observations (Burrage et al 1996) and is evident in figure 3c which shows a 
weak QBO signal extending into the upper stratosphere. It is for this reason that we see 
evidence of a QBO-like variation in dn eq /dz at 50 km in figure 5 and a fairly substantial 
westerly anomaly at the equator in the composite difference of figure 6c, similar to that 
in figure 6b. 


4. Conclusions 

We have presented evidence for a dependence of the North Polar temperature on 
the equatorial wind and its vertical shear in the upper stratosphere in the previous 
September/October. This dependence appears to be related to the timing of onset of the 
westerly phase of the SAO. The evidence provided here strongly suggests an influence 
on polar winter temperatures from the upper stratosphere, in addition to that from the 
lower stratospheric QBO. This upper height region also coincides with the region of the 
atmosphere where the solar cycle is believed to exert its strongest influence. While a 
data correlation study such as this cannot, on its own, establish a causal relationship 
nor provide the mechanism for this relationship (if, indeed, one exists) it nevertheless 
suggests that a better understanding of the equatorial upper stratosphere, its relationship 
to the QBO lower down and to the polar lower stratosphere will be essential if we are 
to fully understand the nature of the equatorial influence on polar temperatures and 
hence on the Earth’s climate. The study also highlights the need for continued high 
quality wind measurements of the equatorial upper atmosphere. The 1990s have been 
an unusual decade, with very few mid-winter warmings, which may be an indicator of 
climate change. Unfortunately, our analysis could not be extended into the 1990s because 
of the lack of adequate equatorial wind data in the region 40-60 km. 


Appendix 

In figures 1 and 4 the significance of eac h poi nt on the correlation plot was tested by 
calculating the value of the statistic: t = r^J where r is the correlation between the 
two time-series and N is the number of data points in the time-series. In the null case of 
no correlation, this statistic is distributed like Student’s t-distribution with N-2 degrees of 
freedom (Press et al. 1992). In order to compensate for the autocorrelation of the wind 
data, the effective sample size N e was substituted for N, where N e = N unless pi > 0 

when N e = N and p\ is the lag-one autocorrelation coefficient (Wilks 1995). 
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‘Annular modes’ are deep, longitudinally symmetric patterns 
which dominate the circulation variability of both hemi- 
spheres from the Earth’s surface to the mesosphere' \ These 
ring-like patterns are characterised by fluctuations in pressure 
of one sign over the polar caps and the opposite sign at low 
latitudes. Their simple structure appears to be a fundamental 
fluid dynamical consequence of the Earth’s size and rotation 
rate 3 . More than two decades of global satellite observations 
allow us to examine the time-height development of the an- 
nular modes. Here we present the first observations that show 
that fluctuations in the annular modes often begin in the 
lower mesosphere and propagate downward, over a period of 
up to three months, to affect surface weather and climate pat- 
terns. These results show that variations in the stratospheric 
annular modes are coupled to the surface, and that in the up- 
per stratosphere, the two hemispheres are coupled. Perturba- 
tions and trends in the circulation of the stratosphere and 
mesosphere, which may be caused by natural variability, in- 
creasing greenhouse gases, ozone depletion, volcanic aerosols, 
or 11-year solar variability, may have a larger influence on 
surface weather and climate trends than has been recognised. 

Variability in the circulation of the stratosphere (-11-48 
km) and lower mesosphere is largest during the winter, when the 
cold, cyclonic polar vortex varies in strength and is disturbed by 
planetary-scale waves. These waves originate in the troposphere, 
transport momentum upward, and interact with the strato- 
spheric flow. In the Northern Hemisphere, where planetary wave 
amplitudes are larger, the process of wave-mean interaction can 
lead to the sudden and dramatic slowing and reversal of the 
westerly polar vortex, even in midwinter (called a major strato- 
spheric warming). The southern vortex is colder, stronger, and is 
less disturbed by waves. 

Theory 4 ’ 5 , observations 6 , and modelling experiments 7 show 
not only that the circulation of the troposphere and stratosphere 
are coupled, but that the influence on the mean flow tends to be 
downward, even while circulation changes are driven by tro- 
pospheric waves that propagate upward.' Stratospheric circula- 
tion anomalies alter the propagation of wave energy from the 
troposphere, and the resulting changes to the flow tend to bur- 
row downward, often into the troposphere. This ‘downward 
control’ 4 has important implications for long-term climate vari- 
ability. Surface weather and climate appear to be sensitive to the 
state of the stratosphere, and the stratosphere is sensitive to 
changes in the concentrations of greenhouse gases, ozone, and 
volcanic aerosols. Annular modes provide a conceptual model 
with which to examine and quantify the observed circulation 
changes which span the troposphere and stratosphere. 

Upward propagating planetary-scale waves are ducted by a 
waveguide which depends on the longitudinally averaged west- 
erly wind, u. Anomalies (defined as deviations from the average 


climatology) in the stratospheric u field act to refract upward 
propagating planetary waves, alter the locations where wave en- 
ergy is deposited, and initiate a positive feedback in which the 
waves are refracted at successively lower altitudes 8 . Variations in 
u tend to occur over deep layers in the vertical, are of opposite 
sign at high and low latitudes 9 ’ 10 , and correspond to changes in 
the annular modes because u varies as the north-south gradient 
of geopotential height surfaces. 

Annular modes provide a better measure of vertical coupling 
to the troposphere than longitudinally averaged fields , such as 17. 
The use of geopotential height surfaces to define the annular 
modes has the advantage of providing a three-dimensional pic- 
ture of vertical coupling. The annular modes are typically de- 
fined within the troposphere as the leading empirical orthogonal 
function (EOF) of slowly varying (e.g. monthly mean) cold- 
season geopotential anomalies. The surface, or 1000-mbar, 
northern annular mode (NAM, Fig. li) is known as the ‘Arctic 
Oscillation’ 1 1 and its southern counterpart is the southern an- 
nular mode (SAM). The NAM pattern is robust, in that it is in- 
sensitive to the details of the calculation (months or data levels 
used, temporal averaging, data set, etc.) 6 . Over the Atlantic sec- 
tor, it is recognised as the North Atlantic Oscillation * (NAO); it 
has been argued that NAO and annular mode are two paradigms 
describing a single phenomenon 13 . Variations in the NAM index 
are related to hemispheric-scale effects not just on weather pat- 
terns, but on extreme weather events such as heavy rainfall, 
snow, wind, and droughts . 

We extend the concept of the annular modes into the lower 
mesosphere (Fig. 1) by defining the annular modes as the lead- 
ing EOF for each of 26 pressure levels from 1000 to 0.316 mbar 
(~57 km). We use November-April data to define the cold sea- 
son in the Northern Hemisphere, but September-December in 
the Southern Hemisphere (the time during which the southern 
vortex is decreasing in size and is susceptible to influence from 
upward propagating planetary-scale waves). The annular modes 
(Fig. 1) account for 20-30% of the (90-day low-pass filtered) 
temporal variance of geopotential height in the troposphere, but 
more than half the variance above the tropopause. The smaller 
variability of the Southern Hemisphere stratospheric circulation 
results in a smaller amplitude of the SAM, approximately 2 A that 
of the NAM. 

We define a daily index of the annular modes for each level 
by regressing daily geopotential anomalies onto the annular 
mode pattern for each level. The resulting indices are similar to 
the principal component time series from the EOF analysis, but 
have daily resolution spanning the entire data record (Fig. 2). 
The troposphere in both hemispheres experiences many fluc- 
tuations in the annular modes that are evanescent above the 
tropopause. However, events during the cold seasons often show 
deep, strong coupling from the lower mesosphere to the surface. 
The tendency is for such events to be initiated at the highest 
levels. Such downward propagation is seen here for the first time 
above 10 mbar and in the Southern Hemisphere. The Northern 
Hemisphere shows more consistent vertical coupling and 
downward propagation of the annular mode, which is to be ex- 
pected because the northern vortex can break down at any time 
during the winter. Stratospheric warmings correspond to a weak 
polar vortex, shown in red 6 . Typically, these events begin in the 
mesosphere and propagate downward over a period of a few 
weeks, often (but not always) reaching the Earth’s surface. The 
negative anomalies (red) are more likely to penetrate into the 
troposphere. 

The time scale for events in the Southern Hemisphere is 
somewhat longer and corresponds to the breakdown of the po- 
lar vortex during spring. As the vortex is reduced in size, begin- 
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ning in August-September, anomalies in the SAM tend to move 
slowly downward over a period of 2-3 months. Vacillations in 
the SAM are not as consistent as the NAM, nor are they confined 
entirely to the cold season. Although the NAM has a very small 
amplitude during the warm season, the SAM in both the tro- 
posphere and stratosphere, has substantial fluctuations through- 
out the year. 

The northern winter of 1998-1999 was the only winter ever 
to experience two major midwinter stratospheric warming 
events (Fig. 3). The alternating, descending positive and negative 
anomalies (shown with no time filtering) provide a striking il- 
lustration of the balance between wave-mean interaction as the 
vortex is diminished during two warming events, and radiative 
cooling (aided by nonlinear wave-wave interactions' 6 ) as the 
vortex rebuilds during the following weeks. The wave-mean 
interaction process appears to be similar to that which drives the 
descending easterly and westerly wind regimes of the quasi- 
biennial oscillation (QBO) in the equatorial stratosphere. Both 
phases of the QBO are forced by a broad spectrum of equatorial 
waves. A similar wave-driven oscillation has been observed in 
the stratosphere of Jupiter' 7 and has been hypothesised to occur 
in the interior of the sun' 8 . In each of these phenomena vertically 
propagating waves produce an orderly, downward propagating 
oscillation from chaotic forcing. The annual cycle allows only 
brief windows during which the modes descend, but ‘perpetual 
January' GCM simulations 7 show irregular, alternating, de- 
scending wind regimes (vacillation cycles). 

Anomalies in 17 are observed' * and simulated to move 
northward as they descend over a period of ~2 months, corre- 
sponding particularly well to stratospheric warmings. Since the 
NAM pattern is much tighter at 10 mbar than at 1 mbar, corre- 
sponding anomalies in 10-mbar 77 are found at higher latitudes. 
For each level we calculated the geostrophic wind anomaly of the 
longitudinally averaged SAM and NAM. As the NAM (Fig. 4b) 
descends, maximum wind anomalies are found at lower levels 
and higher latitudes, beginning at ~50°N and reaching ~70°N in 
the lower stratosphere. Downward propagation of the NAM 
corresponds to northward and downward propagation of zonal 
wind anomalies. Even though the latitudinal structure of the 
NAM is fixed and cannot describe northward propagation at a 
fixed level, the vertical structure of the NAM provides a partial 
explanation for this propagation. In the southern hemisphere 
the scale of the SAM pattern also changes in the vertical, so a 
similar, but smaller, effect is seen (Fig. 4a). 

Ozone depletion and increasing greenhouse gases are ex- 
pected to cool the stratosphere, increase the strength of the 
stratospheric polar vortex, and influence the phase of the Arctic 
Oscillation 22, . We examine the relationship between annular 
mode indices at 10 and 1000 mbar, for each hemisphere (the 
observational record in the Northern Hemisphere (Fig. 5) per- 
mits analysis to 10 mbar since 1958 b ). The correlations (0.54 in 
the Northern Hemisphere, 0.66 in the Southern Hemisphere) do 
not prove cause and effect, but do suggest that a trend in the 
strength of the stratospheric polar vortex would be mimicked at 
the surface, as has been simulated in GCMs 24,i . Previous model- 
ling studies have demonstrated that a trend in the Arctic Oscilla- 
tion may also be caused by land use changes , raising the ques- 
tion of the relative influence on climate of tropospheric and 
stratospheric forcing. The NAM and SAM indices tend to vary 
together in the upper stratosphere above 5 mbar (Fig. 5c) This 
suggests that either the same forcing is acting on both hemi- 
spheres, the hemispheres are linked dynamically, or both. 

Our results indicate that perturbations to the temperature 
and wind fields of the stratosphere and mesosphere can affect 
surface weather and climate. Such perturbations occur on a vari- 


ety of time scales, due to natural variability, volcanic eruptions, 
the 11 -year solar cycle, ozone depletion, and increasing concen- 
trations of greenhouse gases. Wave-mean interaction provides a 
mechanism whereby changes to the flow at high levels and low 
atmospheric densities (the mesosphere contains less than 0.1% 
of the mass of the atmosphere) can be amplified and communi- 
cated downward. We also find evidence of inter-hemispheric 
coupling of the annular modes in the upper stratosphere. 

Methods 

Data 

At 10 mbar and below we use National Centers for Environmental 
Prediction (NCEP) ‘reanalysis’ data. Prior to 1 1 July, 1993 from 6.81 to 1 
mbar we use TOVS geopotential thickness fields added to the NCEP 10- 
mbar field. After 1 1 July, 1993, up to 0.316 mbar we use UKMO thickness 
fields above 10 mbar. 
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Figure 1 Northern (left) and Southern (right) annular mode patterns in 
geopotential height (m). The patterns represent the first EOF of 90-day 
low pass filtered anomalies for November-April in the Northern 
Hemisphere and September- December in the Southern Hemisphere. 
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Figure 2 Time-height development of the NAM (top 3 panels) and SAM (bottom 3 panels) during 1979-1999. The ‘"dices have been 45-day low pass filtered The 
indices are nondimensional, with a value of 1 .0 corresponding to the patterns in Figure 1 The contour interval is 0.5. wrlh values between -0.5 and 0 5 unshaded 





Figure 4 Geostrophic wind anomalies, y (ms 1 ) corresponding to the longitudinally 
averaged annular modes in Figure 1 . a, SAM. b, NAM 


Figure 5 Comparison of NAM and SAM indices. All time series have been 
smoothed with a triangular filter (250-day half width) a, Northern Hemisphere 
1 000 mbar (black) and 1 0 mbar(red) ; b, as In a, but Southern Hemisphere, c , 
Comparison of 1 -mbar NAM (black) and SAM (red) 
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Abstract 

The leading mode of variability of the lower atmosphere circulation in the Northern Hemisphere is a 
largely zonally symmetric mode known as the Arctic Oscillation. We calculate Arctic Oscillation (AO) 
indices on a range of levels from 1000 to lOhPa by means of a principal component analysis of NCEP 
daily geopotential height anomalies. We find the apparent downward propagation of anomalies noted by 
other authors to be statistically significant compared to a red noise model. By examining histograms of 
these indices for each month, we note that the distribution of the index is generally close to Gaussian in 
the troposphere. In the stratosphere, however, the index is negatively skewed in the winter, and 
positively skewed in the spring. We conclude that the positive skewness in April results from the 
coexistence of distinct summer and winter circulation states, and by examining polar stratospheric 
temperatures we conclude that the negative skewness in January is due to the radiatively-determined 
limit on the vortex strength. This radiative limit responds only slowly to anthropogenic forcing, whereas 
changes in planetary wave forcing could have a large impact on the number of warm events. This leads 
to the hypothesis that the vortex strength may respond nonlinearly to anthropogenic forcing, which is 
supported by an observed change in the shape of the histograms of 20-200hPa AO indices in January 
over the past 40 years. In view of the downward propagation discussed, we speculate that this 
nonlinearity extends to the tropospheric response. 
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1. Introduction 

The leading mode of variability in the Northern Hemi- 
sphere troposphere and stratosphere is an approximately 
zonally symmetric pattern known as the Arctic Oscilla- 
tion [Thompson and Wallace , 1998]. This pattern, con- 
ventionally defined as the first empirical orthogonal func- 
tion (EOF) of cold season sea level pressure, has an index 
which is closely associated with the more regional North 
Atlantic Oscillation [ Wallace , 2000]. Thompson and Wal- 
lace [1998] describe how the AO index is closely corre- 
lated with the strength of the stratospheric polar vortex, 
and this relationship is explored more fully by Baldwin 
and Dunkerton [1999], who perform a principal compo- 
nent analysis on the geopotential height field on five levels 
from 1000 to lOhPa. It has long been known that warming 
events propagate downwards in the stratosphere [e.g. Mur- 
ray y I960], but Baldwin and Dunkerton [1999] were the 
first to note such downward propagation in a zonal in- 
dex, with anomalies of both signs traveling from the lOhPa 
level to the surface in around three weeks. However, they 
do not assess the statistical significance of this result. Such 
an analysis is called for since they examine a signal with a 
20-day lag in 90-day low pass filtered data. Since the pub- 
lication of Baldwin and Dunkerton [1999] such downward 
propagation has also been noted by Yamazaki and Shinya 
[ 1999] and P. McCloghrie (pers. comm.) in GCMs. Shin- 
dell et ai [1998] suggest that Arctic Oscillation anomalies 
may propagate downwards by means of an interaction be- 
tween upward propagating planetary waves and the mean 
flow. This suggestion is supported by the work of Kuroda 
and Kodera [1999] who note that zonal wind anomalies 
appear first in the mid-latitude stratosphere, and propagate 
polewards and downwards. They show the Eliasscn-Palm 
flux, and note that it is consistent with the downward prop- 
agation resulting from a wave/mean-flow interaction. 

Many authors have analysed geopotential height datasets 
by computing probability density functions (PDFs) of the 
system in the phase-space of one or more principal com- 
ponents. The field most often analysed is 500hPa geopo- 
tential height in the Northern Hemisphere winter. Both 
Hansen and Sutera [ 1 992] and Coni et ai [ 1 999] note ap- 
parent multimodality in the first principal component of 
this field in a model and observations respectively. By 
contrast Kimoto and Ghil [1993] and Cheng and Wallace 
[ 1993] find the PDF to be unimodal, but significantly dif- 
ferent from a Gaussian, by comparison with a red noise 
model. The skewness of Northern Hemisphere winter tro- 
pospheric geopotential height has been analysed on a grid- 
point basis by Nakamura and Wallace [1991]. They find 
positive skewness to the north of the stormtrack, and nega- 
tive skewness to its south. Baldwin and Dunkerton [1999] 


mention that the skewness of the Arctic Oscillation index 
at lOhPa is slightly negative overall, but they do not con- 
sider how this skewness changes through the year, or in- 
deed whether it is significantly different from zero. We 
estimate PDFs of the daily Arctic Oscillation indices for 
each month on a range of levels between 1 000 and 1 OhPa, 
and test for skewness and departures from Gaussian. The 
reason the issue of multimodality in the climate system is 
considered important is because it could lead to nonlinear- 
ity in the response to an external forcing [ Palmer , 1999]. 
However, we note that multimodality is not a requirement 
for such a response, and a nonlinear stratospheric response 
could lead to a nonlinear tropospheric response even if tro- 
pospheric distributions are Gaussian. 

Histograms of temperature at 30hPa over the North 
Pole for each month were shown by Labitzke [1982]. As 
well as a seasonal cycle in the mean, the winter months 
having the lowest temperatures, these PDFs exhibit a sea- 
sonal cycle in variability, with a variability maximum in 
February, and much less variability in the summer months. 
This is due to the fact that in summer, prevailing east- 
erly winds in the stratosphere prohibit the upward prop- 
agation of planetary waves from the troposphere [ Char - 
ney and Drazin , 1961]. The PDFs also exhibit consider- 
able positive skewness in December and January. Labitzke 
[1977] describes how warmings of the polar stratosphere, 
corresponding to a weakened vortex, are generally asso- 
ciated with a strengthening of the wavenumber one dis- 
turbance to the geopotential height field. Canadian warm- 
ings in particular are associated with a strengthening of the 
Aleutian High and a displacement of the vortex towards 
Siberia. We examine the shape of the vortex when the AO 
index is particularly high and low in this study, to look for 
evidence of such phenomena and to find physical explana- 
tions for the shape of the PDFs. 

The Arctic Oscillation index at the surface has recently 
exhibited an upward trend which is statistically significant 
when compared to a red noise model [Thompson et ai , 
2000], and when compared to control variability in a GCM 
[Gillett et ai , 2000]. There has been an accompanying 
increase in the strength of the stratospheric polar vortex 
[e.g. Graf et ai , 1998; Thompson et ai, 2000], which is 
consistent with a cooling of the polar stratosphere due to 
ozone depletion and greenhouse gases [Perlwitz and Graf, 
1995; Shindell et ai, 1998], These results are supported 
by several GCM modelling studies, which have shown 
that anthropogenic forcing causes an increase in the AO 
index at the surface [Shindell et ai, 1999; Fyfe et ai, 

1 999; Paeth et ai , 1 999] . Corti et ai [ 1 999] note a change 
in the shape of the distribution of the first two principal 
components of 500hPa geopotential height over the past 
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forty years, consistent with a ‘regime transition view of 
climate change [Palmer, 1999]. We look for evidence of 
such a change in the distribution of the AO index in both 
the troposphere and stratosphere, which could signify a 
nonlinear response to external forcing. 

2. Downward Propagation of the Arctic 
Oscillation 

We followed the method of Baldwin and Dunkerton 
[1999] to derive AO indices on a range of levels from 1000 
to lOhPa by means of an EOF analysis of NCEP geopo- 
tential heights (see Appendix for more details). Baldwin 
and Dunkerton [1999] calculate the 20-day lag correlation 
between the 90-day low-pass filtered Arctic Oscillation in- 
dex at lOhPa and successive levels downwards to the sur- 
face, and note an apparent downward propagation. In or- 
der to assess the significance of these results we calculated 
the lagged cross correlation between the unfiltered 10 and 
lOOOhPa indices, and compared our results with a two- 
channel AR( 1) model, constructed as follows: 

x t = Ax t -i + Qz t (1) 

where x t is a 2-element vector representing the AO in- 
dex at the 10 and lOOOhPa levels at time t, A is a 2 x 2 
matrix of AR coefficients, Q is a symmetric 2x2 ma- 
trix, and z t is a 2-element vector of uncorrelated unit vari- 
ance Gaussian noise [ Priestley , 1981] (equivalent to a POP 
model [von Storch et al , 1995]). By post-multiplying 
by xf_ x and taking expectations we obtain the estimator 
A = R(1)R(0) _1 where E[ A] = A. R(s) consists 
of the auto- and cross-covariances of the two indices at 
lag s, estimated from the observed AO indices. By post- 
multiplying by xj and taking expectations, we find that 
QQ t = R(0) - R(1)R(0) -1 R(1) 7 \ If we are to avoid 
building the downward propagation into our noise model, 
since it is its existence that we are to test, then we must 
set the off-diagonal elements of A to be equal. This was 
accomplished here by setting both to the mean of the two. 
Having evaluated A and Q, we generated a 1000-member 
Monte Carlo ensemble using this model. This method 
generates an ensemble of pairs of red noise series, each 
with the same variance and lag-1 autocorrelation as the 10 
and lOOOhPa indices respectively, and with the same lag-0 
cross-correlation. 

In order to use our noise model to assess significance, 
it was important to ensure that it was spectrally realis- 
tic. With daily sampling there were wide discrepancies 
between the AO spectra and red noise, but these discrepan- 
cies were reduced as the sampling interval was increased 


at the expense of a reduction in sample size. We sub- 
sampled at 5 day intervals, which gave reasonable agree- 
ment with a red noise spectrum for those periods important 
for our hypothesis testing of less than about two months 
(figures la and lb), but our results regarding downward 
propagation were insensitive to the exact sampling inter- 
val chosen. A comparison of the observed spectra with 
those of the noise model (figures la and lb) shows re- 
markably good agreement at lOOOhPa at periods up to two 
years, whilst the lOhPa index appears to follow a power 
law other than that of an AR(1) model. This is interest- 
ing, but not of direct relevance to the present investigation, 
and will be pursued elsewhere. The discrepancy between 
the AO index at lOOOhPa and the AR(1) model at peri- 
ods longer than two years is consistent with Thompson et 
al. [2000] ’s findings that the observed trend in the surface 
AO is significant compared with a red noise model. 

The cross-correlation function of the 10 and lOOOhPa 
indices is shown in figure lc. This function is asymmetric, 
having much larger amplitude at negative tag, which cor- 
responds to the lOhPa index leading the lOOOhPa index, 
or downward propagation. The grey band shown in figure 
lc represents the range between the 2.5th and 97.5th per- 
centiles of our noise model. It may clearly be seen that the 
observed lagged cross-correlation is significantly greater 
than that of the noise model for lags of around 15-40 
days, peaking at 20 days, which agrees closely with the 
time for propagation from lOhPa to the surface of around 
three weeks noted by Baldwin and Dunkerton [1999]. In 
order to form a parametric estimate of the significance ol 
the result, we calculated the ratio of the area under the 
cross-correlation function at positive lags to the area at 
negative lags, and compared with the noise model. Us- 
ing this test the downward propagation result is significant 
at the 99% level. 

3. PDFs of daily Arctic Oscillation Indices 

Spectra and cross-correlations only provide informa- 
tion on the first and second moments of the distributions, 
so in order to look in more detail at the statistics ot the 
AO indices, we plotted histograms for each month on 
four levels (figure 2). The Arctic Oscillation is derived 
based on anomalies from the seasonal cycle, so in each 
case the histograms are centred on zero. However, a clear 
seasonal cycle in variance may be seen at all levels. In 
the troposphere the seasonal cycle in variance is relatively 
smooth, with a minimum in July and a maximum in Jan- 
uary. In the stratosphere, by contrast, there is an abrupt 
transition from the variable winter regime to the summer 
regime with very little variance, which generally occurs 
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in April and is known as the final warming [Godson and 
Lee , 1958]. The low variance observed in summer occurs 
because prevailing easterly winds in the stratosphere pro- 
hibit the upward propagation of planetary waves from the 
troposphere [Charney and Drazin, 1961]. Generally our 
histograms appear similar in form to those of the North 
Pole 30hPa temperature shown by Labitzke [1982]. This 
is to be expected since the AO index is anti-correlated with 
polar stratospheric temperatures [Thompson and Wallace , 
2000 ]. 

Several authors have suggested that the first principal 
component of Northern Hemisphere geopotential height 
may be multimodal at 500hPa [Hansen and Sutera , 1992; 
Corti et ai, 1999], even when only monthly data are used. 
However, we find no significant multimodality in these AO 
indices on any level, based on daily data (apparent multi- 
ple maxima in our histograms result from the limited sam- 
ple, and arc not significant). In the troposphere, the PDFs 
arc generally not significantly skewed (when tested against 
a noise model following the approach in section 2), and 
departures from Gaussian are small. In the stratosphere, 
however, the PDFs exhibit negative skewness in Decem- 
ber and January and positive skewness in April. The rea- 
sons for these departures from normality are examined in 
the following sections. Departures from Gaussian were 
generally found to be largest at 70hPa (figure 2), and for 
this reason we show results for this level in the following 
analysis. However the main conclusions hold over all the 
stratospheric levels considered from 150 to lOhPa. 

3.1. April results 

The AO index in April exhibits strong positive skew- 
ness in the stratosphere, most notably at 70hPa (figure 
2). From figure 3 we can conclude that the final warm- 
ing is most likely to occur in April. The final warming is 
associated with a transition to easterly mean flow, which 
prohibits the upward propagation of planetary waves, dra- 
matically reducing stratospheric variability, and render- 
ing summer and winter states physically distinct [Charney 
and DrazifU 1961]. Figure 3 shows the large change in 
variance which is associated with the final warming, and 
which usually occurs in April. We thus conclude that it 
is the superposition of PDFs corresponding to the summer 
and winter states, which we expect to be distinct on phys- 
ical grounds, which leads to the positive skewness of the 
AO PDF in April. 

3.2. January results 

The AO index in January (and less markedly in De- 
cember) exhibits negative skewness in the stratosphere, as 


may be clearly seen in figure 2. However, an examination 
of figure 3 shows that, unlike in April, this non-normality 
is not due to a one-way transition between two distinct 
states: the PDF changes little through the month. Since 
temporal variations in the circulation offered no insight 
into the causes of the skewness, we turned to its spatial 
characteristics, by examining the form taken by the vor- 
tex on days at the extremes of the distribution. We took 
the 5% of the days with the lowest 70hPa index for this 
month, and the 5% of the days with the highest index (as 
indicated by the shading in figure 4a), and plotted a con- 
tour of the 70hPa surface for each day in both cases (fig- 
ures 4b and 4c respectively). We show the full flow rather 
than anomalies to aid physical interpretation. It is apparent 
that on the days with large positive AO index the vortex is 
centred close to the pole, and the pattern of geopotential 
height is almost symmetric, as might be expected. By con- 
trast, when the index is negative, and the vortex is warm, 
its centre is generally displaced south towards Eurasia, but 
there is large variation between individual days. This shift 
of the vortex towards Eurasia which occurs during sudden 
warmings is consistent with a strengthening of the Aleu- 
tian High as discussed by e.g. Labitzke [1977]; Baldwin 
and Holton [1988]. 

The spatial patterns of the strong and weak vortex states 
do not in themselves explain why the PDFs are skewed, 
but the difference in variability within the samples may 
help to do so. Whilst the weak vortex states are very dif- 
ferent to one another, the strong vortex states in January 
are all very similar. For comparison, similar contour plots 
were produced for the strong and weak vortex states in 
April (figures 5b and 5c). A different contour height is 
used because the mean vortex state is weaker. It is notable 
that whilst the general form of the strong and weak vortex 
flows is similar in both months, the strong vortex states in 
January are much more similar to each other than those in 
April. We might thus hypothesise a physical limit to the 
strength of the vortex in January, or correspondingly to 
polar stratospheric temperature, since this is closely cor- 
related with the mode [Thompson and Wallace , 2000]. In 
principle, such a limit exists: in the absence of dynam- 
ical heating mainly by planetary waves, the winter polar 
stratosphere would cool to its radiatively determined state 
[Shine, 1987]. Based on monthly data, Shine [1987] con- 
cludes that mean stratospheric temperatures around the 
North Pole are not strongly influenced by this radiative 
limit, in contrast to temperatures near the South Pole. 
However although the mean temperatures are not close to 
the radiatively determined limit, those on individual days 
may be. We tested this hypothesis by comparing daily 
70hPa North Pole temperatures from the NCEP reanaly- 
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sis with radiatively determined temperatures. These differ 
from radiative equilibrium temperatures, in that they take 
account of the seasonal cycle. The difference is significant 
in the winter stratosphere due to the long radiative cool- 
ing timescales associated with low temperatures. We use 
radiatively determined temperatures updated from Shine 
[1987] using the ozone climatology described in Shine 
[1989] and other small improvements. Note that calcu- 
lated temperatures depend on the set-up of the model used, 
for example, the position of the lower boundary and the 
fact that the ozone climatology is that of a radiatively de- 
termined state. Thus, the radiatively determined tempera- 
tures used here should be regarded as approximate. 

The NCEP polar temperatures at 70hPa show a min- 
imum in January, and as expected, the temperatures are 
strongly anti-correlated with the AO index on the same 
level (r=-0.87 in January, and r=-0.73 in April). His- 
tograms of the 70hPa polar temperature thus show the op- 
posite skewness to those of the AO index (figure 6), and 
the apparent cutoff at low temperatures in January is even 
more marked than the corresponding cutoff at high AO 
index values. The radiatively determined temperature for 
the middle of the month, interpolated to 70hPa at the North 
Pole is shown in each case for comparison. These are of 
course lower than the actual temperatures observed, since 
dynamical effects act to heat the polar winter stratosphere. 
However, in January (figure 6b) the lowest temperatures 
are only about 5K warmer than the radiatively determined 
temperature, which is small compared to the total range of 
around 50K. Note that at the onset of the vortex in Octo- 
ber (figure 6a), when the stratosphere is relatively undis- 
turbed, the mean 70hPa temperature is over 10K warmer 
than the radiatively determined temperature. Thus given 
the long radiative timescale, we would expect the strato- 
sphere to be considerably warmer in January than the ra- 
diatively determined state even if it remained undisturbed 
from October. Further, the mean circulation leads to some 
dynamic heating even in the absence of planetary waves, 
which also prevents the limit from being reached. 

A similar comparison of the 70hPa North Pole temper- 
atures was carried out for April (figure 6c). Here the ra- 
diatively determined temperature is much colder than the 
mean observed temperature (38K versus 2 IK in January), 
so it causes no corresponding cutoff in the PDF. Thus the 
radiatively determined state is indeed likely to act to limit 
negative fluctuations in stratospheric temperature or pos- 
itive AO anomalies, and is likely to be the cause of the 
observed skewness in January. We also found that the 
mean observed temperature in January at 70hPa is closer 
to the radiatively determined temperature than at any of 
the higher levels we considered (the model used to com- 


pute radiatively determined temperatures is not applicable 
to lower levels). This is consistent with the cutoff of AO 
and temperature PDFs being most pronounced here. 

An obvious extension of this work would be to apply 
the same analysis to Southern Hemisphere geopotential 
heights, where conditions are much closer to the radia- 
tively determined state during the winter [Shine, 1987]. 
However, the Southern Hemisphere annular mode exhibits 
much less coupling between the troposphere and strato- 
sphere in this season, and hence the EOF-analysis method 
used in the Northern Hemisphere cannot be applied di- 
rectly. Further, data are only available from 1979, mak- 
ing it difficult to identify changes in the distribution of the 
indices. However, other approaches do allow an annular 
mode to be identified, and its properties will be discussed 
in future work. 

4. Changes in the Arctic Oscillation PDFs 

As mentioned in section 1 , over the past thirty years the 
Northern Hemisphere stratosphere has exhibited a cooling 
trend over the pole, accompanied by an increase in the 
strength of the vortex. In order to examine this effect in 
our data we computed the difference in the means of the 
AO index for each month over the periods 1958-1 977 and 
1978-1997. We assessed the significance of our results by 
comparison with a Monte Carlo ensemble of 1500 series 
generated with an AR(1) model. Noise parameters were 
estimated from the 5-day subsampled AO index. The dif- 
ference in mean AO index is plotted where significant in 
figure 7a. It may be seen that changes in the AO index 
are positive and largest in January both near the surface 
and in the stratosphere. Smaller changes are also found in 
November, December and August. Based on thirty years 
of monthly SLP data, Thompson et ai [2000] also noted 
positive trends in these months, but attributed significance 
to the trends in January, February, March and August. 

Greenhouse gas increases and ozone depletion together 
act to cool the stratosphere radiatively. The former ef- 
fect is dominant in winter, but is only expected to account 
for a change in the radiatively determined state of around 
0.5K/decade in the region we consider, leading to a differ- 
ence between the means of the two halves of the record 
of order IK (K. Shine, pers. comm.). However a second 
mechanism by which vortex strength is effected by an- 
thropogenic forcing is discussed by Shindell et ai [1999]. 
Cooling of the polar stratosphere and warming of the 
tropical upper troposphere causes an increased meridional 
temperature gradient in the region of the tropopause, and 
an associated change in the vertical shear of the zonal 
wind. This change tends to deflect high-latitude plane- 
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tary waves southwards, reducing their heating effect on 
the polar stratosphere, and reducing the number of sudden 
warmings. If this mechanism has been important in deter- 
mining the response of the stratospheric circulation to an- 
thropogenic forcing, we might expect the number of warm 
vortex days in the “tail” to decrease, while the peak at pos- 
itive vortex states stays in much the same place, changing 
the shape of the PDF. Owing to the fact that in January the 
vortex is coldest and the negative skewness is largest, we 
would expect this effect to be strongest in this month. Of 
the levels on which the radiatively determined temperature 
was available, 70hPa had temperatures closest to this limit, 
so changes in the shape of the PDF might be expected to 
be largest here. We thus plot histograms of the daily 70hPa 
AO index and North Pole temperature in January for each 
half of the record, 1958-1977 and 1978-1997 (figure 8). 
Both histograms show that the number of weak vortex 
events (sudden warmings) is lower in the second half of 
the record than in the first, whilst the broad maximum of 
the distribution has moved relatively little. The limiting 
effect is clearest in the polar temperatures, which are most 
directly limited by the radiatively determined state. Such 
a change in the distribution of stratospheric polar temper- 
atures in January has since been noted in a stratosphere- 
resolving GCM in response to stratospheric ozone deple- 
tion (S. Rossier, pers. comm.), which has a similar ef- 
fect on planetary wave propagation to a greenhouse gas 
increase [Shindell et al , 2000]. 

We tested the significance of the change in the shape of 
the AO PDFs by subtracting the mean AO index of each 
half, and computing the Kolmogorov-Smirnov statistic for 
the two distributions (figure 7b). This statistic is just the 
maximum difference between the two cumulative distri- 
butions. We evaluated the significance of these results by 
comparison with our Monte-Carlo noise ensemble. Our 
results indicate that the change in the shape of the PDF in 
January is significant against red noise between 200 and 
20hPa at the 95% level. Thus the shape of the PDF of the 
Arctic Oscillation index is changing in January as well as 
its mean. This is consistent with the hypothesis that the 
main response to forcing has been through the southward 
deflection o( planetary waves, which has reduced the num- 
ber of sudden warmings, whilst the strong vortex states, 
limited by the slowly changing radiatively determined 
state, have changed little. We might thus eventually ex- 
pect the former effect to saturate, so that no sudden warm- 
ings occur in January, leaving the vortex strength corre- 
spondingly less sensitive to further increases in forcing. 
Such a saturation of response has been noted in a GCM by 
Shindell et ai [2000]. This could explain the discrepancy 
between the observed AO sensitivity to forcing and that 


of GCMs with poor stratospheric resolution [Gillett et al . , 
2000], whose unperturbed winter stratospheres are often 
already close to the radiatively-determined state, an issue 
we pursue elsewhere. 

5. Conclusion 

While we find PDFs of the AO index to be Gaussian- 
like in the troposphere, in contrast to the results of Corti 
et ai [1999], we note that nonlinearity in the stratospheric 
response to forcing could make the tropospheric response 
nonlinear through its downward influence. This hypothe- 
sis is supported by our finding that the downward propa- 
gation of AO anomalies in the NCEP reanalysis noted by 
Baldwin and Dunkerton [1999] is statistically significant, 
and that the AO exhibits significant departures from nor- 
mality in the stratosphere. In April PDFs of the AO index 
in the stratosphere are positively skewed, with a narrow 
variance peak corresponding to days on which the final 
warming has occurred, and a tail to positive AO index, cor- 
responding to days when the winter vortex persists. In Jan- 
uary, when the polar stratosphere is dynamically heated by 
planetary wave activity, the observed negative skewness is 
associated with the limit on vortex strength correspond- 
ing to the radiatively determined state. Whilst this state 
changes only slowly in response to anthropogenic forc- 
ing, the southward deflection of planetary waves by the 
changed vertical shear of zonal wind could lead to a large 
reduction in the number of warming events, changing the 
distribution of the index. Using a KS-test and a Monte- 
Carlo ensemble of red noise series, we found that there 
has indeed been a significant change in the shape of the 
PDF between 200 and 20hPa in January between the pe- 
riods 1958-1977 and 1978-1997, consistent with our hy- 
pothesis. As the planetary wave deflection effect saturates 
(when there are no more sudden warmings), we would ex- 
pect a change in the sensitivity of the mean vortex strength 
to anthropogenic forcing. In view of the downward influ- 
ence we observe, we speculate that as greenhouse gas con- 
centrations further increase, the tropospheric Arctic Oscil- 
lation may also show a reduced sensitivity to the forcing, 
with implications for surface climate change. 

Appendix: Calculation of AO Indices 

We use the National Center for Environmental Predic- 
tion/National Center for Atmospheric Research (NCEP/NCAR) 
Reanalysis [Kalnay et al. , 1996] geopotential heights and 
temperatures, available as daily averages on a 2.5° x 2.5° 
grid on 17 levels from 1000 to lOhPa for the years 1958- 
1997. By comparison with other data sources, J. Perlwitz 
(pers. comm.) concluded that the stratospheric data are un- 
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reliable in summer and in the tropics, but that northward of 
20°N they are reliable for the study of large-scale circula- 
tion patterns in the lower stratosphere. However, we note 
that changes in the observations assimilated around 1979 
may influence higher moments of the AO indices, and it 
would be prudent to verify these results using other re- 
analyses as they become available. Following the method 
of Baldwin and Dunkerton [1999], an EOF analysis was 
first performed on 90-day low-pass filtered DJF geopoten- 
tial height anomalies on five levels ( 1 000, 300, 1 00, 30 and 
lOhPa) over the region northward of 20°N. The data on 
each level were normalised separately by dividing by the 
standard deviation, and an equal-area grid was used. Spa- 
tial patterns for each of the 17 levels were then recovered 
by regressing geopotential height data onto the first prin- 
cipal component from the EOF analysis. Indices at each 
level were obtained by regressing the geopotential heights 
onto each of these spatial patterns. These are referred to 
by Baldwin and Dunkerton [1999] as ‘AO signature time 
series’, but we use the term 4 AO indices’ tor simplicity. 

References 

Baldwin, M. P. and T. J. Dunkerton, Propagation of the 
Arctic Oscillation from the stratosphere to the tropo- 
sphere,/ Geophys . Res., 104 , 30937-30946, 1999. 
Baldwin, M. P. and J. R. Holton, Climatology of the strato- 
spheric polar vortex and planetary wave breaking, / 
Atmos . Sci ., 45 , 1 123-1 142, 1988. 

Charney, J. G. and P. G. Drazin, Propagation of planetary- 
scale disturbances from the lower into the upper atmo- 
sphere, / Geophys. Res., 66, 83-109, 1961 . 

Cheng, X. and J. M. Wallace, Cluster analysis of the 
Northern Hemisphere wintertime 500-hpa height field: 
Spatial patterns, / Atmos. Sci., 50, 2674-2696, 1993. 
Corti, S., F. Molteni, and T. N. Palmer, Signature of recent 
climate change in frequencies of natural atmospheric 
circulation regimes, Nature , 398 , 799-802, 1999. 

Fyfe, J. C., G. J. Boer, and G. M. Flato, The Arctic and 
Antarctic Oscillations and their projected changes un- 
der global warming, Geophys. Res. Lett , 26, 1 60 1 — 
1604, 1999. 

Gillett, N. P., G. C. Hegerl, M. R. Allen, and P. A. Stott, 
Implications of changes in the Northern Hemisphere 
circulation for the detection of anthropogenic climate 
change, Geophys. Res. Lett, 27, 993—996, 2000. 
Godson, W. L. and R. L. Lee, High level fields of wind 
and temperature over the Canadian Arctic, Beitr. Phys. 
Atmos., 31 , 40-68, 1958. 

Graf, H. F., I. Kirchner, and J. Perlwitz, Changing lower 
stratospheric circulation: The role of ozone and green- 


house gases, / Geophys. Res., 103, 11251 — 11261, 
1998. 

Hansen, A. R. and A. Sutera, Structure in the phase space 
of a general circulation model deduced from empiri- 
cal orthogonal functions, / Atmos. Sci., 49, 320-326, 

1992. 

Kalnay, E. et al„ The NCEP/NCAR 40-year reanalysis 
project, Bull Amer. Met. Soc., 77, 437-471, 1996. 

Kimoto, M. and M. Ghil, Multiple flow regimes in the 
Northern Hemisphere winter. Part I: Methodology and 
hemispheric regimes, / Atmos. Sci., 50, 2625-2643, 

1993. 

Kuroda, Y. and K. Kodera, Role of planetary waves in 
the stratosphere-troposphere coupled variability in the 
Northern Hemisphere winter, Geophys. Res. Lett, 26, 
2375-2378, 1999. 

Labitzke, K., Interannual variability of the winter strato- 
sphere in the Northern Hemisphere, Mon. Wea. Rev., 
105, 762-770, 1977. 

Labitzke, K., On the interannual variability of the mid- 
dle stratosphere during the Northern winters, / Meteor. 
Soc. Japan, 60, 124-138, 1982. 

Murray, F. W., Dynamic stability in the stratosphere, / 
Geophys. Res., 65, 3273-3305, 1960. 

Nakamura, H. and J. M. Wallace, Skewness of the low- 
frequency fluctuations in the tropospheric circulation 
during the Northen Hemisphere winter, / Atmos. Sci., 
48, 1441-1448, 1991. 

Paeth, H., A. Hense, R. Glowienka-Hense, S. Voss, and 
U. Cubasch, The North Atlantic Oscillation as an in- 
dicator for greenhouse-gas induced regional climate 
change, Climate Dynamics , 15, 953-960, 1999. 

Palmer, T. N., A nonlinear dynamical perspective on cli- 
mate prediction, / Climate, 72,575-591, 1999. 

Perlwitz, J. and H. F. Graf, The statistical connection be- 
tween tropospheric and stratospheric circulation of the 
Northern Hemisphere in winter, / Climate , 8, 228 1— 
2295, 1995. 

Priestley, M. B., Spectral analysis and time series , vol- 
ume 2, 685-687 pp.. Academic Press, London, 1981. 

Shindell, D. T., R. L. Miller, G. Schmidt, and L. Pandolfo, 
Simulation of the Arctic Oscillation trend by green- 
house forcing of a stratospheric model, Nature, 399, 
452^55, 1999. 

Shindell, D. T., D. Rind, and P. Lonergan, Increased polar 
stratospheric ozone losses and delayed eventual recov- 
ery owing to increasing greenhouse-gas concentrations, 
Nature, 392, 589-592, 1998. 


Shindell, D. T., G. A. Schmidt, R. L. Miller, and D. Rind, 
Northern hemisphere winter climate response to green- 
house gas, ozone, solar and volcanic forcing, J. Geo - 
phys. Res ., 2000. submitted. 

Shine, K. R, The middle atmosphere in the absence of dy- 
namical heat fluxes. Quart. J. Roy. Met. Soc 113, 603- 
663, 1987. 

Shine, K. R, Sources and sinks of zonal momentum in the 
middle atmosphere diagnosed using the diabatic circu- 
lation, Quart. J. Roy. Met. Soc., 1 15, 265-292, 1989. 
Thompson, D. W. J. and J. M. Wallace, The Arctic Oscilla- 
tion signature in the wintertime geopotential height and 
temperature fields, Geophys . Res. Lett , 25, 1297-1300, 
1998. 

Thompson, D. W. J. and J. M. Wallace, Annular modes 
in the extratropical circulation Part I: Month to month 
variability,/ Climate, 5, 1000-1016,2000. 

Thompson, D. W. J., J. M. Wallace, and G. C. Hegerl, 
Annular modes in the extratropical circulation Part II: 
Trends, / Climate, 5, 1018-1036,2000. 
von Storch, H., G. Burger, R. Schnur, and J.-S. von Storch, 
Principal oscillation patterns: A review, / Climate , 8, 
377-400, 1995. 

Wallace, J. M., North Atlantic Oscillation / Annular 
Mode: Two paradigms - One phenomenon, Quart. J. 
Roy. Met. Soc., 2000. in press. 

Yamazaki, K. and Y. Shinya, Analysis of the Arctic Oscil- 
lation simulated by a GCM, / Meteor. Soc. Japan, 77, 
1287-1298, 1999. 

M. R. Allen, Space Science and Technology Depart- 
ment, Rutherford Appleton Laboratory, Chilton, Didcot 
OX1 1 0QZ, UK. (e-mail: m.r.allen@rl.ac.uk) 

M. P. Baldwin, Northwest Research Associates, 14508 
NE 20th Street, Bellevue, WA 98007-3713. (e-mail: 
mark@nwra.com) 

N. P. Gillett, Atmospheric, Oceanic and Planetary 
Physics, Clarendon Laboratory, Parks Road, Oxford OX1 
3PU, UK. (e-mail: giIlett@atm.ox.ac.uk) 


This preprint was prepared with AGU’s iATgX macros v5.01. Fite 
aochange pp formatted May 25, 2(XX). 


9 


(a) 





Figure 1 . Power spectra of 5-day subsampled AO index at (a) lOOOhPa, (b) lOhPa, and (c) their cross-correlation function. 
In each case the range between the 2.5th and 97.5th percentiles of the bivariate AR(1) model is shown by a grey band. In c 
negative lag corresponds to the lOhPa index leading the lOOOhPa index. 
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Figure 3. 70hPa AO index for each year 1 958- 1 997 plotted against day number. Note that the index is based on deseasonalised 
geopotential heights. Skewness in January and April is visible as clustering above and below zero respectively. 
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(a) 




Figure 4. a is an enlarged version of the histogram of AO index at 70hPa in January shown in figure 2. The light shaded 
region represents the 5% of days with the lowest AO index, and the 17.75km contour of 70hPa geopotential height on all 
these days is plotted in panel b. The dark shaded region corresponds to the 5% of days with the highest AO index, and the 
associated 17.75km contours of 70hPa height are shown in c. 
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(a) 




(c) 



Figure 5. As figure 4, but for April. The 1 8.2km contour of 70hPa geopotential height is shown in panels b and c. 
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(a) October 



(b) January 



(c) April 



Figure 6. Histograms of daily NCEP 70hPa North Pole temperature over the period 1958-1997 in (a) October, (b) January, 
and (c) April. The corresponding radiatively determined temperatures are shown for comparison. 
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(a) 



Month 

(b) 


KS statistic 



Month 


Figure 7. a shows the change in the mean AO index between the periods 1958-1977 and 1978-1997, only plotted where 
<95% or >95% of those in the noise model, b shows the KS-statistic, a measure of the difference between the AO index 
distributions of the two periods, independent of the change in the mean, and only plotted where significant at the 95% level. 






Figure 8. Histograms of the daily January 70hPa AO index 
and 1978-1997 (grey line). 
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ABSTRACT 

The quasi-biennial oscillation (QBO) dominates the vari- 
ability of the equatorial stratosphere (-16-50 km) and is 
most easily seen as downward propagating easterly and 
westerly wind regimes, with a variable period averaging 
approximately 28 months. From a tluid dynamical perspec- 
tive, the QBO is a fascinating (and perhaps the premier) 
example of a coherent, oscillating mean flow that is driven 
by vertically propagating waves with periods unrelated to 
that of the resulting oscillation. Although the QBO is a 
tropical phenomenon, it affects the stratospheric flow from 
pole to pole by modulating the effects of extratropical 
waves. The QBO affects variability in the mesosphere near 
85 km by selectively filtering waves that propagate upward 
through the equatorial stratosphere, and even appears to 
affect the strength of Atlantic hurricanes. Indeed, study of 
the QBO is inseparable from the study of atmospheric wave 
motions that drive it and are modulated by it. The effects of 
the QBO are not confined to atmospheric dynamics. Chemi- 
cal constituents, such as ozone, water vapor, and methane, 
are affected by circulation changes, in both hemispheres, 
induced by the QBO. There are also substantial QBO sig- 
nals in many of the shorter-lived chemical constituents. 
Through modulation of extratropical wave propagation, the 
QBO has an effect on the breakdown of the wintertime 
stratospheric polar vortices and severity of high-latitude 
ozone depletion. The polar vortex in the stratosphere tends 
to affect surface weather patterns, providing a mechanism 
for the QBO to have a small effect at the Earth’s surface. As 
more data sources (e.g., wind and temperature measure- 
ments from both ground-based systems and satellites) be- 
come available, the effects of the QBO can be more pre- 
cisely assessed. 

This review covers the current state of knowledge of the 
tropical QBO, its extratropical dynamical effects, chemical 


constituent transport, and effects of the QBO in the tropo- 
sphere and above the stratosphere. It is intended to provide a 
broad overview of the QBO and its effects to researchers 
outside the field, as well as a source of information and ref- 
erences for specialists. The history of research on the QBO 
is discussed only briefly, and the reader is referred to several 
historical review papers. The basic theory of the QBO is 
summarized, and tutorial references are provided. 
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1. INTRODUCTION 

1.1. The discovery of the Quasi-Biennial 
Oscillation (QBO) 

The first observations of equatorial stratospheric winds 
were made when it was discovered that debris from the 
eruption of Krakatau (1883) circled the globe from east to 
west in a period of about two weeks; these winds became 
known as the “Krakatau easterlies.” (For a more complete 
review of the discovery of the QBO, and advances in obser- 
vations and theory, the reader is referred to Maruyamci 
[1997] and Labitzke and van Loon [1999]). In 1908 the 
German meteorologist Berson launched balloons from 
tropical Africa and found winds blowing from west to east 
at about 15 km, near the tropopause, which became known 
as the “Berson westerlies.” For nearly 50 years there were 
only sporadic balloon observations [ Hamilton , 1998a] to 
contradict the existence of equatorial stratospheric easterly 
winds overlying westerly winds. Palmer [1954] used upper 
air sounding data, gathered to study nuclear testing at the 
Marshall Islands, to find that the transition between the 
“Berson westerlies” and “Krakatau easterlies” varied from 
month to month and year to year. However, the data were 
insufficient to show any periodicity. By using observations 
from Christmas Island (2.0°N), Graystone [1959] contoured 
two years of wind speeds in the time-height plane, which 
showed gradually descending easterly (italicized terms are 
defined in the Glossary) and westerly wind regimes. (Note: 
The meteorological term “easterlies” describes winds that 
blow from the east, while “westerlies” blow from the west. 
Unfortunately, the terms “ eastward ’ (to the east) and “ west- 
ward ’ (to the west) are used commonly in discussions of 
wave propagation and llow, and have the opposite meaning 
of easterly and westerly. Throughout this review we will use 
both descriptions, as found in the literature.) 

Discovery of the QBO must be credited to the independ- 
ent work of R.J. Reed in the United States and R.A. Ebdon 
in Great Britain. In a paper entitled “The Circulation of the 
Stratosphere” presented at the 40th anniversary meeting of 
the American Meteorology Society, Boston, January 1960, 
Reed announced the discovery, using Canton Island data 
(2.8°S) of “alternate bands of easterly and westerly winds 
which originate above 30 km and which move downward 
through the stratosphere at a speed of about 1 km per 
month.” Reed pointed out that the bands “appear at intervals 
of roughly 13 months, 26 months being required for a com- 
plete cycle.” Ebdon [1960] also used data from Canton Is- 
land, spanning 1954-1959, to show that the wind oscillation 
had an apparent two-year period. Ebdon and Veryard [1961] 
used additional data from Canton Island (January 1954- 
January 1960) at 50 hPa (hectopascals, equal to millbars) to 
show that the wind fluctuated with a period of 25-27 
months, rather than exactly two years. They extended the 
earlier study of Ebdon to include other equatorial stations 


and concluded that the wind fluctuations occurred simulta- 
neously around the equatorial belt, and estimated that the 
wind regimes (when the equatorial winds are easterly or 
westerly) took about a year to descend from 10 hPa to 60 
hPa. Veryard and Ebdon [1961] extended this study to find a 
dominant period of 26 months and observed similar fluctua- 
tions in temperature. 

With the observation of a longer-period cycle starting in 
1963, Angell and Korshover [1964] coined the term “quasi- 
biennial oscillation,” which gained acceptance and was ab- 
breviated “QBO.” Many of the dynamical aspects of the 
QBO are best illustrated by a time-height cross-section of 
monthly-mean equatorial zonal (longitudinal) wind (Figure 
1). Ideally, such a diagram would be of the zonally-averaged 
zonal wind, but the approximate longitudinal symmetry of 
the QBO allows rawinsonde (balloon) observations from a 
single station near the equator to be sufficient. The alternat- 
ing wind regimes repeat at intervals that vary from 22 to 34 
months, with an average period of slightly more than 28 
months. Westerly shear zones (in which westerly winds in- 
crease with height) descend more regularly and rapidly than 
easterly shear zones. The amplitude of -20 ms' 1 is nearly 
constant from 5 hPa to 40 hPa, but decreases rapidly as the 
wind regimes descend below 50 hPa. Rawinsonde observa- 
tions are used up to 10 hPa, and rocketsondes (meteorologi- 
cal rockets) above. The QBO amplitude diminishes to less 
than 5 ms’ 1 at 1 hPa, near the stratopause. The amplitude of 
the QBO is approximately Gaussian about the equator with 
a 12 degree half width and little phase dependence on lati- 
tude within the tropics [Wallace, 1973]. Although the QBO 
is definitely not a biennial oscillation, there is a tendency for 
a seasonal preference in the phase reversal [ Dunkerton , 
1990] so that, for example, the onset of both easterly and 
westerly winds regimes occurs mainly during Northern 
Hemisphere (NH) late spring at the 50-hPa level. Perhaps 
the three most remarkable features of the QBO that any the- 
ory must explain are (I) the quasi-biennial periodicity, (2) 
the occurrence of zonally symmetric westerly winds at the 
equator (conservation of angular momentum does not allow 
advection to create an equatorial westerly wind maximum), 
and (3) the downward propagation without loss of ampli- 
tude. 

1 .2. The search for an explanation of the QBO 

At the time of discovery of the QBO, there were no ob- 
servations of tropical waves and there was no theory 
predicting their existence. The search for an explanation for 
the QBO initially involved a variety of other causes: some 
internal feedback mechanism, a natural period of 
atmospheric oscillation, an external process, or some 
combination of these mechanisms. All these attempts failed 
to explain features such as the downward propagation and 
maintenance of the amplitude of the QBO as it descends 
(energy density grows as the winds descend). 
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Apparently, zonally asymmetric wave forcing is required 
to explain the equatorial westerly wind maximum. Wallace 
and Holton [1968] tried to drive the QBO in a numerical 
model through heat sources or through extratropical plane- 
tary waves propagating into the equator. They showed rather 
conclusively that lateral momentum transfer by planetary 
waves could not explain the downward propagation of the 
QBO without loss of amplitude. They made the crucial re- 
alization that the only way to reproduce the observations 
was to have a driving force (a momentum source) which 
actually propagates downward with the mean equatorial 
winds. 

Booker and Brethertons [1967] seminal paper on the ab- 
sorption of gravity waves at a critical level provided the 
spark that would lead to an understanding of how the QBO 
is driven (see Lindzen [1987] for a historical review of the 
development of the theory of the QBO). It was Lindzen s 
leap of insight to realize that gravity waves could drive a 
QBO with descending shear zones. Lindzen and Holton 
[1968] showed how a QBO could be driven by a broad spec- 
trum of vertically-propagating (propagating to the east or 
west) gravity waves. The crucial step in the understanding 
of the QBO was the realization that vertically propagating 
waves could provide the necessary momentum flux and that 
the periodicity could arise through an internal mechanism 
arising from a two-way feedback between the waves and the 
background flow. The first part of the feedback is the effect 
of the background flow on the propagation of the waves 
(and hence on the momentum fluxes). The second part of 
the feedback is the effect of the momentum fluxes on the 
background flow. Lindzen and Holton [1968] demonstrated 
explicitly, using a two dimensional (2-D) model in which 
the behavior and effect of the waves was parameterized in a 
simple way, that the feedback gave rise to an oscillation 
with a downward propagating pattern of easterly and west- 
erly winds. An important corollary to their work was that 
the period of the oscillation was controlled in part by the 
wave momentum fluxes and hence, a range of periods was 
possible. The fact that the observed oscillation had a period 
close to a subharmonic of the annual cycle was, therefore, a 
pure coincidence. 

It was a bold assertion to ascribe the forcing of the QBO 
to eastward and westward propagating gravity waves, con- 
sidering that only the westward waves had been observed. 
The theory of equatorial waves was first developed during 
the late 1960s, in parallel with the theory of the QBO. The 
solutions included a Rossby mode, and a mode which be- 
came known as the mixed R os sby-gravi ty mode. A third so- 
lution, an eastward-propagating gravity mode in the zonal 
wind only, was called the equatorial Kelvin mode. Maru- 
yama [1967] displayed observations consistent with the 
mixed Rossby-gravity mode. Wallace and Kousky [1968] 
first showed observations of equatorial Kelvin waves in the 
lower stratosphere, and noted that the wave produced an 


upward flux of westerly momentum which could account 
for the westerly acceleration associated with the QBO. 

Lindzen [1970] showed that mixed Rossby-gravity 
waves, through a mechanism involving the Coriolis force 
acting on meridional parcel displacements, contribute a net 
easterly force. Thus, the potential existed for the QBO to be 
driven by a combination of mixed Rossby-gravity waves, 
which would contribute easterly forces, and Kelvin waves, 
which would contribute westerly forces. Holton and Lindzen 
[1972] refined the Lindzen and Holton [1968] model by 
demonstrating, in a one-dimensional (1-D) model, a QBO 
driven by vertically-propagating Kelvin and mixed Rossby- 
gravity waves. The observed amplitudes of these waves, 
though small, were arguably (given the paucity of equatorial 
wave observations) large enough to drive the QBO. The 
Holton and Lindzen mechanism continued to be the ac- 
cepted paradigm for the QBO for more than two decades. 

The conceptual model of the QBO, which formed the ba- 
sis of the Lindzen and Holton [1968] model, was strongly 
supported by the famed laboratory experiment of Plumb and 
McEwan [1978], which used a salt-stratified fluid contained 
in a large annulus. The bottom boundary of the annulus con- 
sisted of a flexible membrane that oscillated up and down to 
produce vertically propagating gravity waves traveling 
clockwise and counterclockwise around the annulus. For 
waves of sufficient amplitude, a wave-induced mean flow 
regime was established that was characterized by downward 
progressing periodic reversals of the mean flow. This ex- 
periment, which remains one of the most dramatic labora- 
tory analogues of a large-scale geophysical flow, showed 
that the theoretical paradigm for the QBO was consistent 
with the behavior of a real fluid system. 

While the observed amplitudes of Kelvin and mixed 
Rossby-gravity waves are sufficient to drive a QBO in ideal- 
ized models of the atmosphere. Gray and Pyle 1 1989] found 
it necessary to increase the wave amplitudes by a lactor ol 
three to achieve a realistic QBO in a full radiative- 
dynamical-photochemical model. Dunkerton [1997] pointed 
out that the observed rate of tropical upwelling (about 1 km 
per month) effectively requires that the QBO wind regimes 
propagate downward much faster than was thought because 
the whole of the tropical stratosphere is moving upward as 
the QBO propagates downward. This fact effectively more 
than doubled the required momentum transport by vertically 
propagating equatorial waves. Observations indicate that 
mixed Rossby-gravity and Kelvin waves cannot provide 
sufficient forcing to drive the QBO with the observed peri- 
ods. Dunkerton reasoned that additional momentum flux 
must be supplied by a broad spectrum of gravity waves — 
similar to those postulated by Lindzen and Holton [1968]. 

Although the QBO is a tropical phenomenon, it influ- 
ences the global stratosphere, as first shown by Holton and 
Tan [1980]. Through modulation of winds, temperatures, 
and the circulation in the meridional plane , the QBO affects 
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(he distribution and transport of greenhouse gases and may 
be a factor in stratospheric ozone depletion. Thus, an under- 
standing of the QBO and its global effects is necessary for 
studies of long-term variability or trends in trace gases and 
aerosols. 


2. AN OVERVIEW OF THE QBO AND ITS GLOBAL 
EFFECTS 

2.1. Zonal wind 

A composite of the QBO in equatorial zonal winds (Fig- 
ure 2) [Paw son et ai , 1993] shows faster and more regular 
downward propagation of the westerly phase , and the 
stronger intensity and longer duration of the easterly phase. 
The mean period of the QBO for data during 1953-1995 is 
28.2 months, slightly longer than the 27.7 months from the 
shorter record in Naujokat [1986]. The standard deviation of 
the mean QBO is also included in Figure 2, showing 
maxima in variability close to the descending easterly and 
westerly shear zones (larger for the westerly phase). This 
mainly reflects deviations in the duration of each phase (as 
seen in Figure 1). 

Dunkerton [ 1990] showed that the QBO may be synchro- 
nized to the annual cycle, demonstrating that the onset of 
the easterly regime at 50 hPa tends to occur during NH late 
spring or summer. His analysis is updated in Figure 3, show- 
ing the onset of each wind regime at 50 hPa. The easterly 
and westerly transitions both show a strong preference to 
occur during April-June. 

The latitudinal structure of the QBO in zonal wind is 
shown in Figure 4, derived from long time series of wind 
observations at many tropical stations [Dunkerton and De- 
lisi , 1985). The amplitude of the QBO is Iatitudinally 
symmetric, and the maximum is centered over the equator, 
with a meridional half width of approximately 12°. Similar 
QBO structure is derived from assimilated meteorological 
analyses, but the amplitude is often underestimated in 
comparison to radiosonde measurements [Pawson and 
Fiorina , 1998; Randel et ai, 1999]. 

Figure 5 provides a dynamical overview of the QBO, its 
sources, and its global dynamical effects, as well as a foun- 
dation for the discussion of the details of the QBO in the 
following sections. The diagram span the troposphere, 
stratosphere, and mesosphere from pole to pole, and shows 
schematically the differences in zonal wind between easterly 
and westerly phases of the QBO. 

Convection in the tropical troposphere, ranging from the 
scale of individual clouds through mesoscale convective 
complexes (spanning more than 100 km) produces a broad 
spectrum of waves (tan wavy arrows), including gravity, 
inertia- gravity, Kelvin, and Rossby-gravity waves (see Sec- 
tion 3). These waves, with a variety of vertical and horizon- 
tal wavelengths and phase speeds, propagate into the strato- 


sphere, transporting easterly and westerly zonal momentum. 
Most of this zonal momentum is deposited at stratospheric 
levels, driving the zonal wind anomalies of the QBO. For 
each wave, the vertical profile of the zonal wind determines 
the critical level at which the momentum is deposited. The 
critical levels for these waves depend, in part, on the shear 
zones of the QBO. Some gravity waves propagate through 
the entire stratosphere and produce a QBO near the 
mesopause known as the mesospheric QBO, or MQBO 
(Section 6). 

In the tropical stratosphere the time-averaged wind 
speeds are small, so the easterly minus westerly composite 
in Figure 5 is similar in appearance to the actual winds dur- 
ing the easterly phase of the QBO. At high latitudes, there is 
a pronounced annual cycle, with strong westerly winds dur- 
ing the winter season. To the north of the equator in the 
lower stratosphere, tropical winds alter the effective 
waveguide for upward and equatorward propagating plane- 
tary-scale waves (curved purple arrows). The effect of the 
zonal wind structure in the easterly phase of the QBO is to 
focus more wave activity toward the pole, where the waves 
converge and slow the mean zonal flow. Thus, the polar 
vortex north of ~45°N shows weaker westerly winds (light 
blue). The high-latitude wind anomalies penetrate the tropo- 
sphere and provide a mechanism for the QBO to have a 
small influence on tropospheric weather patterns (Section 
6). 

2.2. Temperature and meridional circulation 

The QBO exhibits a clear signature in temperature, with 
pronounced signals in both tropics and extratropics. The 
tropical temperature QBO is in thermal wind balance with 
the vertical shear of the zonal winds, expressed for the equa- 
torial /3-plane as 


du _ — R d 2 T 

~dz.~ 7/p dy 2 


(la) 


[Andrews et ai , 1987, eq. 8.2.2], where u is the zonal 
wind, T is temperature, z is log-pressure height (approxi- 
mately corresponding to geometric altitude), y is latitude, R 
is the gas constant for dry air, H = 7km is the nominal (con- 
stant) scale height used in the log-pressure coordinates, and 
P is the latitudinal derivative of the Coriolis parameter. For 
QBO variations centered on the equator with meridional 
scale L, thermal wind balance at the equator is approximated 
as 


(lb) 


du_ -R T 
dz Hb L 2 

The equatorial temperature anomalies associated with the 
QBO in the lower stratosphere are of the order of ±4 K, 
maximizing near 30-50 hPa. Figure 6 compares time series 
(after subtraction of the seasonal cycle) of 30 hPa tempera- 
ture measurements at Singapore with the corresponding 
zonal wind vertical shear at 30 hPa, showing good correla- 
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tion (see also Figure 1). The slope of dii/dz versus tempera- 
ture estimated from regression is consistent with a merid- 
ional scale L ~ 1000-1200 km (-10° of latitude) [Randel et 
al., 1999]. 

Smaller anomalies extend downward, with QBO varia- 
tions of order of ±0.5 K observed near the tropopause [An- 
gell and Korshover, 1964], The QBO temperature anomalies 
also extend into the middle and upper stratosphere, where 
they are out of phase with the lower stratosphere anomalies. 
Figure 7 shows an example of temperature anomalies asso- 
ciated with an easterly phase of the QBO during NH winter 
1994, derived from UKMO stratospheric assimilation data 
extending to 45 km. Although these data probably underes- 
timate the magnitude of the temperature QBO (the retrievals 
average over a layer deeper than the temperature anomaly), 
the out-of-phase vertical structure is a robust feature also 
observed in long time records of satellite radiance meas- 
urements [Randel et al., 1999]. 

Besides the equatorial maximum in QBO temperature, 
satellite data reveal coherent maxima over 20-40 degrees 
latitude in each hemisphere, which are out of phase with the 
tropical signal. This is demonstrated in Figure 8 showing 
regression of stratospheric temperatures over 13-22 km 
(from the Microwave Sounding Unit Channel 4) onto the 
30-hPa QBO winds, for the period 1979-1999. One remark- 
able aspect of the extratropical temperature anomalies is that 
they are seasonally synchronized, occurring primarily dur- 
ing winter and spring in each hemisphere. Nearly identical 
signatures are observed in column ozone measurements 
(Section 5), and this seasonally-synchronized extratropical 
variability is a key and intriguing aspect of the global QBO. 
Because low frequency temperature anomalies are tightly 
coupled with variations in the mean meridional circulation, 
global circulation patterns associated with the QBO are also 
highly asymmetric at solstice (arrows in Figure 7). The tem- 
perature patterns in Figure 8 furthermore show signals in 
both polar regions, which are out ot phase with the tropics, 
and maximize in spring in each hemisphere. Although these 
polar signals are larger than the subtropical maxima, they 
are not statistically significant in this 1979-1999 record 
because of enhanced natural variability in polar regions dur- 
ing winter and spring. 

The modulation by the QBO of zonal-mean wind (Figure 
5) is coupled to modulation of the zonally averaged circula- 
tion in the latitude-height plane. The climatological circula- 
tion is characterized by large-scale ascent in the tropics, 
broad poleward transport in the stratosphere, and compen- 
sating sinking through the extratropical tropopause [Holton 
et al. , 1995]. The transport of chemical trace species into, 
within, and out of the stratosphere is the result of both large- 
scale circulations and mixing processes associated with 
waves. Chemical processes, such as those resulting in ozone 
depletion, depend not only on the concentrations of trace 
species, but may depend critically on temperature. Since the 


QBO modulates the global stratospheric circulation, includ- 
ing the polar regions, an understanding of the effects of the 
QBO on both dynamics and the distribution of trace species 
is an essential component needed to understand global cli- 
mate variability and change. 

Many long-lived trace species, such as N 2 0 and CH 4 , 
originate in the troposphere and are transported into the 
stratosphere through the tropical tropopause. Figure 9 pro- 
vides a summary of the influence of the QBO on the circula- 
tion in the latitude-height plane and the transport of chemi- 
cal trace species. In Figure 9 the contours illustrate sche- 
matically the isopleths of an idealized conservative, long- 
lived, vertically stratified tracer during the NH winter when 
the equatorial winds are easterly near 40 hPa (matching 
Figure 5). Upwelling is reflected in the broad tropical 
maximum in tracer density in the middle to upper strato- 
sphere. The extratropical influences from the QBO result in 
deviations from hemispheric symmetry. 

The heavy arrows illustrate circulation anomalies associ- 
ated with the QBO (the time-averaged circulation has been 
removed), which is here assumed to be easterly at 40 hPa. 
At the equator, the QBO induces ascent (relative to the 
mean tropical upwelling) through the tropopause, but de- 
scent in the middle to upper stratosphere. The lower strato- 
spheric ascent is nearly symmetric about the equator, while 
the circulation anomaly in the middle stratosphere is larger 
in the winter hemisphere. This asymmetry is reflected in the 
asymmetric isopleths of the tracer. In addition to advection 
by the circulation in the latitude-height plane, the tracer is 
mixed by wave motions (approximately on isentropic, or 
constant potential temperature, surfaces). This mixing is 
depicted by the wiggly horizontal arrows. The descent from 
the middle stratosphere circulation anomaly creates a ‘‘stair- 
step” pattern in the tracer isopleths between the equator and 
the subtropics (near 5 hPa). A second stair-step in the mid- 
latitude winter hemisphere is formed by isentropic mixing in 
the region of low potential vorticity gradients which sur- 
rounds the polar vortex, known as the “surf zone” [ McIntyre 
and Palmer, 1983]. Mixing may also occur equatorward of 
the subtropical jet axis in the upper stratosphere, as illus- 
trated by the wavy line near 3 hPa and 10-20°N [Dunkerton 
and O'Sullivan , 1996]. Anomalous transport from the 
Southern Hemisphere (SH) to the NH near the stratopause is 
associated with changes in extratropical planetary wave 
driving (Figure 5). The detailed effects of the QBO on tracer 
transport are discussed in Section 5. 

3. DYNAMICS OF THE QBO 
3. 1 . QBO mechanism 

Since the QBO is approximately longitudinally symmet- 
ric, it is natural to try to explain it within a model that con- 
siders the dynamics of a longitudinally symmetric atmos- 
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phere. In a rotating atmosphere the temperature and wind 
fields are closely coupled and correspondingly both heating 
or mechanical forcing (i.e. forcing in the momentum equa- 
tions) can give rise to a velocity response. Although, as 
noted in Section 1, the current view is that mechanical forc- 
ing, provided by wave momentum fluxes, is essential for the 
QBO, the coupling between temperature and wind fields 
must be taken into account to explain many aspects of the 
structure. 

The essence of the mechanism for the oscillation may be 
most simply demonstrated in a simple representation of the 
interaction of vertically propagating gravity waves with a 
background flow that is itself a function of the vertical 
[Plumb* 1977]. Consider two discrete upward-propagating 
internal gravity waves, forced at a lower boundary with 
identical amplitudes but equal but opposite zonal phase 
speeds. The waves are assumed to be quasi-linear (interact- 
ing with the mean flow, but not with each other), steady, 
hydrostatic, unaffected by rotation, and subject to linear 
damping. The superposition of these waves corresponds 
exactly to a single “standing” wave. As each wave compo- 
nent propagates vertically, its amplitude is diminished by 
damping, generating a force on the mean flow due to con- 
vergence of the vertical flux of zonal momentum. This force 
locally accelerates the mean flow in the direction of the 
dominant wave's zonal phase propagation. The momentum 
flux convergence depends on the rate of upward propaga- 
tion, and hence on the vertical structure in zonal wind speed. 
Even with two equal but opposite waves in a basic state at 
rest, the mean flow is an unstable equilibrium; any small 
deviation from zero will inevitably grow with time. 

Plumb [1977] showed that the zonal wind anomalies de- 
scend in time, as illustrated in Figure 10. Each wave propa- 
gates vertically until its group velocity is slowed and the 
wave is damped as it encounters a shear zone where |I7-c| is 
small. As the shear zone descends (Figure 10a) the layer of 
eastward winds becomes sufficiently narrow that viscous 
diffusion destroys the low-level eastward winds. This leaves 
the eastward wave free to propagate to high levels through 
westward mean flow (figure 10b), where dissipation and the 
resulting eastward acceleration gradually build a new east- 
ward regime that propagates downward (Figures 10c, d). 

The process just described repeats, but with westerly 
shear descending above easterly shear, leading to the forma- 
tion of low-level easterly jet. When the westward jet decays, 
the westward wave escapes to upper levels, and a new east- 
erly shear zone forms aloft. The entire sequence, as de- 
scribed, represents one cycle of a nonlinear oscillation. The 
period of the oscillation is determined, among other things, 
by the eastward and westward momentum flux input at the 
lower boundary and by the amount of atmospheric mass 
affected by the waves. In Plumbs [1977] Boussinesq formu- 
lation the QBO period is inversely proportional to momen- 
tum flux. A quasi-compressible atmosphere results in a sub- 


stantially shorter period due to the decrease of atmospheric 
density with height. If the waves are “underdamped” (able 
to generate an oscillation at all heights) the QBO period is 
once again inversely proportional to momentum flux but is 
also proportional to a density scale height. 

Simple representations such as Plumb's capture the essen- 
tial wave mean-flow interaction mechanism leading the 
QBO. However they cannot explain why the QBO is an 
equatorial phenomenon (notwithstanding its important links 
to the extratropics). Some simple insights on this point come 
from considering the equations for the evolution of a 
longitudinally symmetric atmosphere subject to mechanical 
and thermal forcing. A suitable set of model equations for 
such a longitudinally symmetric atmosphere is as follows: 

du 

— — 2£2sin<|)v = F (2) 


2^2 sin 


du R dT 
dz ciH 3<j) 


= 0 


dT_ 

a7” 



( 3 ) 

( 4 ) 


\ d A i a 

— VCOS <J)H 

acos<|> d ( |> p () dz 


(p„vv)=0 


( 5 ) 


Here <(> is longitude, £2 is the angular velocity of the 


Earth, a is the radius of the Earth and p 0 is a nominal basic 
state density proportional to exp (- z !H), In (4), N 1 is the 
square of the buoyancy frequency (a measure of static 
stability), defined as 
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where T 0 is a reference temperature profile depending only 
on z and k = R/c p , where c p is the specific heat of air at 
constant pressure. Finally u is the longitudinal component of 
wind, T is the temperature deviation from 7o, and v and w 
are respectively the latitudinal and vertical components of 
velocity. 

Equation (2) states that the longitudinal acceleration is 
equal to the applied force, F (which will for the moment be 
taken to be a given function of latitude, height and time), 
plus the Coriolis force associated with the latitudinal veloc- 
ity. The Coriolis force is the important effect of rotation — 
the response to an applied force is not simply an equivalent 
acceleration. Instead, part of the applied force is balanced by 
a Coriolis force — how much depends on how large a latitu- 
dinal velocity is excited. Equation (3) is the thermal wind 
equation coupling the longitudinal velocity field and the 
temperature field, which follows from the assumption that 
the flow is in hydrostatic and geostrophic balance. Equation 
(4) states that the rate of change of temperature is equal to 
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the diabatic heating plus the adiabatic heating associated 
with vertical motion. Here the diabatic heating is repre- 
sented by the term - a T, where a is a constant rate, repre- 
senting long-wave heating or cooling. Equation (5) is the 
mass-continuity equation. 

Equations (2-5) may be regarded as predictive equations 
for the unknowns du /dt ,dT /dt , v and w. They may be 
combined to give a single equation for one of the unknowns, 
with a single forcing term containing the force F. In order to 
simplify the calculations and focus on the phenomena of 
importance, it is convenient to follow Garcia [1987] and 
assume that the time dependence is purely harmonic. Thus 
we write F(<\>, z, t) = Re ) and consider the re- 

sponse in the longitudinal velocity u, assumed to be of the 
form ufaz,t)= Re(u(<|),z>'“' )■ (2-5) can then be trans- 
formed to a single equation 
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The operator acting on u on the left-hand side of the 
equation is elliptic, consistent with the well-known property 
of rotating, stratified systems that localized forcing gives 
rise to a non-local response. For an oscillation with period 
two years, co - 10'Y 1 . The Newtonian cooling rate, a, is, for 
the lower stratosphere, generally taken to be about 5X10 , 
corresponding to a time scale of about 20 days. Hence the 
factor 1 + oc/(/co) appearing in the second term on the left- 

hand side may be approximated by a /(/co). 

A scale analysis of (6) shows that when rotational effects 
are weak, i.e. when sin <j) is small, the dominant balance will 
be between the forcing term and the first term on the left- 
hand side. This implies that the acceleration is equal to the 
applied force. More generally the second term on the left- 
hand side will play a major role in the balance, implying 
that the Coriolis force must be substantially canceling the 
applied force in (2). Following Haynes [1998], a quantita- 
tive comparison of the two terms on the left-hand side of (6) 
may be performed by assuming a height scale D and a lati- 
tudinal scale L for the velocity response. Then, at low lati- 
tudes, noting that sin scales as <j) and hence as Ua , the 
ratio of the second term on left-hand side of (6) to the first is 
4ft 2 L 4 a/(a 2 D 2 A 2 (o). It follows that if L«{aDN/(2£l)) m 
(co/a) 1/4 then the acceleration is approximately equal to the 
applied force. This might be called the “tropical re- 
sponse” — it occurs if the latitudinal scale L is small enough. 


On the other hand if Ly>(aDNl( 2Q)) 1/2 (to/a) 1/4 then the ap- 
plied force is largely cancelled by the Coriolis torque and 
most of the response to the applied force appears as a circu- 
lation in the height-latitude plane. This might be called the 
“extratropical response” (though clearly the scaling would 
require modification if it were to be applied well away from 
the equator). 

The physical reason for the distinction between the tropi- 
cal and extratropical responses is the link between velocity 
and temperature fields in a rotating system, expressed by 
(3), together with the temperature damping implied by (4). 
At high latitudes an applied force, varying on sufficiently 
long time scales, will tend to be cancelled by the Coriolis 
force due to a circulation in the height-latitude plane. This 
circulation will induce temperature anomalies, upon which 
the thermal damping will act, effectively damping the veloc- 
ity response and limiting its amplitude. At low latitudes, on 
the other hand, the force will give rise to an acceleration, 
there will be relatively little temperature response and ther- 
mal damping will have little effect on the velocity response. 
It is as if low-latitude velocities have a longer “memory” 
than high-latitude velocities; anomalies at low latitudes take 
longer to dissipate. Thus the QBO mechanism might be ex- 
pected to work only at low latitudes. The Lindzen and Hol- 
ton (1968) experiments in a 2-D model showed that the 
Coriolis torque reduced the amplitude of the wind oscilla- 
tion away from the equator. Haynes [1998] went further to 
suggest that it might be the transition from the tropical re- 
gime to the extratropical regime that sets the latitudinal 
width of the QBO, rather than, for example, the latitudinal 
scale of the waves that provide the necessary momentum 
flux. Simulations in a simple numerical model where the 
momentum forcing is provided by a latitudinally broad field 
of small-scale gravity waves, designed not to impose any 
latitudinal scale, predicted a transition scale at about 10 de- 
grees. 

To summarize, a long-period oscillation that requires the 
zonal velocity field to respond directly to a wave-induced 
forced is likely to work only in the tropics, since elsewhere 
the force will tend to be balanced by the Coriolis torque due 
to a meridional circulation. For this reason, 1-D models, 
which omit Coriolis torques altogether, can capture the 
tropical oscillation. However, they cannot simulate the lati- 
tudinal structure that arises in part from the increase of 
Coriolis torques with latitude. 

3.2. Waves in the tropical lower stratosphere 

There exists a broad spectrum of waves in the tropics, 
and the details of which waves are driving the QBO are not 
entirely settled. Based on observations of wave amplitudes, 
we now believe that a combination of Kelvin, Rossby- 
gravity, inertia-gravity and smaller-scale gravity waves pro- 
vide most of the momentum flux needed to drive the QBO 
[ Dunkerton , 1997]. All these waves originate in the tropical 
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troposphere and propagate vertically to interact with the 
QBO. Convection plays a significant role in the generation 
of tropical waves; modes are formed through lateral propa- 
gation, refraction, and reflection within an equatorial 
waveguide, the horizontal extent of which depends on wave 
properties (e.g., turning points where wave intrinsic fre- 
quency equals the local inertial frequency). 

Equatorward propagating waves originating outside the 
tropics, such as planetary Rossby waves from the winter 
hemisphere, may have some influence in upper levels of the 
QBO [OrtlancL 1997]. The lower region of the QBO (~20- 
23 km) near the equator is relatively well shielded from the 
intrusion of extratropical planetary waves [O'Sullivan, 
1997]. 

Vertically propagating waves relevant to the QBO are ei- 
ther those with slow vertical group propagation undergoing 
absorption (due to radiative or mechanical damping) at such 
a rate that their momentum is deposited at QBO altitudes, or 
those with fast vertical group propagation up to a critical 
level lying within the range of QBO windspeeds [Dunker- 
ton , 1997]. The height at which momentum is deposited 
depends on the vertical group velocity (supposing for argu- 
ment’s sake that the damping rate per unit time is independ- 
ent of wave properties). Waves with very slow group propa- 
gation are confined within a few kilometers of the tro- 
popause [Li et al . , 1997]. Waves with fast vertical group 
velocity (and phase speeds lying outside the range of QBO 
windspeeds) propagate more or less transparently through 
the QBO. 

Long-period waves tend to dominate spectra of horizontal 
wind and temperature. Higher frequency waves are expected 
to transport momentum to a greater degree than would be 
suggested by a comparison of temperature variance alone. 
We can organize the waves relevant to the QBO into three 
categories: 

(1) Kelvin and Rossby-gravity waves, which are equato- 
rially trapped; periods > 3days; wavenumbers 1-4 

(zonal wavelengths > 10000 km). 

(2) inertia- gravity waves , which may or may not be 
equatorially trapped; periods of ~l-3 days; 
wavenumber -4-40 (zonal wavelengths -1000- 
10000 km). 

(3) Gravity waves; periods <:1 day; wavenumber >40 

(zonal wavelengths -10-1000 km) propagating rap- 
idly in the vertical. (Waves with very short horizontal 
wavelengths (<10 km) tend to be trapped vertically at 

tropospheric levels near the altitude where they are 
forced, and are not believed to play a significant role 
in middle atmosphere dynamics.) 

The observations reviewed below suggest that intermedi- 
ate and high-frequency waves help to drive the QBO. How- 


ever, uncertainties remain in the wave momentum flux spec- 
trum, with regard to actual values of flux and the relative 
contribution from various parts of the spectrum. Although 
the momentum flux in mesoscale waves is locally very 
large, it is necessary to know the spatial and temporal distri- 
bution of these waves in order to assess their role in the 
QBO. Available observations are insufficient for this pur- 
pose. For intermediate-scale waves, it is unclear what frac- 
tion of the waves are important to the QBO without a more 
precise estimate of their phase speeds, modal structure, and 
absorption characteristics. Twice-daily rawinsondes provide 
an accurate picture of vertical structure but have poor hori- 
zontal and temporal coverage. Their description of horizon- 
tal structure is inadequate, and temporal aliasing may occur, 
obscuring the true frequency of the waves. 

The QBO, in principle, depends on wave driving from the 
entire tropical belt, but the observing network can only 
sample a small fraction of horizontal area and time. Thus it 
is uncertain how to translate the information from local ob- 
servations of intermediate and small-scale waves into a use- 
ful estimate of QBO wave driving on a global scale. Ulti- 
mately satellite observations will provide the needed cover- 
age in space and time. These observations have already 
proven useful for planetary-scale equatorial waves and 
small-scale extratropical gravity waves with deep vertical 
wavelength. Significant improvement in the vertical resolu- 
tion of satellite instruments and their ability to measure or 
infer horizontal wind components will be necessary, how- 
ever, before such observations are quantitatively useful for 
estimates of momentum flux due to intermediate and small- 
scale waves in the QBO region. 

3.2.1. Kelvin and Rossby-gravity waves 

Kelvin and Rossby-gravity waves were detected many 
years ago, using rawinsonde observational data, by Yanai 
and Maruyama [1966] and Wallace and Kousky [1968]; 
these discoveries were important to the development of a 
modified theory of the QBO by Holton and Lindzen [1972]. 
For reviews of early equatorial wave observations see Wal- 
lace, [1973]; Holton , [1975], Cornish and Larsen , [1985], 
Andrews et al . , [1987], and Dunkerton, [1997]. Interpreta- 
tion of disturbances as equatorial wave modes relies on a 
comparison of wave parameters (e.g., the relation of hori- 
zontal scale and frequency), latitudinal structure (e.g., sym- 
metric or antisymmetric about the equator) and phase rela- 
tionship between variables (e.g., wind components and tem- 
perature) with those predicted by theory. The identification 
of equatorial modes is relatively easy in regions with good 
spatial coverage so that coherent propagation may be ob- 
served. 

Long records of rawinsonde data from high-quality sta- 
tions have been used to derive seasonal and QBO-related 
variations of Kelvin and Rossby-gravity wave activity near 
the equator [Maruyama, 1991; Dunkerton, 1991b, 1993; 
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Shiotani and Horinouchi, 1993; Sato et ai , 1994; Wikle et 
ai, 1997]. The QBO variation of Kelvin wave activity ob- 
served in fluctuations of zonal wind and temperature is con- 
sistent with the expected amplification of these waves in 
descending westerly shear zones. Annual variation of 
Rossby-gravity wave activity is observed in the lowermost 
equatorial stratosphere and may help to explain the observed 
seasonal variation of QBO onsets near 50 hPa [Dunkerton, 
1990]. 

Equatorially trapped waves have been observed in tem- 
perature and trace constituent data obtained from various 
satellite instruments. Most of these studies dealt with waves 
in the upper stratosphere relevant to the stratopause semian- 
nual oscillation; a few, however, also observed waves in the 
equatorial lower stratosphere relevant to the QBO (e.g., 
Salby et al [1984]; Randel [1990]; Ziemke and Stanford 
[1994]; Canziani et ai [1995]; Kawamoto et ai [1997]; 
Shiotani et ai [1997]; Mote et ai [1998]; Canziani and Hol- 
ton [1998]). It is difficult to detect the weak, shallow tem- 
perature signals associated with vertically propagating equa- 
torial waves, and satellite sampling usually recovers only 
the lowest zonal wavenumbers (e.g. waves 1-6). Neverthe- 
less, satellite observations are valuable for their global view, 
complementing the irregular sampling of the rawinsonde 
network. 

2-D modeling studies [Gray and Pyle [1989], Dunkerton , 
1991a; 1997] showed that Kelvin and Rossby-gravity waves 
are insufficient to account for the required vertical flux of 
momentum to drive the QBO. The required momentum flux 
is much larger than was previously assumed because the 
tropical stratospheric air moves upward with the Brewer- 
Dobson circulation. When realistic equatorial upwelling is 
simulated in the model, the required total wave flux for a 
realistic QBO is 2-4 times as large as that of the observed 
large-scale, long-period Kelvin and Rossby-gravity waves. 
Recent 3-D simulations (e.g., Takahashi and Boville [1992]; 
Hayashi and Golder [1994]; Takahashi [1996]) described in 
Section 3.3.2, confirm the need for additional wave fluxes. 
Therefore it is necessary to understand better from observa- 
tions the morphology of inertia-gravity and gravity waves 
and their possible role in the QBO. 

3.2.2. Inertia-gravity waves 

Eastward propagating equatorial inertia-gravity waves are 
seen in westerly shear phases of the QBO, while westward 
propagating waves are seen in easterly shear phases. Obser- 
vational campaigns using rawinsondes have provided data 
with high temporal and vertical resolution, so that analysis 
is possible both for temporal and vertical phase variations. 

Cadet and Teitelbaum [1979] conducted a pioneering 
study on inertia-gravity waves in the equatorial region, ana- 
lyzing 3-hourly rawinsonde data at 8.5°N, 23.5°W during 
GATE (Global Atmospheric Research Project Atlantic 
Tropical Experiment). The QBO was in an easterly shear 


phase. They detected short vertical wavelength (<1.5 km) 
inertia- gravity wave-like structure having a period of 30-40 
hours. The zonal phase velocity was estimated to be west- 
ward. 

Tsuda et ai [1994a,b] conducted an observational cam- 
paign focusing on waves in the lower stratosphere at Watu- 
kosek, Indonesia (7.6°S, 112.7°E) for 24 days in February to 
March 1990 when the QBO was in a westerly shear phase. 
Wind and temperature data were obtained with a temporal 
interval of 6 hours and vertical resolution of 150 m. Figure 
1 1 shows a time and height section of temperature fluctua- 
tions with periods shorter than 4 days. Clear downward 
phase propagation is observed in the lower stratosphere 
(above about 16 km altitude). The vertical wavelength is 
about 3 km and wave period is about 2 days. Similar wave 
structure was seen also for zonal («) and meridional wind 
(v) fluctuations. The amplitudes of horizontal wind and 
temperature fluctuations were about 3 ms' 1 and 2K, respec- 
tively. 

Based on hodographic analysis assuming that these fluc- 
tuations are due to plane inertia-gravity waves, Tsuda et ai 
[1994b] showed that most waves propagated energy east- 
ward and upward in the lower stratosphere. Similar charac- 
teristics are observed in their second campaign in Bandung, 
Indonesia (107.6°E, 6.9°S) during another eastward shear 
phase of the QBO (November 1992 to April 1993) [Shimizu 
et ai, 1997]. 

Statistical studies of equatorial inertia-gravity waves have 
been made using operational rawinsonde data at Singapore 
(1.4°N, 104.0°E). Maruyama [1994] and Sato et ai [1994] 
analyzed the year-to-year variation of 1-3 day wave activity 
in the lower stratosphere using data from Singapore span- 
ning ten years. Extraction of waves by their periods has 
merit since the ground-based wave frequency is invariant 
during the wave propagation in a steady background field. 
The QBO can be considered sufficiently steady for these 
purposes for inertia-gravity waves having periods shorter 
than several days. 

Maruyama analyzed the covariance of zonal wind and the 
time derivative of temperature for 1-3 day components and 
estimated the vertical flux of zonal momentum per unit den- 
sity u'w' using the following relation derived from the 
thermodynamic equation for adiabatic motions: 



where T is temperature, t is time, u and w are the zonal and 
vertical components of wind velocity, c the ground-based 
horizontal phase speed ( c =c-u ) the intrinsic horizontal 
phase speed, u the background wind speed, and the overbar 
indicates a time average. Since C is not obtained from the 
observational data, this estimate is possible only when u is 
small enough to assume c/c ~ 1 . Maruyama showed that the 
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momentum flux u'w ' is largely positive and that the magni- 
tude is comparable to that of Kelvin waves in the westerly 
shear phase of the QBO. 

Sato et al. [1994] examined interannual variation of the 
power and cross spectra of horizontal wind and temperature 
fluctuations in the period range of 1-20 days at Singapore. 
They found that spectral amplitudes are maximized around 
the tropopause for all components in the whole frequency 
band, although the altitudes of the tropopause maxima are 
slightly different. The T and u spectra are maximized around 
10-day period, corresponding to Kelvin waves. In the lower 
stratosphere, the wave period shortens, e.g., 9 days at 20 km 
to 6 days at 30 km. On the other hand, v spectra are maxi- 
mized around 5 days slightly below the tropopause corre- 
sponding to Rossby-gravity waves. The Rossby-gravity 
wave period also becomes shorter with increasing altitude in 
the lower stratosphere, consistent with the analysis of Dun - 
kerton [1993] based on rawinsonde data at several locations 
over the tropical Pacific. An important fact is that spectral 
amplitudes are as large at periods shorter than 2-3 days as 
lor long-period Kelvin waves and Rossby-gravity waves. 

The activity of inertia-gravity and Kelvin waves is ob- 
served to be synchronized with the QBO. Figure 12 shows 
power and cross spectra as a function of time averaged over 
the height region 20-25 km in the lower stratosphere. A 
lowpass filter with cutoff of 6 months was applied in order 
to display the relation with the QBO more clearly. Dominant 
peaks in the power spectra of T and u are observed in the 1- 
3 day period range during both phases of the QBO and 
around 10-day period in the westerly shear phase of the 
QBO. The latter peak corresponds to Kelvin waves. 

The quadrature spectra Qti,(w) correspond to the covari- 
ance ot zonal wind and time derivative of temperature. Thus 
large negative values observed around 10-day period in the 
westerly shear phase show the positive u\v' associated with 
Kelvin waves [Maruyama, 1991, 1994], Such tendency is 
not clear at shorter periods in the quadrature spectra. Clear 
synchronization with the QBO is seen in the cospectra 
C Tu { o>) in the whole range of frequencies. Positive and nega- 
tive values appear alternately in the westerly and easterly 
shear phases, respectively, though the negative values are 
weak around 1 0-day period. This feature cannot be ex- 
plained by the classical theory of equatorial waves in a uni- 
form background wind [ Matsuno , 1966], which predicts that 
the covariance of T and u should be essentially zero. 

Dunkerton [1995] analyzed theoretically and numerically 
the covariance of T and u for 2-D (plane) inertia-gravity 
waves in a background wind having vertical shear, and de- 
rived the following relation 



waves. This theory was extended to 3-D equatorially 
trapped waves [ Dunkerton , 1998]. According to (8), the 
covariance is proportional to the vertical shear and vertical 
flux of horizontal momentum, or radiation stress. The sign 
of the covariance is determined by the vertical shear, inde- 
pendent of the horizontal and vertical direction of inertia- 
gravity wave propagation. This is qualitatively consistent 
with the observation in Figure 12c. 

Sato and Dunkerton [1997] estimated momentum fluxes 
associated with 1-3 day period waves directly and indirectly 
based on the quadrature and cospectra of T and u compo- 
nents at Singapore obtained by Sato et al. , [1994]. Unlike 
Kelvin waves which propagate only eastward, inertia- 
gravity waves can propagate both eastward and westward. 
Thus the net momentum flux estimate from quadrature spec- 
tra as obtained by Maruyama [1994] may be a residual after 
cancellation between positive and negative values. On the 
other hand, cospectra correspond to the sum of absolute 
values of positive and negative momentum fluxes. Using an 
indirect estimate of momentum fluxes from cospectra and a 
direct one from quadrature spectra, positive and negative 
parts of momentum fluxes can be obtained separately. 

The direct estimate for Kelvin waves (5-20 day period) is 
2-9xl0~ 3 m 2 s~ 2 and accords with the indirect estimate 
within the estimation error, supporting the validity of the 
indirect method. The result for 1-3 day periods is shown in 
Figure 13, assuming plane inertia-gravity waves. The indi- 
rect estimate of 1-3 day period components for westerly 
shear is 20-60xl0” 3 m 2 s ~ 2 , while the direct estimate is 
only 0 - 4x lO’Vi 2 ^* 2 . For easterly shear, the indirect esti- 
mate is 10 - 30 x 10 3 mV 2 , while the direct estimate is 
almost zero. The discrepancy between the indirect and direct 
estimates indicates a large cancellation between positive and 
negative momentum fluxes. 

There is ambiguity in the indirect estimate according to 
the assumed wave structure. If equatorially trapped modes 
are assumed, the values should be reduced by 30-70%. On 
the other hand, if there is aliasing from higher frequency 
waves (wave periods shorter than 1 day for twice-daily data 
such as rawinsonde data at Singapore), the actual momen- 
tum flux should be much larger than shown in Figure 13. 
Even considering these ambiguities, it appears that interme- 
diate-frequency period gravity or inertia- gravity waves have 
significant momentum flux compared to Kelvin and 
Rossby-gravity waves. 

Bergman and Salby [1994] calculated equatorial wave ac- 
tivity propagating into the stratosphere based on high- 
resolution imagery of the global convective pattern and 
some simple assumptions about the relation of cloud varia- 
tions to the properties of the waves that would be generated. 
Figure 14 shows the geographical distribution of the vertical 
component of Eliassen-Palm flux that they derived. The 
components with periods shorter than 2 days have large 
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Eliassen-Palm flux, compared with longer period waves. 
The short-period wave generation is large over the African 
and American continents and in a broad area from the Indian 
Ocean to western tropical Pacific. 

The analysis of Bergman and Salby does not provide a 
quantitative estimate of actual wave fluxes, but supports the 
idea that intermediate and small-scale waves contribute sig- 
nificantly to the QBO. From their analysis we note, First, 
that the contribution is large relative to that of planetary- 
scale equatorial waves, by a factor of about 2.5, and second, 
that most of the activity in smaller-scale waves is associated 
with zonal phase speeds lying in the range of QBO wind- 
speeds [ Dunkerton , 1997]. 

The excitation of inertia-gravity waves by deep tropical 
convection occurs either through a process of self- 
organization, in which waves and convection support one 
another, or (more simply) as a result of irregular, seemingly 
random activity as convective elements impinge on a strati- 
fied layer above. The self-organization of waves and con- 
vection occurs mainly on longer horizontal and temporal 
scales as illustrated (e.g., Takayabu et al [1996], Wada et al 
[1999], and Wheeler and Kiladis [1999]). 

3.2.3. Gravity waves 

Deep convection is also a dominant source for high- 
frequency gravity waves near the equator. Numerical simu- 
lations of convection that include stratospheric altitudes 
[Fovell et ai, 1992; Alexander et al, 1995; Alexander and 
Holton , 1997] generate high-frequency gravity waves that 
appear prominently above storm clouds (Figure 15). Theory 
predicts a close association of high-frequency waves in the 
lower stratosphere and the storms that generate them be- 
cause the energy propagation direction is most vertically 
oriented for the highest frequency waves, 

co = N cos 0 (9) 

where 0 is the angle between the vertical and the lines of 
constant phase or direction of the group velocity vector. In 
these simulations, the high-frequency waves carry signifi- 
cant momentum flux suggesting they could play a promi- 
nent role in driving the QBO [Alexander and Holton, 1997]. 

Observational studies have identified large-amplitude, 
high-frequency gravity waves in the stratosphere directly 
above deep tropospheric convection [Rottger, 1980; Larsen 
et ai , 1982; Sato, 1992, 1993; Sato et al, 1995]. Sato 
[1993] was able to estimate the vertical momentum flux 
carried by the waves. These midlatitude observations found 
significant values of 0.03 Pa, an order of magnitude larger 
than the estimated time-average, zonal mean flux needed for 
tropical gravity waves to be important to the QBO [Dunker- 
ton, 1997]. Note that momentum flux is properly measured 
in units of Pascal (Pa), equal to air density times the product 
of velocity components. In most QBO literature the density 
term is ignored, and the resulting “flux is described in units 
of mV 2 . Near the tropical tropopause the density is about 


0.1 in MKS units, providing an easy conversion between the 
two definitions of flux. 

High-frequency gravity waves have also been detected in 
aircraft observations in the lower stratosphere. NASA's ER- 
2 aircraft flies at altitudes up to ~20 km and has participated 
in numerous campaigns that included tropical flights in the 
stratosphere. Observations of winds and temperatures on- 
board have been used to detect gravity waves. Pfister et al. 
[1986, 1993a,b] detected short horizontal wavelength 
waves < 150 km in temperature and horizontal wind fields 
associated with cumulus convection over Panama and 
Northern Australia. Pfister et al [1993a,b] proposed a 
“convective topography” mechanism for generating these 
waves and used a model to estimate the vertical momentum 
flux that could be generated by this mechanism and the im- 
pact such waves might have on the QBO momentum 
budget. The calculated effect was small (< 10%) compared 
to estimates of planetary-scale wave driving. Other forcing 
mechanisms may be active, however, and the effects of 
high-frequency waves would be under-represented in these 
calculations [ Dunkerton , 1997]. The estimates are also un- 
certain because of unknown geographical distribution and 
occurrence frequency of convectively-forced gravity waves 
based on only a few case studies. 

Alexander and Pfister [1995] used observations of both 
the horizontal and vertical winds together to estimate mo- 
mentum flux along an ER-2 flight path over deep convec- 
tion north of Australia. The momentum flux emphasizes 
shorter-period waves. Very large values -0.1 Pa were ob- 
served over the deepest highest clouds. More extensive cor- 
relations of these data with cloud top temperature were re- 
ported by Alexander et al [1999] and are shown in Figure 
16. These results suggested that large values of momentum 
flux are correlated with deep convection and have similar 
magnitudes to waves generated in 2-D simulations ol tropi- 
cal convection [Alexander and Holton , 1997]. While the 
geographical and seasonal distributions of such fluxes are 
uncertain, the results nonetheless suggest a potentially im- 
portant role for high-frequency gravity waves in the QBO. 

The gravity waves observed in the low-latitude (~12°S) 
rawinsonde analyses described by Allen and Vincent [1995] 
showed a seasonal cycle suggesting convection as an impor- 
tant source during the December-February monsoon season. 
These results represented only one year of observations; 
however, subsequent analysis of data from Cocos Island, 
also at ~12°S spanning 6 years, showed a similar correlation 
with the monsoon season, but it is modulated by QBO 
winds [Vincent and Alexander, 1999]. The months with peak 
momentum flux were found to coincide with the strongest 
easterly winds. At these times the waves also propagate pre- 
dominantly eastward. Parallel theoretical calculations of 
gravity wave propagation and interaction with the back- 
ground flow qualitatively support the convective topography 
mechanism because of the anisotropy observed in wave 
propagation direction, and because the altitude where the 
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waves are generated was inferred to be very high, close to 
the tropopause. Karoly et ai [1996] also observed a correla- 
tion between inertia-gravity wave activity and deep convec- 
tion in tropical sounding data. 

3.3. Numerical models of the QBO 

3.3. 1 . One-dimensional models 

1-D models, in which wind and wave fields are functions 
of height only have been used to explore many aspects of 
QBO behavior relevant in the real atmosphere. Sometimes 
the one-dimensionality may be justified, as in Holton and 
Lindzen [1972], by integrating the dynamical equations in 
latitude to derive evolution equations for the zonal flow 
integrated across the tropics. This approach may be accept- 
able for Kelvin waves, which are known to depend primar- 
ily on the near-equatorial wind and which generate a simple 
profile of the mean-flow forcing. Other equatorial waves, 
however, are significantly affected by latitudinal shear and 
produce more complicated latitudinal profiles of mean-flow 
forcing [Andrews and McIntyre , 1976; Boyd, 1978; Dunker- 
ton , 1983a]. High-frequency gravity waves propagate 

mostly vertically, and their interaction with the mean flow is 
more easily described in a ID model. For these waves the 
latitudinal dependence of sources is important, however, 
owing to the seasonally varying distribution of tropical con- 
vection [Allen and Vincent , 1995], For these and other rea- 
sons the utility of 1-D models is limited and their value lies 
in their simplicity rather than their realism. Plumb [1977] 
(discussed in Section 3.1) included non-rotating gravity 
waves in order to demonstrate a number of basic properties 
of QBO behavior. 

1-D models have been used to investigate several aspects 
of QBO wave forcing, including the effect of scale- 
dependent radiative damping of the waves Hamilton [1981], 
effects of laterally propagating Rossby waves Dunkerton 
[1983b], the effect of self-acceleration of wave phase speed 
and of wave saturation [Tanaka and Yoshizawa , 1987], up- 
welling of the tropical Hadley circulation [Saravanan, 
1990]. Interannual variations of forcing, as might be due to 
ENSO may cause long-term variations of QBO period by 
Geller et al. [1997], 1-D models are also useful for 
interpreting the results of 2-D models [e.g., Dunkerton 

m Two-dimensional Models 

Many important aspects of the QBO behavior can be con- 
sidered only using models that represent latitude as well as 
height variation. This applies self-evidently to questions 
concerning latitudinal structure, but also to wider questions 
concerning the interaction of the QBO with the annual cycle 
and the effect of the QBO on tracer distributions (to be dis- 
cussed in more detail in Section 5). 

The first detailed 2-D model analysis of the latitudinal 
structure, including the circulation in the latitude-height 
plane (the meridional circulation), was that by Plumb and 


Bell [1982b], who assumed that the wave momentum fluxes 
were due to equatorial Kelvin waves and mixed Rossby- 
gravity waves. For the wind field at any instant they calcu- 
lated the height-latitude structure of the waves on the basis 
of a linear, steady-state calculation [Plumb and Bell , 1982a]. 
These momentum and heat fluxes were used to force the 
longitudinally symmetric dynamical equations including 
thermal damping through a Newtonian cooling term, with 
the structure of the waves being recalculated at each time 
step. Their attempt to simulate the QBO in this manner was 
broadly successful, except that they were limited to cases 
where the amplitude of the oscillation was only about half 
of that observed. Dunkerton [1985] and Takahashi [1987] 
used different strategies for calculating the wave momentum 
fluxes and succeeded in simulating oscillations of realistic 
amplitude. The 2-D models showed explicitly that the QBO 
meridional circulation was in the sense of sinking at the 
equator in westerly shear zones and rising at the equator in 
easterly shear zones. In westerly shear zones the thermal 
wind equation (2) implies a maximum in temperature at the 
equator which is maintained against thermal damping by 
adiabatic heating due to sinking motion. The opposite holds 
in easterly shear zones. The pattern of shear zones and 
meridional circulations is shown schematically in Figure 17. 

The models mentioned above were focused on the equa- 
torial regions and did not include a realistic seasonal cycle 
in winds or temperatures. A realistic QBO was achieved by 
Gray and Pyle [1989] in a full radiative-dynamical model 
(which therefore had a realistic seasonal cycle) only by in- 
creasing their parameterized wave momentum forcing by a 
factor of three larger than could be justified by Kelvin and 
Rossby-gravity waves alone. This additional forcing was 
required in order that the QBO wind regimes would propa- 
gate downward despite climatological upwelling in the trop- 
ics [ Dunkerton , 1997]. This was additional evidence of the 
requirement for forcing from waves other than Kelvin and 
Rossby-gravity waves. Mengel et ai [1995] obtained a 
QBO-Iike oscillation in a 2-D middle atmosphere model in 
which eddy momentum transport was due solely to the 
Hines [1997] gravity wave parameterization. Although the 
QBO simulated in this model was weak and sensitive to 
vertical diffusion, the parameterization’s striking ability to 
reproduce the entire phase structure of observational equato- 
rial oscillations [Burrage et ai , 1996] emphasizes the im- 
portance of gravity waves that have a wide range of phase 
speeds and whose amplitude grows with height. 

The 2-D simulations by Gray and Pyle [1989] demon- 
strated that the influence of the QBO extends to all latitudes. 
For example, the QBO circulation with rising or sinking at 
the equator is compensated by an opposing circulation off 
the equator, which gives rise to a temperature anomaly in 
the subtropics and midlatitudes of opposite sign to that at 
the equator (see also Plumb and Bell [1982b] and Dunker- 
ton [1985].) A significant interhemispheric asymmetry in the 
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timing and amplitude of the subtropical anomalies was also 
present, due to the interaction of the QBO with the seasonal 
cycle [Gray and Dunkerton , 1990], The QBO-induced me- 
ridional velocity in the winter hemisphere and, correspond- 
ingly, the vertical velocity in the subtropics winter hemi- 
sphere are substantially larger than those in the summer 
hemisphere, particularly above -25 km [Jones et al 1998; 
Kinnersley , 1998], probably due in part to the asymmetric 
subtropical angular momentum gradients under solstice 
conditions. 

The meridional circulation affects chemical tracers such 
as ozone and gives rise to strong QBO signals in such trac- 
ers at all latitudes, with significant interhemispheric asym- 
metry (see Section 5). However, there may also be a signifi- 
cant feedback of the ozone QBO on the dynamics of the 
QBO, since changes in ozone have radiative implications 
and, in particular, have a direct effect on short-wave heating. 
The effect of including the coupling between QBO ozone 
anomalies and heating rates tends to reduce the heating rate 
that would otherwise be calculated from a given temperature 
anomaly in the lower stratosphere [Hasebe, 1994, Li et al., 
1995, Kinnersley and Pawson, 1996, and Randel et al, 
1999]. Hasebe [1994] argued that this effect must be taken 
into account in order to explain the observed phase relation 
in the lower stratosphere between the QBO signals in ozone 
and wind. In the upper stratosphere the ozone heating en- 
hances the QBO vertical velocity. 

Plumb and Bell [1982b] noted that the advective effects 
of the meridional circulation can account for the observed 
asymmetry in the descent of easterlies and westerlies (with- 
out any need for asymmetry in the waves providing easterly 
and westerly momentum fluxes). Downward advection of 
momentum associated with the westerly shear zone en- 
hances the descent of the westerlies, while the upward ad- 
vection of momentum associated with the easterly winds 
inhibits the descent of the westerlies. An additional advec- 
tive effect, since the vertical velocities are maximal near the 
equator, is to narrow the region of strongest westerly accel- 
eration and to broaden that of strongest easterly acceleration 
[ Hamilton , 1984, Dunkerton and Delisi , 1985, Dunkerton , 
1985, Takahashi , 1987; Dunkerton, 1991]. 

3.3.3. Three-dimensional models (including 
mechanistic and GCMs) 

In the 1-D and 2-D models discussed above, the waves 
that contribute to the driving of the QBO must be param- 
eterized. 3-D models offer the possibility of explicit simula- 
tion of the waves, without any need tor simplifying assump- 
tions that allow parameterization. In “mechanistic” models 
waves are artificially forced (e.g. by imposed heating fields 
or imposed lower boundary perturbations). One particular 
issue where mechanistic models have given important in- 
sight concerns whether the easterly momentum forcing 
needed to account for the equatorial QBO can be supplied 


by Rossby-gravity waves. A 3-D mechanistic model simula- 
tion by Takahashi and Boville [1992] in which a Kelvin 
wave and a Rossby-gravity wave were forced at the lower 
boundary, gave a good representation of the QBO in the 
lower stratosphere. However, the amplitudes of the Kelvin 
wave and, particularly, the Rossby-gravity wave were con- 
siderably stronger than the observed values. This added to 
the evidence that a much broader spectrum of waves was 
needed. 

In general circulation models (GCMs) the waves are gen- 
erated spontaneously in the model (although that is not to 
say that the generation processes are necessarily realistic). 
For various reasons the ability to simulate a realistic QBO is 
seen as a stringent requirement on a GCM. Firstly, the large 
apparent contribution of the QBO to interannual variability 
in the entire middle atmosphere, not just low latitudes (see 
Section 4) means that, ultimately, if a GCM that includes the 
middle atmosphere is to be regarded as realistic then it must 
represent the QBO. Secondly, since the QBO is believed to 
be partially driven by waves excited by cumulus convection 
(on a whole range of scales), the ability or inability of a 
model to simulate the QBO might also have implications for 
the realism of the simulated tropical tropospheric circulation 

The QBO was not simulated in GCMs even if the Kelvin 
and Rossby-gravity waves had realistic amplitudes (CCM2: 
Boville and Randel [1992]; SKYHI: Hayashi and Colder 
[1994]). Until relatively recently no GCM had successfully 
simulated the QBO, or an analogous long-period wave- 
diven oscillation (the term “QBO-like” is used). It appears 
that simulation of the subtle QBO mechanism puts demand- 
ing requirements on a GCM. Here we present an overview 
of the successful simulations and discuss the main features 
upon which a successful simulation depends. 

The first realistic simulation of the QBO in a GCM was 
made by Takahashi [1996]. He used the CCSR/NIES GCM 
with horizontal resolution of T21 (triangular truncation at 
total wavenumber 21, equivalent to a grid spacing of about 
600 km, or 5° of latitude) and with a vertical grid spacing of 
500 m in the stratosphere. This was an unprecedentedly fine 
vertical resolution which allowed waves with small vertical 
wavelengths to propagate vertically and interact with the 
mean flow; most previous GCMs had a vertical grid spacing 
of 2 km or more in the stratosphere. 

All GCMs, for numerical reasons, include some kind ot 
horizontal diffusion or its equivalent. In this case, to obtain 
a QBO it was necessary also to reduce the coefficient of the 
fourth order horizontal diffusion by one order of magnitude 
from its standard value. The model then produced a QBO- 
like oscillation with a period of 1.5 years. Other simulations 
have been made by Horinouchi and Yoden [1998] (aqua- 
planet GCM, period 1.1 years), Hamilton et al. [1999] 
(SKYHI, period 1 year), and Untch [1998] (ECMWF, realis- 
tic period). Figure 18 shows the most realistic simulation to 
date [ Takahashi , 1999], with a period of 2.3 years, using a 
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horizontal resolution to T42. The horizontal diffusion coef- 
ficient had to be reduced by a factor of four from its stan- 
dard value. 

Although several GCMs have produced simulations of 
the QBO, there is no simple set of criteria that will guaran- 
tee a successful simulation. Further, attempts to produce 
more realistic simulations are difficult (as well as time con- 
suming) because the simulations depend on the subtle inter- 
action of several factors. The simulation of a QBO in a 
GCM requires fine vertical resolution in the stratosphere, a 
small diffusion coefficient, moderate to high horizontal 
resolution, and a convection scheme that generates suffi- 
cient waves to drive the QBO. 

The vertical grid spacing in the stratosphere must be suf- 
ficient to resolve the waves and their interaction with the 
mean flow. The values used in simulations of the QBO 
range from 700 m [Horinouchi et al , 1996] to 1500 m 
[Untch, 1998]. Fine horizontal resolution is usually not re- 
quired. The Horinouchi et al. [1996] simulation used a hori- 
zontal resolution of only T21. In a T63 simulation, Untch 
[1998] found that a QBO developed, but then disappeared 
due to a long-term westerly drift in the upper stratosphere. 
The drift was eliminated and the QBO persisted in T159 
simulation. Hamilton et al [1999] found that a horizontal 
resolution of 2° x 2.4° was needed. 

The development of an oscillation of the zonal averaged 
flow can be prevented by the diffusion that smoothes the 
meridional structure. All the QBO simulations have diffu- 
sion time scales longer than the period of the simulated 
QBO-like oscillation. Waves propagating upward into the 
stratosphere are also damped by the diffusion. Takahashi 
[1996] showed that a decrease in the diffusion coefficient 
increased the power of waves in the stratosphere while hav- 
ing little effect on the troposphere. 

As in the real atmosphere, it appears that a broad spec- 
trum of waves supplies the necessary forcing in these simu- 
lations. However, For example, Takahashi et al [1997] sug- 
gested that the easterly acceleration of the model's QBO was 
due to gravity waves as well as Rossby waves from the NH 
winter and Rossby-gravity waves, while westerly accelera- 
tion was due to Kelvin and gravity waves. In contrast, in the 
Takahashi [1999] model gravity waves were the dominant 
forcing for the QBO. 

Horinouchi and Yoden [1998] also did a thorough wave 
analysis. The frequency distribution of the momentum flux 
and its magnitude were roughly consistent with the observa- 
tional estimates of Sato and Dunkerton [1997] for Singapore 
(see Section 3.2). In particular, the Kelvin and Rossby- 
gravity waves played a small role. 

The tropospheric source of the waves, which is domi- 
nated by latent heat release due to cumulus convection, is a 
crucial characteristic. Since most GCMs roughly reproduce 
the climatological mean precipitation, the low-frequency 
components of the latent heating will be similar. However, 


transient characteristics of cumulus convection, which are 
important for the wave excitation, differ greatly among 
models. 

All of the QBO simulations employed the moist convec- 
tive adjustment scheme except for the model of Untch 
[1998], which used the Tiedtke [1989] scheme. The moist 
convective adjustment scheme tends to produce intermittent 
grid-scale cumulus convection. The Takahashi et al [1996] 
T21 model also produced a QBO-like oscillation when the 
scheme was replaced by the prognostic Arakawa-Schubert 
scheme [Nagashima et al, 1998], which also tends to pro- 
duce highly transient grid-scale heating pulses. However, 
some cumulus parameterizations produce little transient 
cumulus convection. For example, convection due to the 
Zhang and McFarlane [1995] scheme in the NCAR CCM3 
(which has not successfully simulated the QBO) results in a 
momentum flux to the stratosphere that is quite weak, even 
though the time-averaged precipitation is very realistic. 

A final factor that may affect the ability of a GCM to 
simulate a QBO is the Brewer-Dobson circulation tropical 
upwelling. As noted in Section 3.3.1, strong upwelling tends 
to slow the descent of the QBO. The upwelling in the Taka- 
hashi [1999] model, which had a realistic period, was a little 
weaker than the observational estimates of Mote et al 
[1996]. The Takahashi [1997] and Horinouchi and Yoden 
[1998] models, which had shorter periods, had unrealisti- 
cally weak upwelling. It is possible that in some models 
unrealistically strong upwelling may prevent the simulation 
of a QBO. 

The most important factor in reproducing the QBO is 
probably the use of a fine vertical resolution to resolve 
equatorial gravity waves. The horizontal diffusion should be 
weak enough not to prevent the evolution of the mean flow 
oscillation and not to be the primary damping mechanism 
for the waves. Transient characteristics of tropical cumulus 
convection are also important since they determine the exci- 
tation of the waves. Despite the recent success, those mod- 
els that reproduced the QBO might have overly active cu- 
mulus convection and hence large amplitude of gravity 
waves with resolved scales. Thus, a parameterization of 
subgrid-scale gravity waves due to convection might be 
necessary to produce the QBO with the same waves that 
drive it in the real atmosphere. 

4. DYNAMICAL EFFECTS IN THE EXTRATROPICAL 
STRATOSPHERE 

Any connection between the equatorial QBO and the ex- 
tratropical atmosphere must be viewed in the context of the 
seasonal cycle and variability of the extratropical strato- 
sphere. The zonal-mean circulation in the stratosphere un- 
dergoes a much stronger seasonal cycle than the tropo- 
sphere, with an actual reversal of winds in the extratropical 
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stratosphere. During the winter season the high-latitude 
stratosphere cools, forming a deep, strong vortex. The 
strong westerlies are replaced by easterlies with increasing 
solar heating in the spring and summer. 

In both hemispheres the smoothly varying seasonal cycle 
described above is modified by the effects of planetary-scale 
waves (herein, simply planetary waves), which are forced 
largely by land-sea contrasts and surface topography. These 
waves propagate vertically and meridionally into the winter 
stratosphere (Figure 5), but are strongly trapped in the mean 
easterly winds of the summer hemisphere [Chcirney and 
Drazin , 1961]. 

The NH has much greater land-sea contrast and larger 
mountain ranges than the SH, resulting in larger amplitude 
tropospheric planetary waves. Consequently, the northern 
winter stratosphere tends to be much more disturbed by 
planetary waves than the southern stratosphere. Large am- 
plitude waves can rapidly disrupt the northern polar vortex, 
even in midwinter, replacing eastward winds with westward 
winds and causing the polar stratosphere to dramatically 
warm. Such events are called major stratospheric warmings. 
The transition from westerlies to easterlies in the springtime 
usually occurs in conjunction with a planetary wave event, 
and is called the final warming. In the NH, the timing of the 
final warming is highly variable, and tends to occur during 
March or April. In the SH the final warming occurs in No- 
vember-December, with much less interannual variability 
[Waugh and Randel , 1999]. 

In the NH, the planetary wave amplitudes are just large 
enough for midwinter sudden warmings to occur during 
some years but not others. Thus, the northern stratosphere is 
sensitive to the effects of vertically propagating planetary 
waves, resulting in large interannual variability in the 
strength of the polar vortex. It appears that this sensitivity to 
the upward and equatorward propagation of planetary waves 
allows the equatorial QBO to influence the polar strato- 
sphere by modulating the flux of wave activity (EP flux) 
[Dunkerton and Baldwin , 1991]. 

The definitive identification of an extratropical QBO sig- 
nal has been difficult due to the shortness and limited spatial 
coverage of data sets. In the NH data up to 10 hPa appear to 
be reliable beginning in the 1950s. Above the 10-hPa level, 
and in the SH lower stratosphere, the lack of rawinsonde 
coverage has limited the production of reliable gridded data 
to the period beginning in the late 1970s, when satellite 
temperature retrievals began. Most of the literature on the 
extratropical influence of the QBO has focused on the NH 
simply because the data record is longer and more reliable. 
Part of the difficulty in identifying a NH QBO signal is that 
the QBO accounts for only a fraction of the variance. In 
addition to the variability of tropospheric forcing, other sig- 
nals, such as the 11 -year solar cycle, volcanoes, and sea 
surface temperature anomalies appear to influence the vari- 
ability of the extratropical stratosphere. 


Holton and Tan [1980, 1982] presented strong empirical 
evidence that the QBO influences the extratropical northern 
stratosphere by using gridded data for 16 NH winters 
(1962-1977) to form easterly and westerly phase compos- 
ites of 50 hPa geopotential. They showed that geopotential 
height at high latitudes is significantly lower during the 
westerly phase of the QBO. They also found a statistically 
significant modulation of the springtime zonal wind in the 
SH. Naito and Hirota [1997] found a robust QBO signal 
during November-December, but a strong relation to the 11- 
year solar cycle during January-February, suggesting that 
solar influence modifies the signal during late winter. 

4. 1 . Mechanism for Extratropical Influence 

The modulation of the polar vortex by the QBO (now 
called the Holton-Tan effect) is relatively easy to under- 
stand, at least in general terms. Typically, the dominant di- 
rection of wave activity propagation for tropospheric plane- 
tary waves is upward and equatorward, and vertical propa- 
gation is limited to the waves with the largest spatial scales 
(primarily waves 1 and 2) [Charney and Drazin, 1961]. In 
the high-latitude stratosphere, these waves distort the vortex 
from zonal symmetry, and if the amplitudes are large 
enough, the vortex can be displaced from the pole or dis- 
rupted so that easterlies replace westerlies near the pole. 
Concomitant with large wave amplitudes is the “ breaking ” 
of planetary waves [McIntyre and Palmer , 1983, 1984] lead- 
ing to the erosion of the vortex and weaker westerly winds. 

The state of the zonal-mean circulation in the NH winter 
stratosphere depends on the degree to which planetary 
waves affect the polar vortex. In the NH wave amplitudes 
are large enough that the time-mean vortex is much smaller 
than it would be in the absence of the waves. By contrast, in 
the SH, wave amplitudes are much smaller (because both 
land-sea contrasts and topography are smaller) and the win- 
tertime climatological zonal winds are approximately twice 
those in the NH. The vertical and meridional propagation of 
planetary waves depends on the latitude-height structure of 
the zonal-mean wind, which may be thought of as refracting 
the waves as they propagate out of the troposphere. Quasi- 
stationary waves cannot propagate in easterly winds, and the 
phase of the QBO (in the tropics and subtropics) alters the 
effective wave guide and position of the boundary between 
easterly and westerly zonal-mean winds (the critical line). 

If the mean flow in the tropical stratosphere is westerly, 
planetary waves are able to penetrate into the tropics and 
even across the equator without encountering a critical line. 
By contrast, when the mean flow in the tropical stratosphere 
is easterly, planetary waves encounter a critical line on the 
winter side of the equator. Thus, when there are easterly 
winds in the tropics, the effective wave guide for planetary 
wave propagation is narrower and the wave activity at mid 
and high latitudes of the winter hemisphere is stronger. 
Stronger planetary waves in high latitudes lead to greater 
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wave-induced drag on the mean flow, reduced westerly 
winds, and hence a warmer winter stratospheric pole. 

4.2. Observations of Extratropical Influence 

The effect of the QBO on the strength of the northern 
winter stratospheric vortex may be seen by comparing com- 
posites of extratropical zonal-mean wind made during east- 
erly and westerly phases of the QBO. The QBO phase must 
be defined precisely, and typically the equatorial mean zonal 
wind at a particular level is used. Holton and Tan [1980] 
defined the QBO phase using equatorial winds at the 50 hPa 
level, but other authors have used 45, 40, and 30 hPa. The 
motivation for picking a particular wind level is to optimize 
the extratropical signal. In the following figures, the QBO 
phase was defined [Wallace et al. , 1993; Baldwin and Dun - 
kerton , 1998a], by using the first two empirical orthogonal 
functions (EOFs) of the vertical variations of equatorial 
winds. The QBO definition was adjusted to optimize the 
extratropical signal in either the NH or SH. The resulting 
definition of the QBO phase is very similar to the 40 hPa 
equatorial wind for the NH composites and the 25 hPa wind 
for the SH composites. 

Figure 19 illustrates the difference in zonal mean wind 
between easterly and westerly QBO composites, using 
NCEP analyses for the period 1978-1996. The difference is 
formed by calculating separate averages of the January wind 
data when the QBO is at 40 hPa easterly and westerly, and 
then taking the difference between these averages. The 
zonal winds were derived from the geopotential fields using 
the balance method [Robinson, 1986; Hitehman et al, 1987; 
Randel, 1 987] . Although this method works well in the ex- 
tratropics, it is not possible to derive accurate winds in the 
tropics, where winds have been interpolated between 10°N 
and 10°S, making the QBO too weak. The northern signal is 
dominated by a modulation of the polar vortex which ex- 
tends from the surface to the 1 hPa level. Differences are of 
the opposite sign south of ~40°N and blend into the upper 
branch of the tropical QBO. 

The features in Figure 19 are not confined to the middle 
stratospheric levels where the equatorial QBO is defined. 
Rather, the most prominent features are found above 10 hPa. 
Correlation analysis of zonal-mean wind in the latitude- 
height plane [Baldwin and Dunkerton , 1991] suggests that 
the upper level features, including those in the SH near 
30°S, represent modulation of the diabatic circulation by the 
QBO. It is reasonable to suppose that this circulation is 
weakly modulated by the Holton-Tan oscillation. For exam- 
ple, during the easterly phase of the QBO, polar tempera- 
tures are relatively high, implying a stronger cross- 
equatorial flow and a weak westerly anomaly near 30°S. 
The feature near 30°S may be interpreted not only as a di- 
rect effect of the QBO, but a remote effect of the Holton- 
Tan wave driving communicated by the mean meridional 
circulation. This behavior is also seen in numerical models 


(see Section 4.3). 

The statistical significance of the NH winter QBO signal 
has been addressed by several authors, and was treated in 
detail by Baldwin and O'Sullivan [1995], The details of 
such an analysis (definition of the QBO, selection of winter 
months, level of data) are critical to the outcome of statisti- 
cal tests. For example, the effect of the QBO is observed to 
be large in December and January, but weaker in February. 
Using 1964-1993 NCEP data for DJF at levels up to 10 hPa, 
they showed that the effect of the QBO (defined by the 40- 
hPa Singapore wind) is statistically significant at the 0.001 
level using 10 hPa geopotential, as measured by the field 
significance test of Barns ton and Livezey [1987]. In zonal- 
mean wind composites the statistical significance of the 
effect of the QBO increases with height, at least through 10 
hPa, with much higher significance at 10 hPa than at 30 
hPa. 

The NH dipole pattern illustrated in Figure 19 is not 
unique to the influence of the QBO - it represents the lead- 
ing mode of variability of the northern winter stratosphere 
[Nigam, 1990; Dunkerton and Baldwin , 1992]. During win- 
ter, the QBO appears to excite the “northern annular mode 
(NAM\ also called the Arctic Oscillation) [Thompson and 
Wallace , 1998, 1999; Baldwin and Dunkerton , 1999], NH 
QBO composites (of geopotential, wind, temperature, etc.) 
tend to reflect the degree to which the QBO excites the 
NAM. One phase of the NAM is represented by a cold, 
strong vortex, while the opposite phase is represented by a 
weaker vortex and higher polar temperatures. The dipole in 
Figure 19 is much more prominent in the stratosphere, and 
weaker in the troposphere. The connection to the strato- 
sphere is limited to the winter season, but the tropospheric 
NAM is observed during all seasons. The tropospheric as- 
pects of the NAM, and the tropospheric patterns associated 
with the QBO will be discussed in Section 6.2. 

Figure 20 illustrates the 5-hPa seasonal development of 
the NH zonal wind composite difference. The extratropical 
signal begins during autumn in midlatitudes, and reaches a 
high-latitude maximum during January. The late winter 
(February and March) composite difference is insignificant 
north of 40°. The abrupt diminution of the signal indicates 
that the QBO modulates the strength of the northern winter 
polar vortex until midwinter, but has little effect on the tim- 
ing of the final warming. 

Major stratospheric warmings are defined as a reversal of 
the zonal-mean zonal wind at 10 hPa, 60°N to easterly, and 
higher temperature at the pole than the zonal average at 10 
hPa, 60°N. The Holton-Tan effect implies that major warm- 
ings should be more common when the QBO is easterly. 
Unfortunately, such a simple measure is not robust because 
the definitions of both QBO phase (EAV) and warming 
(Yes/No) are arbitrary. Using NCEP reanalyzed data from 
1958-1999, and a 40-hPa definition of the QBO phase, 
there were six westerly warmings and ten easterly warm- 
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ings. However, if the "Berlin” hand-analyzed data are used 
with a 45-hPa definition of the QBO, and major warmings 
are defined synoptically, there were ten westerly warmings 
and eleven easterly warmings [K. Labitzke, personal com- 
munication, 2000]. The disparity between these results illus- 
trates that this procedure is too sensitive to the definitions of 
QBO phase and warmings. Composites, as in Figures 19- 
20, are more robust, can be used in both hemispheres, and 
quantify the Holton-Tan effect in both latitude and time. 

The SH polar vortex is much stronger, longer-lived, and 
more quiescent than its NH counterpart. During winter, 
planetary waves generally do not disrupt the southern vortex 
in the lower to middle stratosphere. It is not surprising that 
the observations indicate that the QBO does not signifi- 
cantly modulate the strength of the Antarctic vortex in win- 
ter. However, as shown in Figure 21, the QBO modulates 
the strength of winds in midlatitudes during late autumn, as 
in the NH. Unlike the NH, the wind modulation is apparent 
only in midlatitudes throughout the winter and early spring. 
The striking difference between Figures 21 and 20 is that in 
the SH, the largest influence of the QBO occurs during late 
spring (November), at the time of the final warming. In the 
SH the October vortex is of the same magnitude as the NH 
vortex in January. Since planetary wave amplitudes are 
much smaller in the SH, the effect ot the QBO is seen only 
at the vortex periphery until the vortex is relatively small. 
The 5-hPa observations from Figures 20 and 21 may be 
compared to the (8-hPa) modeling results of O Sullivan and 
Dunkerton [1994], discussed in the next section. 

4.3. Model Simulations of Extratropical Influence 

It is hard to formulate a simple quantitative model for this 
mechanism (analogous, say, to the Holton-Lindzen model 
[Holton and Lindzen, 1972] for interaction of the equatorial 
mean flow with vertically-propagating waves). The main 
complications are (1) the planetary waves propagate both 
vertically and meridionally, (2) the effects of critical lines 
on planetary wave propagation are not easy to predict theo- 
retically. In the absence of a simple theory, the effects of the 
QBO on extratropical waves and mean circulation have 
been studied extensively in detailed numerical simulation 
experiments with models of varying complexity. 

In purely zonally-symmetric models the response of the 
flow to near-equatorial momentum forcing is strongly con- 
fined to tropical and subtropical latitudes [Plumb, 1982]. 
This also appears to be the case for 2-D models in which the 
eddy fluxes are parameterized as a simple linear diffusion 
(e.g., Gray and Pyle [1989]). 

The effect of the QBO on sudden warmings was investi- 
gated in numerical experiments by Dameris and Ebel [1990] 
and Holton and Austin [1991]. Both studies used fairly brief 
integrations of 3-D mechanistic models forced by idealized, 
rapidly growing midlatitude lower boundary perturbations. 
They found that the development of the high-latitude strato- 


spheric flow can be very strongly influenced by the tropical 
winds in the lower stratosphere, although this depends on 
the strength of the imposed wave forcing. Holton and Austin 
found that for weak wave forcing the high latitude flow was 
largely unaffected by the tropical winds, but that this sensi- 
tivity to tropical winds increased as the strength of the wave 
forcing was raised. Within some range of the amplitude of 
the forcing, the model developed a sudden warming when 
there were easterly winds in the tropical stratospheric initial 
conditions, but not when there were westerly winds. As the 
wave forcing was increased further, sudden warmings oc- 
curred regardless of the state of the tropical winds, and the 
sensitivity of the extratropical flow evolution to tropical 
winds diminished. 

Simplified model studies have the advantage of allowing 
the relevant parameters (notably the strength of the wave 
forcing from the troposphere) to be varied in a controlled 
manner. O'Sullivan and Salby [1990] and Chen [1996], used 
models with rather fine horizontal resolution, but limited to 
a single layer in the vertical. Their experiments were run 
with simple linear relaxation of the zonal-mean state in or- 
der to simulate the effects of radiative transfer in constrain- 
ing the mean flow in the real atmosphere, and included an 
imposed wave-1 lower boundary forcing in the winter ex- 
tratropics. The results showed that the models simulated 
high latitude effects of the tropical QBO in the same sense 
as observed. OSullivan and Young [1992] and O’Sullivan 
and Dunkerton [1994] used a global 3-D mechanistic model 
forced with a specified wave-1 perturbation at the 10-km 
lower boundary. The model parameterized radiative effects 
with a linear relaxation of temperature to an imposed radia- 
tive state. Their simulations also included a seasonal cycle. 
Experiments were repeated with a range of amplitudes oi 
the extratropical wave forcing, for initial conditions in the 
tropics representative of easterly and westerly QBO phases. 
The flow in the NH polar region in these simulations was 
found to be nearly unaffected by the tropical wind through 
November. In December, January and February, however, 
there was a significant effect with the time-mean polar vor- 
tex in the stratosphere being stronger when the tropical 
winds were westerly. The contrast in the winter-mean strato- 
spheric polar vortex strength between east and west QBO 
phases depended strikingly on the amplitude of the wave 
forcing adopted and on whether the parameterized gravity 
wave drag was included. 

In a GCM the tropospheric forcing of the stratosphere 
(and its interannual variation) is generated self-consistently 
within the model. Hamilton [1998b] used a GCM run for a 
continuous 48 year period with a time-varying tropical mo- 
mentum forcing that produced a 27-month QBO in the 
equatorial zonal wind with shear zones of realistic magni- 
tude and descending at a realistic speed. Figure 22 shows 
the January composite zonal mean wind for 20 years with 
westerly 40 hPa equatorial winds minus that for 20 years 
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with easterly 40 hPa winds). The results show the tendency 
for a weaker polar vortex in the east QBO phase. The gen- 
eral features (outside of the tropics) of the composite differ- 
ence in Figure 22 compare well with the observed January 
composite shown as Figure 19. Figure 23 shows the easterly 
and westerly phase composite 28-hPa North Pole tempera- 
ture for each month from October through April in the 
model integration. The systematic effects of the tropical 
QBO on NH polar stratospheric temperature in the model 
seem to be largely confined to the period December through 
March. The midwinter differences in Fig. 28 are ~4-5°C, 
which are comparable to those seen in the observed easterly 
minus westerly phase composite of Dunkerton and Baldwin 
[1991]. 

Niwano and Takahashi [1998] have investigated the 
stratospheric extratropical variability in a version of the 
Japanese Center for Climate System Research model, which 
does spontaneously produce a QBO-like oscillation in the 
tropics with period of about 1.4 years. They analyzed a 14- 
year integration and calculated January-March composites 
based on the five most easterly and five most westerly QBO 
phases as judged by the equatorial zonal wind averaged be- 
tween 7 and 50 hPa. The results show that on average the 
NH polar vortex was weaker in the easterly QBO phase, by 
-15 ms' 1 near70°N at 1 hPa. 

In summary, a wide range of models have been applied to 
the study of influence of the QBO on the extratropical 
stratosphere. Thus far all studies seem to have focused ex- 
clusively on the NH in late autumn through early-spring. 
The published model studies are unanimous in showing at 
least some tendency for the strength of the polar vortex 
through the stratosphere to be positively correlated with the 
equatorial zonal wind near 40 hPa. These many theoretical 
model results thus lend credence to the reality of the Holton- 
Tan e I fee t. All the model studies are idealized to one degree 
or another, but the one that is most complete (and hence 
bears detailed comparison with observations) is that of 
Hamilton [1998b]. The results from this study are generally 
in reasonable agreement with available observations, sug- 
gesting that even GCMs that fail to simulate a tropical QBO 
can have a realistic representation of the wave-modulated 
interactions between the low and high-latitude stratosphere. 

4.4. Interaction of the QBO with other low- 
frequency signals 

The extratropical QBO signal may be identified statisti- 
cally in a long data record, but it is only part of the large 
interannual variability of the northern winter stratosphere. 
Several other signals contribute to variability or interact 
with the QBO signal to produce other frequencies of vari- 
ability in the observed data record. 

Since upward-propagating planetary waves require west- 
erly mean winds to penetrate the stratosphere, the ex- 
tratropical influence of the QBO should be limited largely to 


the period between late autumn and early spring. Baldwin 
and Tung [1994] showed that the QBO modulates the ex- 
tratropical annual cycle signal so that the signature of the 
QBO, rather than simply a single spectral frequency peak of 
~28 months, includes two additional spectral peaks at the 
annual frequency plus or minus the QBO frequency. These 
studies demonstrated that the “3-peak QBO” spectrum 
[Tung and Yang , 1994a] can be expected from the Holton- 
Tan effect acting to modulate the annual cycle. They consid- 
ered a harmonic with the period of the QBO that acts to 
modulate a signal consisting of an annual mean plus a sinu- 
soid with an annual period. The combined signal of the 
QBO in the extratropics together with the annual cycle can 
be represented mathematically as 

(A + B sin £2j2f) x sin £2 qbo^ = A sin £2 qbo^ + 

B/ 2 cos (£2j2 - £2QBo)f “ Bl 2 cos (1T2 i2+£2qbo)^ 

where £2| 2 and £2q B0 denote the annual and tropical QBO 
frequencies, respectively. With a QBO period of 30 months 
(the average QBO period during 1979-92, used by Tung and 
Yang), the last two terms of the above equation represent 
variations with periods of 20 and 8.6 months. The 30, 20, 
and 8.6-month spectral peaks were found in ozone [Tung 
and Yang, 1994a, b] ; angular momentum and Eliassen-Palm 
flux [Baldwin and Tung , 1994]; and isentropic potential 
vorticity [Baldwin and Dunkerton , 1998a]. 

Studies using time series of stratospheric temperatures 
[Salby et ai, 1997; Marquardt , 1999] suggest that low fre- 
quency variability of the middle and upper stratosphere in- 
cludes a biennial mode with a period of exactly 24 months. 
Such a purely biennial signal cannot be the result of quasi- 
biennial forcing. Salby et al. as well as Baldwin and Dun- 
kerton [1998a] speculated that a biennial mode might 
propagate into the stratosphere from the upper troposphere. 
It is unclear why a biennial mode, which may be found in 
the troposphere, would be amplified to become important in 
the polar stratosphere. Baldwin and Dunkerton could find 
no explanation for such a biennial mode, and noted that the 
statistical significance of the biennial spectral peak is not 
high — the mode may simply be an artifact of using a short 
(32 year) data record that happened to have biennial vari- 
ability. They noted that the Holton-Tan mechanism would 
tend to make polar anomalies in PV change sign from year 
to year. This tendency, together with random chance, could 
account for the observed biennial mode. If this interpreta- 
tion is correct, then the biennial mode, in all likelihood, will 
not continue. 

Several researchers have considered that remote effects 
from El Nino Southern Oscillation (ENSO) could influence 
the extratropical stratosphere. Such an influence could mas- 
querade as a QBO signal, or at least be difficult to separate 
from a QBO signal. Wallace and Chang [1982] were unable 
to separate the effects of ENSO and the QBO on the tropical 
stratosphere in 21 winters of NH 30-hPa geopotential. Van 
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Loon and Labitzke [1987] also found that the phases of the 
QBO and ENSO tended to coincide. By removing cold and 
warm ENSO years (keeping years only with weak ENSO 
anomalies) they displayed results similar to Holton and Tan. 
Subsequent observational studies (e.g., Hamilton [1993]; 
Baldwin and O'Sullivan, [1995]), and modeling [Hamilton, 
1995] show a consistent picture in which the influence of 
ENSO on the zonal-mean structure of the vortex is largely 
confined to the troposphere. In the lower stratosphere, 
ENSO appears to modulate the amplitudes of large-scale 
stationary waves. 

Decadal variability, possibly related to the 11 -year solar 
cycle, clearly exists in data records which began in the 
1950s. Labitzke [1987] and Labitzke and van Loon [1988] 
studied the observed late-winter NH circulation classified 
by both the level of solar activity and the QBO phase. Naito 
and Hirota [1997] found a robust QBO signal during early 
winter, but they found a strong relation to the solar cycle 
during late winter. Figure 24 summarizes the solar/QBO 
results as scatter plots of mean 30-hPa geopotential heights 
during January and February above the North Pole vs. 10.7 
cm solar radio flux (a proxy for the 11 year cycle in solar 
activity). The data set can be grouped into four categories 
based on the QBO phase and solar activity level. In years 
with low solar activity, the polar winter vortex tends to be 
disturbed and weak when the QBO is easterly, but deeper 
and undisturbed when the QBO is westerly. In years with 
strong solar activity, however, westerly phases of the QBO 
are associated with disturbed winters, whereas easterly 
phases of the QBO are accompanied with deep and undis- 
turbed polar vortices. Hence, the QBO acts as predicted by 
Holton and Tan [1980] in years with low solar activity, but 
appears to reverse its behavior during years with high solar 
activity. Only two cases do not fit into this scheme: 1989 
and 1997. 

It is the subject of active debate whether or not decadal 
variability is caused by the 11 -year solar cycle, but there is 
increasing evidence through modeling that the solar cycle 
has a significant influence on winds and temperatures in the 
upper stratosphere. Over the 11 -year solar cycle, the solar 
“constant” (i.e., the radiative energy input into the earth's 
atmosphere summed over the entire spectrum) varies by less 
than 0.1% [Willson et ai, 1986]. Variability in the UV re- 
sponsible for most of ozone heating is less than 1% 
[ Rottman , 1999]. The variability rises to 8% only at wave- 
lengths shorter than 200 nm, but these wavelengths may 
affect indirectly the ozone chemistry through enhanced pro- 
duction of odd oxygen, which in turn could affect middle 
atmospheric heating rates and dynamics. 

Following earlier solar modeling [Haigh, 1994; 1996; 
1999] and solar-QBO modeling [Rind and Balachandran, 
1995; Balachandran and Rind, 1995], Shindell et ai [1999] 
used a troposphere-stratosphere-mesosphere GCM with 
interactive ozone and realistic values ol UV torcing to show 


that ozone changes amplify irradiance changes to affect cli- 
mate. Circulation changes introduced in the stratosphere 
penetrated downward, even reaching the troposphere. The 
modeling studies find a more intense Hadley circulation 
during solar maximum conditions. They concluded that the 
observed record of geopotential height variations in the NH 
are, in part, driven by solar variability. 

Figure 20 showed that the observed QBO modulation of 
zonal wind in the NH middle stratosphere is essentially over 
by February, and the observations show decadal variability 
coherent with the solar cycle during January-February. The 
possibility exists for the QBO to dominate early winter, 
while solar influence is manifest during late winter. Due to 
the strong absorption of ozone in the UV occurring in the 
upper stratosphere and mesosphere, a solar influence on the 
thermal structure in these regions of the atmosphere is plau- 
sible. This, in turn, might affect the strength of the planetary 
wave driven “extratropical pump” [Holton et ai, 1995]. A 
mechanism involving downward propagation of strato- 
spheric anomalies, through modification of planetary wave 
propagation from below, is discussed in Section 6.2. 

Salby and Callaghan [2000] showed that the QBO west- 
erlies below 30 hPa vary with the solar cycle, as do the east- 
erlies above 30 hPa. Changes in the duration of wind re- 
gimes were found to introduce a systematic drift in the QBO 
phase during northern winter. 

Various hypotheses have been proposed to explain the 
observed stratospheric decadal variability without reference 
to the solar cycle. These hypotheses rely on the QBO inter- 
acting with other signals. Teitelbaum and Bauer [1990] and 
Salby and Shea [1991] argued that the wintertime 11-year 
variability is a byproduct of the analysis procedure which 
involves the stratification of the data into years with respect 
to the QBO. Gray and Dunkerton [1990] showed the possi- 
bility of an 11 -year cycle arising from the beating of the 
QBO with the annual cycle. Salby et ai [1997] and Baldwin 
and Dunkerton [1998a] suggested a modulation of the tropi- 
cal QBO by a biennial extratropical signal (which exists, but 
has not yet been explained) would result in a period of 1 1 
years. These would provide an explanation of the 11 -year 
variability without referring to the solar variability, although 
the observed in-phase relation to the solar cycle remains 
unexplained. 

The origin of the observed biennial or decadal variability 
does not affect the reality of the influence of the QBO on 
the extratropical stratosphere. The debate centers on the 
reality of both the solar and biennial signals, and whether or 
not the observed decadal variability could arise in the ab- 
sence of solar influence. An alternative explanation was 
given by Marquardt [1999]. By analyzing observed data 
rather than analytical functions, he showed that the biennial 
signal can be understood as a natural consequence of the 
modulation of Tung and Yang's extratropical QBO by an 11- 
year signal. 
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5. EFFECTS OF THE QBO ON CHEMICAL 
CONSTITUENTS 

5.1. Background 

There is a substantial body of evidence for the influence 
of the QBO on chemical constituents in the atmosphere. 
Initial evidence came from ground-based observations of 
column ozone at two subtropical stations reported by Funk 
and Garnham [1962], which were shown by Ramanathan 
[1963] to be associated with the stratospheric wind oscilla- 
tion. In an examination of historical ozone data, Angell and 
Korshover [1964] showed that a QBO signal is evident in 
Shanghai data (31°N) in the 1930s. Subsequent information 
about the temporal, latitudinal and vertical structure of the 
ozone QBO has come primarily from satellite observations, 
due to their global nature and better temporal and spatial 
sampling. Although no single satellite has been operational 
for the whole period, the large-scale pattern and evolution of 
the ozone QBO, including its height structure, has been well 
characterized over this period, often with overlapping meas- 
urements from two or more instruments. 

The first simulation of the ozone QBO was carried out by 
Reed [1964] using a simplified linearised model. However, 
it was not until 1986 that the QBO was studied in a full pho- 
tochemical model. Ling and London [1986] included the 
QBO variation in zonal wind in a 1-D radiative dynamical 
photochemical model of the stratosphere. This was soon 
followed by a 2-D simulation [Gray and Pyle, 1989] ena- 
bling the modeling of the latitudinal structure and interac- 
tion with the annual cycle, and subsequently by 3-D simula- 
tions [Hess and O’Sullivan , 1995] which included a better 
representation of the wave-driven transport. Subsequent 
studies employing both 2-D and 3-D models have increased 
our understanding of the mechanisms of the ozone QBO, 
and these are described in more detail in the following sec- 
tions. 

In their 2-D QBO simulation, Gray and Chipperfield 
[1990] also noted QBOs in many of the other trace gases 
carried in the model, some of which were substantial. Some 
of these model predictions were confirmed by the subse- 
quent analysis of SAGE II N0 2 measurements [Zawodny 
and McCormick , 1991] and, more recently, by measure- 
ments of CH 4 , H 2 0, HF, HC1 and NO from the Halogen 
Occultation Experiment (HALOE - e.g. [Luo et al. , 1997, 
Ruth et al. , 1997, Randel et al, 1998, Dunkerton , 1999]). 
Additionally, there is a well-documented modulation of the 
distribution of volcanic aerosol in the lower stratosphere 
(e.g. [Trepte and Hitchman, 1992, Hitchman et al, 1994, 
Grant et al, 1996, Choi et al, 1998]) which has been asso- 
ciated with the QBO. 

5.1.1. Ozone: the Equatorial Anomaly 

The close association of variations in equatorial column 
ozone with the zonal wind QBO is illustrated by Figure 


25(a) which shows a time-series of SBUV and SBUV/2 
equatorial ozone anomalies together with a reference QBO 
wind time series. The signal varies between ±10 DU, ap- 
proximately ±4% of the background total column amount. 
The reference QBO wind time series was calculated by mul- 
tiplying the observed equatorial stratospheric wind profile 
with the weighting profile shown in Figure 25b. The latter 
profile was derived empirically to optimize the fit to column 
ozone, neglecting volcanic periods. Note the excellent cor- 
respondence between the observed ozone anomaly and the 
reference wind series. Positive anomalies are present when 
the zonal winds in the lower stratosphere are westerly, while 
negative anomalies correspond to easterlies. The strongest 
correlations with ozone are achieved with the weighting 
profile biased towards the winds around 20-30 hPa winds 
rather than the 40-50 hPa reference wind normally used in 
correlations with the extratropical zonal wind [Dunkerton 
and Baldwin, 1991]. 

The variable period of the equatorial ozone QBO is 
clearly evident in Figure 26, which shows the latitude time- 
series of the column ozone anomaly. Note that the equatorial 
QBO signal is not synchronized with the annual cycle since 
there is no apparent preferred season in which the anomalies 
change sign or reach their maximum amplitude. 

A mechanism to explain the equatorial column ozone 
QBO anomaly was first suggested by Reed [1964]. The tim- 
ing of the maximum westerly vertical wind shear at a par- 
ticular level corresponds to the warmest phase of the tem- 
perature QBO on the equator. This is therefore the time of 
maximum diabatic cooling, which will induce relative sink- 
ing of air parcels through isentropic surfaces. This vertical 
motion occurs in a region of the atmosphere where the 
ozone mixing ratio increases with height and where the life- 
time of ozone is rapidly changing. Below about 28 km the 
chemical lifetime is relatively long compared with dynami- 
cal processes and ozone may be considered to be a long- 
lived tracer. Above 28 km its chemical lifetime shortens 
considerably. The relative descent of air through this region 
produces an increase in the total column of ozone, since at 
levels above 28 km ozone is replaced by chemical produc- 
tion on relatively short timescales. Thus, the maximum col- 
umn ozone occurs when the column has been displaced far- 
thest downwards into the lower stratosphere. This will be 
after the descent of the westerly shear zone i.e. at around the 
time of maximum westerlies in the lower stratosphere. The 
converse is true of an easterly shear situation. Mass continu- 
ity also requires there to be a return arm to this circulation in 
the subtropics with upwelling associated with westerly 
equatorial shear, and downwelling associated with easterly 
equatorial shear. 

While this conventionally accepted mechanism accounts 
for a large component of the variability, there are neverthe- 
less additional important factors contributing to the column 
anomaly. Figure 27 shows the height time-series of the QBO 
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ozone density anomaly (DU/km) from the SAGE II data set. 
A regression analysis has been applied in order to isolate the 
QBO variability. Ozone density can be used to visually de- 
termine contributions to the column ozone anomaly (simply 
a sum in the vertical of the ozone density anomalies). The 
QBO dominates the ozone variability at the equator, with 
alternating positive and negative anomalies which propagate 
downwards with time [Zawodny and McCormick, 1991, 
Hasebe, 1994, Randel and Wu , 19961. There are two regions 
of maximum ozone perturbation: in the lower stratosphere 
(20-27 km) and the middle stratosphere (30-37km). The 
anomalies at these two levels are approximately a quarter 
cycle out of phase. There is, therefore, a small contribution 
to the column from the region above 28 km, which also in- 
fluences the timing of the maximum column ozone anomaly. 
The equatorial QBO column ozone signal is simulated rea- 
sonably well by models (e.g. [Ling and London , 1986, Gray 
and Pyle , 1989, Tung and Yang , 1994, Chipperfield et al , 
1992, 1994, Kinnersley and Tung , 1998, Jones et al, 1998, 
Hess and O'Sullivan, 1995, Nagashima et al, 1998]), par- 
ticularly when the observed winds are used to force a realis- 
tic zonal wind QBO period [Gray and Ruth, 1993] and 
when the effect of the ozone anomaly itself on the heating 
rate is included [Hasebe, 1994, Li et al, 1995, Huang, 
1996]. 

Although the origin of the lower level anomaly below 28 
km in Figure 27 is conventionally assumed to be due to the 
vertical advection of ozone as proposed by Reed and de- 
scribed above, there is evidence that vertical advection may 
not be the only contribution to the lower level anomaly. If 
only vertical advection were responsible for the ozone 
anomalies one would expect a quadrature relationship be- 
tween the anomalies in ozone and vertical velocity and 
hence also with temperature (since perturbations in tempera- 
ture and vertical velocity are inversely proportional). Obser- 
vations of ozone and temperature, however, indicate that 
they are approximately in phase; Jones et al. [1999] have 
suggested that horizontal, as well as vertical, advection by 
the induced meridional circulation may contribute to the 
equatorial ozone QBO, thus explaining the lack of quadra- 
ture. 

The ozone QBO anomaly above 28 km has been shown 
to be controlled by changes in the photochemical sources 
and sinks of ozone, primarily via transport-induced varia- 
tions of NO y (the total reactive nitrogen reservoir) [Chipper- 
field et al, 1992, 1994, Politowicz and Hitchman, 1997, 
Jones et al, 1998]. 

5.1.2. Ozone: Subtropical and Higher Latitudes 

A QBO signal in the subtropics that extends to mid and 
high latitudes is clearly evident in Figure 26. There is a 180° 
phase change at around 15° in each hemisphere with the 
higher latitude anomaly extending at least to 60°N but with 
its maximum at approximately 30-40° latitude. Broadly, this 


concurs with the presence of a return arm of the local equa- 
torial QBO circulation described above, with ascent (de- 
scent) in the subtropics associated with westerly (easterly) 
equatorial shear. However, there are two significant depar- 
tures from the signature one would expect from this simple 
explanation. Firstly, the theoretical equatorial QBO circula- 
tion is confined to low latitudes and cannot explain the pres- 
ence of an ozone QBO signal poleward of about 30°. Sec- 
ondly, the timing of the subtropical anomalies is such that 
they are not symmetric about the equator. The subtropical 
and higher latitude anomaly maxima and minima in the two 
hemispheres are approximately six months apart and coin- 
cide with the local late winter/spring. This timing is con- 
firmed in Figure 28 in which we show a regression fit of the 
TOMS column ozone amounts to the 30 hPa Singapore 
winds [Randel and Cobb , 1994]. On average, the subtropi- 
cal regression anomalies reach their maximum in March and 
August in the Northern and Southern Hemispheres respec- 
tively. However, note from Figure 27 that occasionally there 
is a “missed” subtropical anomaly e.g. in 1981, 1986 and 
1991 in the NH and in 1993 in the SH. There is therefore a 
change in the “period” of the column ozone QBO as one 
moves to higher latitudes [Hilsenrath and Schlesinger, 
1981], with a phase relationship between the equatorial and 
subtropical anomalies which is constantly changing and is 
more complicated than a symmetric QBO circulation would 
imply. 

The timing of the subtropical and high latitude anomalies 
in late winter/spring is highly suggestive of a modulation of 
the ozone QBO by the annual cycle [Bowman, 1989, Hamil- 
ton, 1989, and Lait et al, 1989]. Early proposed mecha- 
nisms for the seasonal synchronization of the subtropical 
anomalies centered on a modulation of the low latitude 
ozone anomalies once they had been produced by the classic 
symmetric QBO circulation. For example, Holton [1989] 
proposed that transport of the equatorial anomaly by the 
poleward winter circulation could explain the seasonality, 
while Gray and Dunkerton [1990] suggested that downwel- 
ling during the winter would preserve the subtropical ozone 
anomaly while upwelling in summer would destroy it. On 
the other hand, Hamilton [1989] suggested the possibility of 
a modulation of the seasonally varying eddy transport of 
ozone into the subtropics. However, more recently, Jones et 
al, [1998] and Kinnersley [1999] have shown that there is a 
significant modulation of the meridional circulation induced 
by the QBO itself. This seasonal dependence of the circula- 
tion occurs primarily because of nonlinear horizontal advec- 
tion of zonal momentum in the tropics and subtropics by the 
mean circulation, which is highly asymmetric during sol- 
stice periods. This results in a strongly asymmetric QBO 
circulation in which the winter hemisphere half of the circu- 
lation is substantially reinforced and the summer hemi- 
sphere circulation is weakened (Figure 9). Thus, the classic 
symmetric QBO circulation may not be present except, per- 
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haps, in the very low stratosphere where horizontal advec- 
tion by the mean circulation is thought to be weak and at 
equinox when the mean circulation is at its weakest. The 
asymmetry in subtropical ozone anomalies therefore arises 
primarily through its formation by an asymmetric QBO cir- 
culation rather than by the subsequent disruption of a sym- 
metric ozone pattern. 

The “missed” subtropical and mid-latitude anomalies in 
1981, 1986 and 1991 in the NH and in 1993 in the SH (Fig- 
ure 26b) are thought to be due to the timing of the equatorial 
QBO relative to the annual cycle. The formation of a sig- 
nificant winter subtropical anomaly requires not only a 
strong vertical wind shear at the equator so that a strong 
QBO circulation is induced but also the presence of strong 
background horizontal advection which will strengthen the 
winter side of the induced QBO circulation as described 
above. These conditions need to last for a month or two to 
allow the ozone distribution to respond to the induced circu- 
lation. Hence, if either of these requirements are not present 
for a sufficient length of time, a significant subtropical 
anomaly is unlikely to form in that year. 

Similarly, if the timing and duration of the equatorial 
wind QBO is such that it persists in the same phase for two 
successive winters of either hemisphere, then two anomalies 
of the same sign will occur in successive winters in the sub- 
tropics of that hemisphere. The latter is evident in the SH in 
1983/84 and 1988/89. This phenomenon can be thought of 
as a nonlinear interaction between the annual cycle and the 
QBO cycle [Gray and D unker ton, 1990], which results in a 
low frequency modulation of the amplitude of the subtropi- 
cal and mid-latitude ozone anomaly. Over the period of data 
used in Figure 28 this causes the SH regression anomaly to 
be coincidentally larger on average than the NH anomaly. 
Over a longer time span both anomalies would presumably 
be of similar size. 

A typical latitude-altitude cross-section of the modeled 
QBO in ozone from the Jones et at., [1998] model is shown 
in Figure 29 for winter solstice in the NH with an eastward 
QBO wind maximum at about 26 km in the tropics. In both 
the tropics and subtropics the QBO in ozone consists of two 
maxima centered in the lower and middle stratosphere with 
the tropical and subtropical anomalies approximately 180° 
out of phase, in agreement with observations. Also in 
agreement with observations, the modeled anomalies are 
large in the winter hemisphere and small in the summer 
hemisphere in both the photochemically and dynamically 
controlled regions. This is due to the asymmetry in the 
QBO-induced circulation as illustrated by the corresponding 
mass streamfunction also shown in Figure 29. The asym- 
metric ozone anomaly arises directly through advection of 
ozone at the lower levels and indirectly though the advec- 
tion of NO y at the upper levels. 

The Jones et ai [1998], model did not include extratropi- 
cal influences such as the QBO modulation of planetary 


wave breaking. While it produces a good simulation in the 
tropics and subtropics, the modeled ozone anomalies do not 
extend as far poleward as suggested by observations. The 
upper level anomalies inferred from observations maximize 
between 10°-40°, whereas those in the lower stratosphere 
extend to at least 60° [Randel and Wu , 1996], The mecha- 
nism for the poleward extension of the dynamically con- 
trolled ozone anomaly is not well understood, although it 
likely involves an interaction between the planetary waves 
and the equatorial QBO. The modulation of planetary wave 
forcing by the equatorial wind QBO results in a stronger 
large-scale mean circulation in easterly phase years. 
Stronger downwelling in the winter mid-latitudes will pro- 
duce a relatively larger column ozone anomaly in easterly 
years than westerly years, as observed [Tung and Yang , 
1994]. On the other hand, the extension of the ozone anoma- 
lies to mid and high latitudes in the models of Gray and 
Pyle [1989] and Kinnersley and Tung [1998] is a result of 
the seasonally varying eddy transports in their models, 
which transfer the subtropical anomalies to higher latitudes. 
Figure 30 shows the correlation between the simulated 
ozone anomaly of Kinnersley and Tang [1998] and the ob- 
served anomaly between 1978 and 1993. At mid latitudes, 
therefore, a number of factors and feedback processes are 
found to contribute to the final ozone QBO and the QBO 
may actually be responsible for a larger part of the observed 
ozone anomaly than a simple correlation with the wind 
QBO would suggest. 

In addition to the QBO signals in the tropics and mid- 
latitudes, analysis of TOMS data and other long-term re- 
cords suggest a further region of QBO influence in the win- 
ter polar regions [Oltmans and London , 1982, Garcia and 
Solomon , 1987, Bowman , 1989, Lait et ai , 1989, Randel 
and Cobb , 1994]. The polar ozone QBO is approximately in 
phase with mid-latitudes and is seasonally synchronized in 
the same way, with maximum amplitude in springtime. Ob- 
servational evidence for the polar ozone QBO is less statis- 
tically significant than that in the tropics or midlatitudes, at 
least partly because of the high level of interannual variabil- 
ity in the springtime vortex associated with planetary wave 
forcing from the troposphere [Kinnersley and Tung , 1998]. 
There is also a suggestion of a feedback loop between the 
QBO in temperature, the formation of polar stratospheric 
clouds and hence with the underlying chemical destruction 
that gives rise to the Ozone Hole [Poole et ai, 1989, 
Manzini et ai, 1991, Bute hart and Austin , 1996]. 

5. 1 .3. QBO Anomalies in Other Trace Species 

The QBO influences a wide range of other trace gases in 
the atmosphere, including methane, water vapor, volcanic 
aerosol and many of the shorter lived species such as N0 2 
and N 2 0 5 . Volcanic aerosol distributions following major 
eruptions in equatorial latitudes have graphically illustrated 
the different circulation patterns associated with the two 
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QBO phases, as shown in Figure 31 from Trepte and Hitch- 
man [1992]. In a descending westerly phase (Figure 31a) the 
aerosol shows a distinctive “double peak” with relative 
maxima in the subtropics and a minimum at the equator in 
the region 20-50 hPa structure compared with the descend- 
ing easterly phase (Figure 31b) which has a single equato- 
rial peak. The heavy arrows indicate the approximate sense 
of the induced QBO circulations. The near-symmetry of the 
double peak in the aerosol distribution, in contrast to the 
emphasis on hemispheric asymmetry in the previous sec- 
tion, is probably due to being so low down in the equatorial 
stratosphere and hence relatively unaffected by the influence 
of the asymmetrical mean circulation. 

Figure 32 shows the interannual anomalies in H 2 0 over 
the equator from HALOE observations [Randel et al., 
1998]. The QBO anomaly in H 2 0 ascends slowly with time, 
at approximately the rate of ascent of the background mean 
upwelling, in contrast to the slow descent of the ozone 
anomaly. In general the global variations in H 2 0 anomaly 
mimic those observed in CH 4 , but with opposite sign. In the 
upper stratosphere (35-45 km) these are approximately in- 
phase with the QBO winds near 30 km. The CH 4 and H 2 0 
variations over this region cancel to a large degree, such that 
there are much smaller variations in H 2 0+2*CH 4 (the vari- 
able part of total hydrogen); this cancellation confirms the 
production of water vapor via methane oxidation at these 
levels and is evidence that these stratospheric anomalies 
arise via variations in transport. 

One region where differences between CH 4 and H 2 0 
anomalies are observed is in the tropical middle strato- 
sphere, where H 2 0 (and H 2 0+2*CH 4 ) show QBO variations 
over approximately 25—35 km, but none are observed for 
CH 4 . These patterns are equatorially centered (approx. 
15°N-S) and are highly correlated with the QBO zonal 
winds near 20 hPa. One possible mechanism for the tropical 
FI 2 0 signal is that the QBO temperature variations modulate 
tropical tropopause temperatures (e.g., Reid and Gage 
[1985]) and hence the magnitude of water vapor entering the 
lower stratosphere [Mote et al., 1996]. However, the 
HALOE H 2 0 anomalies do not show strong coherence be- 
tween 100 and 30 hPa levels in Figure 32; and the variations 
at the lowest levels do not exhibit a strong QBO behavior. 
The mechanism(s) of the middle stratosphere QBO signal in 
H 2 0 await(s) clarification. 

The measurements of long-lived tracers also show the 
large effects of the QBO over midlatitudes. Figure 33 shows 
the January and April CH 4 distributions from HALOE in 

1993 (westerly phase) and 1994 (easterly phase). In January 

1994 the isolines form a distinctive “staircase” pattern be- 
tween the tropics and NH midlatitudes that is very different 
from the January 1993 behavior; this structure is also seen 
in MLS and CLAES data [Dunkerton and O’Sullivan, 1996, 
O'Sullivan and Dunkerton, 1997, Gray and Russell, 1999, 
Gray, 2000]. Gray [2000] has shown that this large asym- 


metry in QBO anomalies between the hemispheres is con- 
sistent with the asymmetric QBO circulation of Jones et al., 
[1998] and Kinnersley [1999] discussed earlier. The detailed 
staircase pattern in Figure 33c nevertheless results from a 
complicated interaction of advection by the local QBO cir- 
culation and the effects of isentropic mixing at midlatitudes. 

In April 1993 (Figure 33b) the distribution displays a dis- 
tinct “double peak” feature near 0.3-5 hPa, significantly 
higher than that displayed by the volcanic aerosol in Figure 
31. This double peak results from vertical advection by the 
circulation associated with westerly shear of the semiannual 
oscillation {SAO) [Gray and Pyle, 1986, Sassi et al., 1997]. 
The HALOE observations show a distinct QBO variation in 
the amplitude of this SAO double peak, with a prominent 
double peak feature in westerly phase years (e.g. Figure 
33b) but a barely discernible one in easterly phase years 
(e.g. Figure 33d) [Ruth et al. 1997, Randel et al. 1998], This 
is counter-intuitive, since during an easterly QBO phase 
there is enhanced vertical propagation of eastward propagat- 
ing waves and hence a stronger SAO westerly wind shear. 
However, in a 2-D modeling study, Kennaugh et al., [1997] 
showed that this increased eastward wave forcing also 
caused the westerly SAO phase to descend much more rap- 
idly during an easterly QBO phase. Hence, the SAO circula- 
tion that produces the double peak does not remain at any 
one level for long enough for the tracers to respond to its 
presence. 

In the January CH 4 distributions, there is also a distinct 
QBO signal in the steepness of the isolines in the subtropics 
at around 30 hPa. In 1993 (Figure 33a) they slope gently 
from equator to mid-latitudes but in 1994 (Figure 33c) the 
isolines are almost vertical [Gray and Russell, 1999]. This 
difference is not so apparent in the aerosol distributions of 
Figure 31, probably due to the unusual distribution of tracer 
gradients immediately after a volcanic eruption, but it is 
more apparent in composite aerosol measurements [Hitch- 
man et al., 1994], Again, it is not clear whether this is due to 
advection by the QBO circulation or to the QBO influence 
on the equatorward extent of Rossby wave mixing and 
hence on the sharpness of the PV and tracer gradients at the 
subtropical edge of the surf zone. Some studies have sug- 
gested a sensitivity of isentropic mixing to the QBO [Dun- 
kerton and Baldwin, 1991, O'Sullivan and Young, 1992; 
O'Sullivan and Chen, 1996; O'Sullivan and Dunkerton, 
1997] and there is evidence of this in satellite measurements 
[Grant et al., 1996]. On the other hand, Waugh [1996] used 
analyzed winds and contour advection techniques to make 
quantitative estimates of mixing in the subtropics, but found 
no sensitivity to the QBO. This may be due to a lack of reli- 
able wind data in the subtropics. Gray and Russell [1999] 
have pointed out that the strong QBO signature is in the 
steepness of the isolines and not in the isentropic gradients, 
which suggests that advection by the QBO circulation is 
important in setting up this feature. This is supported by the 
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modeling of Jones et al„ [1998] and Kinnersley [1999] 
which reproduce some aspects of this steepening using ad- 
vective transport only. 

The QBO anomalies in the shorter lived species N0 2 , as 
measured by HALOE, are shown in Figure 34. The general 
pattern of the anomalies in NO are similar to those of N0 2 
and reflect QBO-induced changes in the abundance of NO y 
in the stratosphere [Gray and Chipperfield , 1990; Jones et 
a/., 1998]; the anomalies in NO y are produced through the 
influence of vertical advection on the vertical gradient in 
NO y ( Chipperfield et ai, 1994; Politowicz and Hitchman , 
1997; Jones et ai , 1998]. Below about 5 hPa, the amplitude 
of the QBO in N0 2 is larger than in NO as a result of the 
influence of the QBO anomalies in temperature on the ratio 
of the abundance of NO to N0 2 , through the reaction NO + 
O^ - > N0 2 + 0 2 [Gray and Chipperfield , 1990]. The reason 
for a large QBO signal in NO above 5 hPa has not yet been 
investigated. However, it is likely due to the dominant con- 
tribution of NO to NOy at these altitudes. 

The pattern of QBO anomaly in HC1 observed by 
HALOE (not shown) is similar to that in NO and N0 2 . In 
their modeling study Gray and Chipperfield [1990] found a 
similarity between the QBO anomaly in NO y and Cl y (the 
total inorganic chlorine reservoir). The anomalies were pro- 
duced by the influence of vertical advection on the positive 
vertical gradients in Cl y and NO y in the lower and middle 
stratosphere. In the upper stratosphere, where the vertical 
gradient in Cl y is weak, the QBO anomalies in HCI are also 
weak, providing supporting evidence for this mechanism. 


6. THE QBO ABOVE AND BELOW THE 
STRATOSPHERE 

6. 1 . Mesospheric QBO 

Rawinsonde measurements of equatorial winds (to -30 
km, since the 1950s) and rocketsonde wind observations 
from stations near 8°S and 8°N (to -60 km, since the 1960s) 
(as shown in Figure 1). Satellite measurements of equatorial 
winds in the stratosphere and mesosphere from the High 
Resolution Doppler Imager (HRDI) on the Upper Atmos- 
phere Research Satellite (UARS), beginning in November, 
1991, provide equatorial wind data from 10-40 and 50-115 
km. The HRDI data allowed the discovery of a QBO in the 
upper mesosphere [Barrage et ai , 1996], called the meso- 
spheric QBO (MQBO). Monthly mean HRDI equatorial 
winds are shown in Figure 35 (top panel) confirming the 
QBO up to 40 km, and the mesospheric SAO from -55-85 
km. By removing the annual and semi-annual harmonics 
(similar to removing the seasonal cycle in the top panel of 
Figure 1) a mesospheric QBO centered near 85 km becomes 
apparent. The wind variations have been confirmed by radar 
observations at Christmas Island (2°N) during the HRDI 
period. The HRDI data show that the MQBO extends out to 


±30° latitude, with a 180° phase difference with the strato- 
spheric QBO (-40 hPa). 

The HRDI data record is too short to reliably confirm that 
the MQBO is linked to the stratospheric QBO. Garcia et ai 
[1997] suggested that the easterly phase of the mesospheric 
SAO is usually stronger only when deep easterlies are pre- 
sent. This relationship holds during 1992-1995, but excep- 
tions are found when Christmas Island radar data from 
1990-1991 are examined. Mesospheric westerlies do not 
show marked interannual variability, and they are not corre- 
lated with the QBO. 

The possibility of a connection to the QBO is strength- 
ened by both modeling and theoretical evidence. Mayr et ai 
[1997] used a 2-D model to simulate oscillations in the 
equatorial stratosphere and mesosphere resulting from verti- 
cally propagating gravity waves. The modeled QBO was not 
confined to the stratosphere, but showed a QBO in the upper 
mesosphere, similar to that observed in HRDI data and 
Christmas Island radar data. The theoretical explanation 
involves selective critical layer absorption or wave filtering 
of small-scale gravity waves as they traverse the underlying 
winds in the stratosphere, together with the complimentary 
wave breaking at higher altitude levels in the upper meso- 
sphere. This process also generates the SAO in the upper 
mesosphere. 

The amplitudes of the various equatorial oscillations in 
zonal wind, as a function of height, are summarized in Fig- 
ure 36. The annual cycle (dotted line) is relatively small in 
the stratosphere (-5 ms* 1 ). The stratospheric QBO is shown 
from 16-40 km, peaking near 20 ms 1 at about 25 km. The 
amplitude in the troposphere is small. The amplitude of the 
QBO between 40 and 70 km is not shown due to uncertainty 
and difficulty defining what part of the variability is related 
to the QBO (see Figure 1 ). 

6.2. Effect of the QBO on the extratropical 
troposphere 

In Section 4 it was shown that the QBO, by modulating 
the wave guide for vertically propagating planetary waves, 
affects the circulation of the extratropical winter strato- 
sphere. This modulation is more easily seen in the NH, 
where wave amplitudes are larger and the stratospheric cir- 
culation is disrupted by sudden warmings in which the usual 
westerly middle stratospheric winds become easterly. Figure 
19 showed that modulation of the zonal wind by the QBO in 
the NH in January appears to extend below the tropopause. 
Angell and Korshover [1975] showed a strong correlation 
between Balboa 50-hPa zonal wind and the displacement of 
the northern vortex at 300 hPa, near the tropopause. The 
surface signature of the QBO was first examined by Holton 
and Tan [1980], who showed the difference between 1000- 
hPa geopotential for the two phases of the QBO. An update 
of Holton and Tan’s calculation, for 1964-1996 data, is 
shown as Figure 37. The pattern is characterized by modula- 
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tion of the strength of the polar vortex, and anomalies of the 
opposite sign at low- to mid-latitudes. The pattern in Figure 
37 is essentially similar to that shown by Holton and Tan. 
Hamilton [1998], in a 48-year GCM simulation with an im- 
posed QBO, found that the difference, between QBO 
phases, of the strength of the upper tropospheric polar vor- 
tex was statistically significant. 

Disturbances in the extratropical stratosphere affect only 
-25% of the vertical column of air, and there is increasing 
evidence, from observations, numerical models, and concep- 
tual models, that stratospheric anomalies do influence the 
troposphere. It is not necessary to limit the discussion to the 
influence of the QBO, but to think of any circulation anom- 
aly in the stratosphere (e.g. due to solar influence, the QBO 
phase, a volcanic eruption). Boville [1984] showed, using a 
GCM, that a change in the high-latitude zonal wind struc- 
ture introduced changes in the zonal-mean flow all the way 
down to the boundary condition that corresponds to the 
Earth’s surface, as well as in wave structures. He concluded 
that the degree of trapping of the planetary waves in the 
troposphere is determined by the strength and structure of 
the stratospheric mean zonal wind, resulting in sensitivity ot 
the troposphere to the stratospheric zonal wind structure. 
Boville [1986] explained further that when the high-latitude 
winds, especially in the lower stratosphere are strong, it 
tends to inhibit the vertical propagation of wave activity into 
the polar stratosphere. If the winds are weak, wave activity 
can propagate more effectively into the polar stratosphere. 
The process was found, in a GCM, to be tightly coupled to 
tropospheric generation of vertically-propagating planetary 
waves. 

Kodera et al. [1990] used both observations and GCM 
output, to show that anomalies in the mid-latitude upper 
stratosphere (1 hPa) in December tend to move poleward 
and downward, reaching the troposphere approximately two 
months later. In general, these effects can be understood in 
terms of the modification of the zonal-mean zonal wind, 
which acts as a wave guide for planetary wave propagation. 
Stratospheric anomalies tend to induce changes in wave 
propagation at lower levels, which affect the convergence of 
the waves, which further modifies the zonal-mean flow. 
Over time, the net effect appears to be downward (and 
poleward) movement of the anomalies. 

A complementary approach to understanding the down- 
ward link to the troposphere has been to examine “modes of 
variability.” Such modes may be thought of as patterns 
which tend to recur and which account for a large traction of 
variability; patterns should be robust and be found through 
different analysis schemes. For example, the NH winter 
zonal wind tends to vary in a dipole pattern (e.g., Figure 
19). This coupled mode of variability between the northern 
winter stratosphere and troposphere was discussed by Ni- 
gam [1990], who examined rotated empirical orthogonal 
functions (EOFs) of zonal-mean wind. Nigam’s result 


showed that the dominant mode of variability in zonal-mean 
wind appears as a deep north-south dipole, with a node near 
40-45°N (similar to Figure 19). The poleward part of the 
dipole represents fluctuations in the strength of the polar 
vortex. 

Coupling between the stratosphere and troposphere was 
further explored by Baldwin et al. [1994], who examined 
geopotential patterns in the middle troposphere that were 
linked to the stratosphere. Using singular value decomposi- 
tion (also called maximum covariance analysis) between 
500 hPa geopotential and mean zonal wind, they showed 
that the leading mode had a strong dipole signature in mean 
zonal wind, extending from the surface to above 10 hPa. 
The north-south dipole mode accounts for a large fraction of 
the variance in zonal wind, and is found by a variety of tech- 
niques. 

The leading mode of variability of the northern ex- 
tratropical troposphere/stratosphere is characterized by a 
deep, zonally-symmetric or “annular” structure [Thompson 
and Wallace, 1999]. This dipole mode in zonal-mean zonal 
wind is coupled to a horizontal wave structure of geopoten- 
tial anomalies in the troposphere. The surface pattern re- 
sembles the North Atlantic Oscillation, but is more symmet- 
ric in longitude. Thompson and Wallace [1998] showed that 
the surface pattern corresponds to the leading EOF of win- 
tertime monthly-mean sea level pressure. The entire mode, 
throughout the troposphere and stratosphere, is known as 
the Northern Annular Mode (NAM). The surface NAM pat- 
tern is also know as the Arctic Oscillation [Thompson and 
Wallace, 1998], and is broadly similar to the QBO signature 
shown in Figure 37, suggesting that the QBO may act to 
modulate the deep structure of the NAM. It is now becom- 
ing clear that all of the studies of modes of NH variability 
produce patterns that are essentially slight variants of the 
NAM. The NAM represents a dominant, robust, naturally 
occurring mode of variability, and if the QBO phase can 
affect the NAM in the stratosphere, it can be expected that 
there would be a surface signature of the QBO. 

The NAM is closely linked to stratospheric sudden warm- 
ings [Baldwin and Dunkerton, 1999], and every major 
warming shows a clear signature in the magnitude of the 
NAM. This relationship can be expected because both phe- 
nomena, in the stratosphere, relate to the strength of the po- 
lar vortex. As the strength of the stratospheric polar vortex 
changes, the surface NAM signature tends to vary. Baldwin 
and Dunkerton examined this relationship and demonstrated 
that large, sustained variations in the strength ol the strato- 
spheric polar vortex tend to propagate downward to the 
earth’s surface. The time for propagation from 10 hPa to the 
surface was found to be variable, but averaged about three 
weeks. They also examined the relationship between the 
QBO and the NAM, which was found to be strongest during 
December in the middle stratosphere, and weaker as winter 
progressed. The QBO appears to be a factor that influences 
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the NAM. That influence tends to appear as modulation of 
the strength of the polar vortex, from the lower mesosphere 
to the Earth’s surface [Baldwin and Dunkerton , 2000]. 

6.3. Effects of the QBO on the tropical 
troposphere 

Because the QBO has its maximum amplitude over the 
equator, it is logical to inquire whether this oscillation has 
any effect on the underlying tropical troposphere. In this 
regard, it is important to keep two things in mind. First, the 
zonal wind and temperature anomalies of the QBO do not 
penetrate significantly below the tropopause (Figure 1). The 
temperature QBO at the tropopause is small relative to the 
annual cycle. Second, it is known that the tropical tropo- 
sphere has a quasi-biennial oscillation of its own, uncorre- 
latcd with the stratospheric QBO [ Yasunari , 1985; Gutzler 
and Harrison, 1987; Kawamura, 1988; Lau and Sheu, 1988; 
Moron et al. , 1995; Shen and Lau , 1995]. Unlike the latter, 
the “tropospheric QBO” is irregular in time, asymmetric in 
longitude, and propagates slowly eastward, with largest am- 
plitude near the Maritime Continent. 

Although Yasunari [1989] suggested that the tropospheric 
oscillation is coherent with the stratospheric QBO, his re- 
sults and those of other authors tend to disprove the claim. 
For example, the Hovmoller plot of biennially filtered upper 
tropospheric winds has an irregular variation in longitude 
and time, with at least two apparently distinct oscillations in 
the Pacific and Atlantic sectors, neither of which correlate 
well with the stratospheric QBO. The oscillation is a bit too 
fast over the Atlantic, and bit to slow over the Pacific. On 
the basis of a longer record, some authors view the strato- 
spheric and tropospheric QBOs as completely unrelated 
[ Barnett , 1991; Xu , 1992]. As far as linear correlations are 
concerned, tropospheric and stratospheric QBOs do not ex- 
hibit a consistent phase relationship over several decades. 
Their morphologies are so different that it is difficult to see 
any obvious connection. 

On the other hand, a more subtle relationship (either 
nonlinear or multivariate) might exist between these phe- 
nomena. There is evidence that ENSO warm events affect 
the rate ot QBO westerly descent [Maruyama and Tsuneoka, 
1988]. This effect is dynamically plausible [Dunkerton, 
1990; Geller et al. , 1997] but would not result in any linear 
correlation. Perhaps in a similar way, the stratospheric QBO 
influences the underlying troposphere, its effect mixed with 
that of other phenomena or constrained to operate nonline- 
arly only at certain times and places. 

The most promising connection between the stratospheric 
QBO and tropical troposphere is found in the interannual 
variation of Atlantic hurricane activity [Gray, 1984a, b; 
Shapiro, 1989; Hess and Eisner, 1994; Landsea et al., 1998; 
Eisner et al, 1999]. Strong hurricanes originating in the 
tropical Atlantic occur significantly more frequently in sea- 
sons when the overlying QBO is westerly or becoming 


westerly near 50 hPa. The reverse is true in the opposite 
phase of the QBO. The stratospheric QBO remains one of 
several predictors of Atlantic hurricane activity in seasonal 
forecasts issued by W. Gray and collaborators at Colorado 
State University: http://typhoon.atmos.colostate.edu/ 

forecasts/. 

It is unclear whether the QBO has a similar influence on 
typhoons in the western Pacific [Chan, 1995; Baik and 
Paek, 1998; Lander and Guard, 1998]. The dynamics of 
hurricane formation are somewhat different in the two oce- 
anic regions. 

A convincing explanation of the QBO's effect on hurri- 
cane activity has not been given. Suggestions were made 
concerning the effect of lower stratospheric vertical wind 
shear on penetrative convection associated with strong 
storms [Gray et al, 1992a, b], the effect on lower strato- 
spheric static stability [Knaff, 1993], and the effect of QBO 
winds on the position of critical levels for tropical easterly 
waves [Shapiro, 1989]. Evidence supporting a role of the 
QBO in hurricane activity is derived from multiple regres- 
sion in which the predictors are chosen subjectively from 
experience. The possibility that other forms of quasi- 
biennial variability might equally well explain the Atlantic 
hurricane connection has not been explored [Shapiro, 1989]. 

Other apparent effects of the QBO in the troposphere in- 
clude the remarkable finding of Chao [1989] that the Earth’s 
length of day has an interannual variation coherent with the 
stratospheric QBO’s angular momentum. This result reflects 
the fact that atmospheric angular momentum is intimately 
connected to the rotation rate of the Earth. In a similar vein, 
Del Rio and Cazenave [1994] discuss a possible effect on 
polar motion. Fontaine et al [1995] found that contrasting 
precipitation regimes in west Africa are associated with the 
stratospheric QBO. Collimore et al [1998] showed a corre- 
lation, albeit imperfect, between the QBO and deep convec- 
tive activity in regions of strong convection. In the realm of 
much smaller signals, Hamilton [1983] found quasi-biennial 
variability in the amplitude of the semidiurnal surface pres- 
sure oscillation. These and other published and unpublished 
evidences of the QBO's effect on the tropical troposphere 
remain tantalizing enough for further study, and demonstrate 
that the stratospheric QBO should be properly simulated in 
models of the tropical atmosphere. 

7. CONCLUSIONS 

In a paper summarizing work on the then recently dis- 
covered QBO, Reed [1967] stated: 

“Perhaps a simple explanation will soon be found, 
and what now seems an intriguing mystery will be 
relegated to the category of a meteorological freak. 

Or perhaps the phenomenon will prove to have a 
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greater significance than we now might envisage , ei- 
ther because of some intrinsic property it possesses 
or because of its effect on other related areas of re- 
search. ” 

With the benefit of more than 3 decades of QBO related 
research, we can now say assuredly that the QBO is more 
than a meteorological “freak.” Indeed, as we have tried to 
demonstrate in this review, it does have importance beyond 
that envisaged in the 1960s, both for its inherent fluid dy- 
namical characteristics and for its relevance to issues of 
global atmospheric climate and chemistry. 

The QBO is a spectacular demonstration of the role of 
wave, mean-flow interactions in the fluid dynamics of a 
rotating stratified atmosphere. As has been elegantly argued 
by McIntyre [1993], what makes the dynamics of a rotating 
stratified atmosphere special is the ubiquitous occurrence of 
wave motions, and the fact that wave propagation and wave 
refraction are generally accompanied by a transport of mo- 
mentum. The QBO could not exist if it were not for momen- 
tum transfer by wave propagation and refraction. The de- 
pendence of wave refraction on the mean How provides the 
mechanism whereby wave-induced momentum fluxes in the 
equatorial stratosphere can produce a feedback onto the 
mean flow. In the QBO not only do the oscillating waves 
interact with the mean flow to produce a flow rectification, 
but also the rectified flow itself oscillates on a period com- 
pletely different from that of the driving waves. 

Figure 1 shows that the QBO, which in the 1960s was re- 
garded by some as likely to be merely a transient phenome- 
non, is a very persistent feature of the circulation of the 
equatorial stratosphere. We now have direct observations of 
20 full cycles of the oscillation, and there is indirect evi- 
dence covering a much longer time. By study of long-term 
variations in the solar semidiurnal tidal signal in surface 
pressure at equatorial stations (which is sensitive to the 
zonal winds in the stratosphere) Hamilton [1983] and Teitel- 
baum et al [1995] have argued that the QBO must have 
existed for at least the past 120 years. 

This robust nature of the QBO suggests that similar phe- 
nomena should be present on other planets with rotating 
stratified atmospheres and equatorial convection zones. In- 
deed an analogous oscillation, the quasiquadrennial oscilla- 
tion (QQO), has been documented in the equatorial atmos- 
phere of Jupiter [Leovy et al. t 1991, Friedson, 1999The ob- 
served meridional scale of the QQO on Jupiter ~7° latitude 
is about 1/2 that of the terrestrial QBO. For the parameters 
given by Friedson, this is consistent with the transition scale 
discussed in Section 3.1, provided that the vertical scale of 
the forcing is set at 12 km, rather than the 4-km value ap- 
propriate for the terrestrial QBO. ]. McIntyre [1994] has 
suggested that a similar oscillation may occur in the solar 
interior. 

The possibility of broader implications of the QBO for 


other areas of research, as suggested by Reed in the above 
quotation, has certainly proved to be true. As discussed in 
Section 6,2, the influence of the QBO on interannual cli- 
mate variations in the extratropical troposphere and the 
stratosphere is a major subject of current interest. Attempts 
to better understand and predict trends and variability of 
atmospheric ozone also require careful consideration of di- 
rect and indirect effects of the equatorial QBO on the ozone 
layer (see Section 5). Thus, models of interannual climate 
variability and of global stratospheric chemistry both should 
include the effects of the QBO either explicitly or through 
some parameterization. 

Unfortunately, simulation of the QBO remains a great 
challenge for general circulation models. Many such models 
are currently being used for prediction of climate trends and 
variability associated with human induced changes in the 
concentrations of various greenhouse gases. Yet, as dis- 
cussed in Section 3.3.2, most such models are unable to 
spontaneously generate a realistic QBO. The atmosphere, 
however, has no such difficulties. A skeptic might argue that 
the absence of such a robust global-scale dynamical phe- 
nomenon is evidence that the models are still far from real- 
ity in some important respects. The arguments presented 
above suggest that gravity waves generated by equatorial 
convection are essential to forcing of the QBO. This in turn 
suggests that among other things, better parameterization of 
the dynamics of convective systems is required if GCMs are 
to routinely reproduce the fascinating wave, mean-flow 
feedback interactions that result in the equatorial QBO. 


GLOSSARY 

[3-plane: An approximation of the Coriolis parameter in 
which / = / 0 + (3}\ where beta is a constant. The Corio- 
lis parameter varies linearly in the north-south direc- 
tion. 

Boussinesq limit: A simplification in which density is 
treated as a constant except where it is coupled to 
gravity in the buoyancy term of the vertical momen- 
tum approximation. 

Coriolis Parameter: / = 2Q sin ([), where Q. is the rota- 
tion rate of the earth, and (() is latitude. 

Critical level or line: For a wave propagating on a mean 
flow, the point at which the wave’s phase speed equals 
the mean flow speed in the direction of wave propaga- 
tion. 

Easterly: From the east. 

Eastward: To the east. 
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Eliassen-Palm (EP) flux: The Eliassen-Palm flux [Elias - 
sen and Palm , 1961; Edmon et al , 1980] is a measure 
of the propagation of wave activity in the latitude- 
height plane. The convergence of the EP flux is a 
measure of the (retarding) force exerted by the waves 
on the zonal-mean flow. 

GCM: General Circulation Model 

Gravity waves: Oscillations usually of high frequency 
and short horizontal scale, relative to synoptic-scale 
motions, that arise in a stably stratified fluid when par- 
cels are displaced vertically. 

Inertia-gravity waves: Low-frequency gravity waves 
that are substantially affected by the Coriolis force. 

Kelvin Waves: At the equator, eastward propagating 
waves with negligible meridional velocity component 
and Gaussian latitudinal structure in zonal velocity, 
geopotential, and temperature, symmetric about the 
equator. 

Log-pressure height: A vertical coordinate which is 
proportional to the logarithm of pressure. It is ap- 
proximately equal to physical height. 

Meridional plane: The latitude-height plane. 

Mixed Rossby-gravity waves: Westward propagating 
waves with meridional velocity maximum on the equa- 
tor, symmetric about the equator. Zonal wind, geopo- 
tential, and temperature are antisymmetric. 

MLS: Microwave Limb Sounder 

MQBO: Mesospheric quasi-biennial oscillation 

NAM: Northern Annular Mode. 

NCEP: National Centers for Environmental Prediction 
(Washington, DC, USA) 

Phase of the QBO: Easterly or westerly as defined by 
the equatorial winds at a precise level. Historically, the 
level has been 50 to 30 hPa, with 40 hPa being typical. 

Planetary-scale waves: Tropical or extratropical distur- 
bances with low zonal wave number (1-3), e.g., equa- 
torial Kelvin waves or Rossby waves in the winter 
stratosphere. 

Planetary wave breaking: Rossby waves owe their exis- 
tence to meridional gradients in potential vorticity, 
which provide a restoring force to allow propagation 
of the waves. In the wintertime stratosphere, the polar 
vortex provides steep meridional gradients of PV, sur- 
rounded by smalt gradients. At the edge of the vortex, 
there is overturning (in the meridional direction of PV 
by Rossby waves as their amplitudes become large. 
This “wave breaking” results in long, drawn-out 


tongues of PV and irreversible mixing of PV at small 
scales. 

Rossby waves: See planetary-scale waves. 

SAO: Semi-annual oscillation, with a period of six 
months. 

Surf zone: See Rossby wave breaking. 

Wave guide: The path of propagation of these waves is 
analogous to the propagation of light through a me- 
dium with a variable index of refraction. For planetary- 
scale waves, the index of refraction depends on the 
zonaliy-averaged zonal wind. 

Westerly: From the west. 

Westward: To the west. 

Zonal: Longitudinal. The zonal wind is positive from the 
west. 
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Figure 1 . Time-height section (top) of the monthly-mean zonal wind component (ms 1 ), with the seasonal cycle 
removed, for 1964-1990. Below 31 km equatorial radiosonde data are used: Canton Island (2.8°N, 1964-August 
1967), Gan/Maledive Islands (0.7°S, September 1 967-December 1975), and Singapore ( 1 ,4°N, January 1976-1990). 
Above 31 km rocketsonde data from Kwajalein (8.4°N) and Ascension Island (8.6°S) are shown. The contour- 
interval is 6 ins' 1 with the band between -3 and +3 unshaded. Red represents positive (westerly) winds. After Gray ct 
al. , 2000. In the bottom panel the data are bandpass filtered to retain periods between 9 and 48 months. 
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Figure 1. Time-height section (lop) of the monthly-mean zonal wind component (ms 1 ), with the seasonal cycle 
removed, for 1964-1990. Below 31 km equatorial radiosonde data are used: Canton Island (2.8°N, 1964-August 
1967), Gan/Maledive Islands <0.7°S, September 1967-December 1975), and Singapore ( L4°N, January 1 976- J 990). 
Above 31 km rocketsonde data from Kwajalein (8.4°N) and Ascension Island (8.6°S) are shown. The contour 
interval is 6 ms * with the band between -3 and +3 unshaded. Red represents positive (westerly) winds. Alter Gray cl 
al h 2000. in the bottom panel the data arc bandpass Filtered to retain periods between 9 and 48 months. 
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Figure 3. Histograms of the number of transitions (zero crossings) at 50 hPa grouped by month. Individual years 
are listed in the boxes. Easterly to westerly transitions are displayed in the top panel, while westerly to easterly 
transitions are shown in the bottom panel. After Pawson et al„ [1993]. 



Figure 4. Harmonic analysis of 30-hPa zonal wind, showing the amplitude of the annual cycle (squares) 

^ StarS> ’ at ’ d residual deseasoned component (circles). Symbols show individual radiosonde^tation 
amplitudes. Solid lines are based on binned data. From Dunkerton and Delisi [19851. 



Figure 5. Dynamical overview of the QBO during the northern winter. The propagation of various tropical waves 
are depicted by tan arrows, with the QBO driven by upward-propagating gravity, inertia-gravity. Kelvin, and Rossby- 
gravity waves. The propagation of planetary-scale waves (purple arrows) is shown at mid-to high latitudes. Black 
contours indicate the difference in zonal-mean zonal winds between easterly and westerly phases of the QBO wheie 
the QBO phase is defined as easterly by the 40-hPa equatorial wind. Easterly anomalies are light blue and westerly 
anomalies are pink. In the tropics the contours are similar to the observed wind values when the QBO is easterly. The 
mesospheric QBO (MQBO) is shown above ~80 km, while wind contours between -50-80 km are dashed due to 
observational uncertainty. 
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Figure 5. Dynamical overview of the QBO during the northern winter. The propagation of various tropical waves 
arc depicted by tan arrows, with the QBO driven by upward-propagating gravity, inertia-gravity, Kelvin, and Rosshy- 
gravilv waves. The propagation of planetary-scale waves (purple arrows) is shown at mid-lo high latitudes. Black 
contours indicate the difference in zonal-mean zonal winds between easterly and westerly phases of the QBO where 
the QBO phase is defined as easterly by (he 40-hPa equatorial wind. Easterly anomalies arc light blue and westerly 
anomalies are pink. In the tropics the contours are similar to the observed wind values when the QBO is caslcily. The 
mesospheric QBO (MQBO) is shown above -BO km, while wind contours between -50-80 km are dashed due to 
observational uncertainty. 
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Figure 7. Cross sections of QBO anomalies in February 1994. Temperature anomalies are contoured (+/- 0.5, 1.0, 
1.5,...K) and components of the residual mean circulation (v, w) as vectors (scaled by an arbitrary function of 
latitude. From Randel et al. [1999]. 
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Figure 8. QBO regression using 13-22-km MSU temperature data for 1979-1998. Stippling denotes regions 
here the statistical fits are different front zero at the 2 o level. Updated from Rondel and Cobb 1 1994 1. 




Figure 9. Overview of tracer transport by QBO wind anomalies and mean advcction. Contours illustrate 
schematically the isopleths of a conservative tracer during the northern winter when the QBO is in its east phase at 
40 hPa (matching Figure 5). Tropical up welling causes the broad maximum in tracer density in the middle to upper 
equatorial stratosphere, while the QBO causes deviations trom hemispheric symmetry. Red arrows near the equator 
depict circulation anomalies from the QBO. The circulation anomaly in the equatorial lower stratosphere is 
approximately symmetric, while the anomaly in the upper stratosphere is much stronger in the winter hemisphere. 
The descent near the equator (~5 hPa) and ascent to the north (-5 hPa, 10°N) combine to produce a “stair-step 
pattern. A second stair step is formed in midlalitudes by horizontal mixing. 
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Figure 9. Overview of tracer transport by QBO wind anomalies and mean advcction. Contours illustrate 
schematically the isopieths of a conservative tracer during the northern winter when the QBO is in its east phase at 
40 hPa (matching Figure 5). Tropical upwclling causes the broad maximum in tracer density in the middle to upper 
equatorial stratosphere, while the QBO causes deviations from hemispheric symmetry. Red arrows near the equator 
depict circulation anomalies from lhe QBO, The circulation anomaly in the equatorial lower stratosphere is 
approximately symmetric, while the anomaly in the upper stratosphere is much stronger in the winter hemisphere. 
The descent near the equator (-5 hPa) and ascent to the north (~5 hPa, 1(TN> combine to produce a **siair-sicp*‘ 
pattern. A second stair step is formed in midlatitudcs by horizontal mixing. 
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Figure 11 A timc-heieht section of short-period t< 4 days) northward component (top) and temperature 
component (bottom) at Watukosek. Indonesia (7.6°S. 1 12.7°E) for 24 days in February to March 1990. From Tsuda c, 

al. < 1994b) 






Figure 12. Power spectra for (a) T and (b) u fluctuations at Singapore as a function of time, averaged for a height 
region of 20-25 km. (c) Cospectra and (d) quadrature spectra of T and u components are also shown. The thick solid 
line represents a QBO reference time series. From Sato and Dunkerton [1997]. 
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Figure 13 Momentum flux estimates for short-period (1-3 days) component in the (a) westward and (b) eastward 
shear phases. Left panels show indirect estimates corresponding to a sum of absolute values of positive and negative 
momentum fluxes. Right panels show direct estimates corresponding to net momentum fluxes. From Sato [ 1997]. 




Figure 14 . Geographical distribution of vertical component of Eliassen-Palm flux of equatorial waves with 
periods (a) longer than 5 days, (b) of 2-5 days, and (c) shorter than 2 days estimated based on high-resolution 
imagery of the global convective pattern constructed from six satellites. Shown here is the frequency-integrated 
absolute value of the Eliassen-Palm flux. Units are arbitrary and contour increments are linear. From Bergman and 
Salby[ 1994]. 
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Figure 15. Stratospheric gravity waves above a simulation of tropical convection. The shading shows vertical 
velocities. The scale runs from ± 1.2 ms ' to enhance the wave perturbations in the stratosphere, but maximum values 
in the troposphere exceed ± 5 ms 1 . Potential temperature contours at I0K intervals (thin lines) and the outline ol the 
storm cloud (heavy lines) are also shown [Alexander and Holton . 1997], 
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Figure 16. (a) Gravity wave momentum flux in the stratosphere vs. cloud-top temperature below for oceanic 
flights north of Australia during Jan-Feb 1987. (b) Shows the directional distribution of the momentum fluxes in 
these observations. 
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Figure 17. Schematic latitude-height sections showing the mean meridional circulation associated with the 
equatorial temperature anomaly of the QBO. Solid contours show temperature anomaly isotherms, dashed contours 
are zonal wind isopleths. Plus and minus designate signs of zonal wind accelerations driven by the mean meridional 
circulation, (a) Westerly shear zone, b) easterly shear zone. After Plumb and Bell [1982b]. 







Figure 18. lime-height section of the zonal-mean zonal wind over the equator simulated by Takahaslit 1 1999], 
c time coordinate ranges from day 0 ( l January with a 360-day model year) to day 1830 The contour interval is 6 
; Red and blue shading represents westerly and easterly winds, respectively. 







Figure 18 . Time-height section of the zonal-mean zonal wind over the equator simulated by Takahaxhi [1999], 
The time coordinate ranges from day 0 ( 1 January with a 360-day model year) to day 1 830. The contour interval is 6 
ms' 1 Red and blue shading represents westerly and easterly winds, respectively. 






Figure 19. January latitude-height zonal-mean wind difference between the average of all years with westerly 
QBO and those with easterly QBO. The phase of the QBO is optimized for the NH, and is defined using the EOF 
technique of [Baldwin and Dunkerton, 1998b], and is nearly equivalent to using 40-hPa equatorial winds. The 
contour interval is 2 ms 1 and negative values are shaded. From Baldwin and Dunkerton [1998b]. 
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Figure 20. NH 5-hPa month-latitude zonal-mean wind difference between the average of all years with westerly 
QBO and those with easterly QBO. The phase of the QBO is defined as in Figure 19 and the data are monthly means. 
The contour interval is 2 ms 1 and negative values are shaded. From Baldwin and Dunkerton [1998b]. 
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Figure 21. SH 5 -hPa month-latitude zonal-mean wind difference between the average of all years with westerly 
QBO and those with easterly QBO. The calculation is as in Figure 20, except the phase of the QBO is optimized for 
the SH, which is nearly equivalent to using 25 hPa equatorial winds. From Baldwin and Dunkerton [1998b]. 



Figure 22. Westerly minus easterly phase composite of the zonal-mean zonal wind in January from the 48 year 
GCM integration described in Hamilton [1998]. The composite was made based on the 20 Januarys with most 
westerly equatorial winds at 40 hPa and the 20 Januarys with most easterly equatorial 40 hPa winds. Results are 
shown here to 1 hPa but the model domain actually extends to 0.01 hPa. 
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EQUATORIAL OZONE ANOMALIES 3 , 



YEAR 


Figure 25. (a) Time-series of equatorial ozone anomaly in percent (solid line) from SBUV and SBUV/2 together 
with a reference QBO wind time series (dotted line) compiled by multiplying the observed winds at Singapore by the 
weighting profile shown in (b). From Handel and Wu [1996]. 
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Figure 26. Latitude-time sections of column ozone anomalies from combined SBUV-SBUV/2 data, (a) Full 
anomalies defined as deseasonalized and detrended over 1979-94. and (b) the QBO component derived by seasonally 
varying regression analysis. Data in all panels have been multiplied by cos (latitude) to account for area weighting. 
Contour interval is 3 DU, with zero contours omitted and positive values shaded. Diagonal hatching denotes 
unreliable data. Vertical lines denote January of each year. From Randel and Wu [1996]. 
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Figure 27. Height - time series of interannual anomalies ozone density (DU/km) derived using a regression 
analysis to isolate the QBO variation. Contour intervals are 0.3 DU/km with zero contours omitted and positive 
values shaded. From Randel and Wu [ 1996], 
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Figure 28. Latitude - time section of QBO-associated regression fit of zonal mean column ozone (DU) from the 
TOMS instrument to the 30 hPa Singapore winds for the period 1979-1994. Stippling denotes regions where the 
statistical fits are different from zero at the 2c level. Hatched regions denote the polar night, where no ozone data are 
available. Updated from Randel and Cobb [1994], 
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Figure 29. QBO anomaly in ozone (DU/km) in January from the 2-d model of Janes at ul. 1 1998) together with 
the corresponding mass streainiunctioti. Contour intervals are 0.2 DU/km. with zero contours omitted. Positive 
(negative) ozone anomaly values are lightly (darkly) shaded. Solid contours indicate positive strcamfunction. 




Figure 30 . Correlation of detrended TOMS monthly mean ozone column with the ozone column from the model 
of Kinnevsley and Fung j 1998] (forced by the observed Singapore wind) covering the period November 1978 April 
1993 (contour interval 0.1). 
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Figure 31. Height-latitude cross-sections of observed aerosol extinction ratio for two 40 day periods 
representative of the two phases of the QBO. (a) Westerly shear phase centered on 11 November 1984 (contour 
interval 2.5). (b) Easterly shear phase centered on 4 October 1988 (contour interval 0.5). From Trepte and Hitchman 
[19921. 



Figure 32. Height-time sections of interannual anoi 


The contour interval is 0.1 ppmv with zero contour omi 


e$ in HjO over the equator from the HALOE instrument 
I. From Randel et aL , [1998]. 
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Figure 34. Height - time cross-section of the QBO anomalies in NO, in the 10°S-10°N latitude band derived 
from HALOE observations. The contour interval is 0.2 ppbv, and zero contours omitted. These anomalies have been 
derived from regression analyses, following Rondel and Wu [1996], 
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*»gure 35. HOI ) HRDI measurements of the zonally averaged zonal eastward (westerly) winds in the tropical 
stratosphere and mesosphere (left) from 1992 to 1998. The lower figure shows the QBO from 20-40 km as 
descending easterly (green to blue) and westerly (red to yellow) winds. In the lower mesosphere (60-80 km) the wind 
structure ,s dominated by the semiannual oscillation (SAO). a twicc-yearlv reversal ,n the winds (BOTTOM) 

.nesosphtfre m km^M T™' “I'f (UPPei ' figUre) Sh ° WS lhat the of the QBO extends into the 

sphut (80 kni). Mesospheric w,nd changes occur coincident wiih the change of the QBO winds near 30 km 

This coupling between the mesosphere and the stratosphere is believed to be caused by small-scale upward , novine 
atmospheric waves as indicated by wiggly arrows. From the UARS brochure, modified from the original provided 
by M. Burrage and D. Ortland. Courtesy M. Schoeberl. g P 



Figure 36, Vertical distribution of the amplitude of the MQBO, MSAO, SSAO* SQBO and annual component at 
the equator. MQBO is based on the UARS/HRDI observation [Burrage et ai t 1996], SAO is based on the rocket 
observation at Ascension Island [Hirota, 1978] and the annual component is after CIRA - 1986. 
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Figure 37. Difference in 1000-hPa geopotential composites (meters) between westerly and easterly OBO 
composites. December-February monthly-mean NMC data for 1964-1996 are used. V 
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